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Abstract 



I derive a general effective theory for hot and/or dense quark matter. After introducing general 
projection operators for hard and soft quark and gluon degrees of freedom, I explicitly compute 
the functional integral for the hard quark and gluon modes in the QCD partition function. 
Upon appropriate choices for the projection operators one recovers various well-known effective 
theories such as the Hard Thermal Loop/ Hard Dense Loop Effective Theories as well as the 
High Density Effective Theory by Hong and Schafer. I then apply the effective theory to cold 
and dense quark matter and show how it can be utilized to simplify the weak-coupling solution 
of the color-superconducting gap equation. In general, one considers as relevant quark degrees 
of freedom those within a thin layer of width 2Aq around the Fermi surface and as relevant gluon 
degrees of freedom those with 3-momenta less than Agi. It turns out that it is necessary to choose 
Aq <C Agi, i.e., scattering of quarks along the Fermi surface is the dominant process. Moreover, 
this special choice of the two cut-off parameters Aq and Agi facilitates the power-counting of 
the numerous contributions in the gap-equation. In addition, it is demonstrated that both the 
energy and the momentum dependence of the gap function has to be treated self-consistently in 
order to determine the imaginary part of the gap function. For quarks close to the Fermi surface 
the imaginary part is calculated explicitly and shown to be of sub-subleading order in the gap 
equation. 



3 



4 



Contents 



1 Introductioril 

1.1 Quark matter and strong interactions! 



1.2 Color supcrconductivitvl 



1.3 Search for an effective approacbl 



2 A general effective action for quark matter 



2.1 Deriving the effective action 
"2.1.1 Setting the .stagel 



2.1.2 Integrating out irrelevant quark rnodej 

2.1.3 Integrating out hard gluon modes! . . . 



2.1.4 Tree-level effective action 



2.2 Recoverv of known effective theoriea 
" 2.2.1 HTT;/HDT: effective actinnl 



■2.2.2 High-densitv effective theory 
2.3 Towards a general effective Theori\ 



2.3.1 Power counting ouark loops at large and small T 



3 Application to color superconductivity 



3.1 Calculation of the QCD gao parameter 



3.1.1 CJT formalism for the effective theory! 



3.1.2 Dvson-Schwinger equations for relevant quarks and soft gluond 

3.1.3 Solution of the gap equation to sub-leading ordeil 



3.2 The imaginary part of the gap functioni 



3.2.1 Solving the complex ga.p eqiiatiol 



3.2.2 Estimating the order of magnitude of A and fil 



3.2.3 Estimating aud n\8] and 6A 



3.2.4 Reproducing Reo!>fei.- + ir^.k) to subleading order 



3.2.5 Calculating lm6(eA. + ir^.k) near the Fermi surface! 



4 Summary and Outlook! 



A Matsubara sums in quark loop j 



CONTENTS 



List of Figures 



1.1 The dependence of a, on the considered energy scale O 

1.2 The phase diagram of strongly interacting matter.! . . 



1.3 Feynman diagram for one gluon exchange among two quarks 

1.4 The separation of scales.l 



2.1 The full propagator for irreleyant auarksJ 



2.2 The diagrammatic symbol for the factor (1 + aAGnoo) 



-1 



2.3 The term -^ ,08-^, 



2.4 The graphical re presentation of the term Tr^ In Q22 ( ^-^^j) 



2.5 The term A^.Tr^ 



2.6 The fermionic contribution to the term Ao.7\ 



2.7 The term Ao.7m\ 

2.8 The term Ayll-^^A? according to Eq. 

2.9 The term /lon«/^ol 



2.10 The fermionic contribution to the term A^H 



2.11 The term /toH-,)/! 



2.12 The three- and four-gluon vertices in S a[A^ 



^1 



2.13 The term ^'i ^'i in the effectiye action 

2.14 The term Tr^ ln^22^[74i] in the effectiye action (f2.4H|) 

2.15 The term Tr^ In A^oH^i , ^1, ^'i] in the effe ctiye a ction ((TIH)) 



2.16 The term ^TrA^^^Tb in the effective action (12.461) 



2.17 A particular patch covering the Fermi surface.1 



2.18 The integration region in the quark-quarkhole contribution in Eq. (|2.99|) 



(EM 



3.1 Diagrammatic representation of F 

3.2 Diagrammatic representation of Fn 

3^3 The contour C i n Ea. (13^35]) .1 

3. 4 Same as in Fig. but for magnetic hard gluon exchange 



3.5 Evaluating the Matsubara sum for HDL-resummed gluon propagators 
3''^^^^ontourci^umventin^^l^^ .1 
3. 7 The contour C in Eq. 



3.8 Deforming the contour C 



3.9 



Hard gluon exchange with momentum p > A^vi ^ u. 
3.10 The integration regions of f and u>t in Eo. (|3.12'fj)T 



8 



LIST OF FIGURES 



4.1 Momentum regime of quarks near the Fermi surface and soft gluons 



117 



List of Tables 



3.1 Estimates f OT^A^^ and Ar<,,^ at different energy scales and C <^ M.I 



3.2 Estimates for and An„t at different energy scales and C M 



98 
98 



7U 



2J Estimates for A \^^ and A^..^.. at different energy scales and ( « M\ 101 



3.4 Estimates for .A^'*,, and -Ap^i. at different energy scales and ( < M .I 101 

^^^_^^imates f ar_B^^j_anA_Brn^^^ energy scales_aiid_^ _^_M.1 105 



3.6 Estimates for B„'^^ and 13n,,+ at different energy scales and C ^. M 



3.7 Estimates for B^^^ and Bp^u at different energy scales and C M 



105 
108 



M Estimates for and gp^io at different energy scales and ( < M.I 108 



9 



10 LIST OF TABLES 



Chapter 1 

Introduction 



1.1 Quark matter and strong interactions 

In contrast to philosophers who search for meaning in nature, physicists investigate the properties 
of nature and search for the laws that inanimate nature obeys. Apparently, the properties of 
any given piece of (known) matter depend on its temperature and its density. It is therefore 
not at all a naive question to ask: "What happens to some piece of matter if I heat and squeeze 
it. . . further and further?" One correct answer could be: "You may use QCD to describe it." 
At least this applies to "normal" matter made of atoms. Atoms have a nucleus, which is 
composed of neutrons and protons, which in turn consist of quarks. Quarks are fermions and 
interact with each other by exchanging gluons, which mediate the strong interaction. The charge 
corresponding to the strong interaction is called color and the quantum field theory describing 
this interaction Quantum Chromodynamics (QCD). Generally, all particles that interact strongly 
are composed of quarks and called hadrons. Those hadrons, which are composed of a quark and 
an antiquark, are called mesons and those composed of three quarks, as the neutrons and protons 
mentioned above, baryons. While all quarks carry color charge, all hadrons are in total color 
neutral. 

As long as the temperature and density of the considered system are not too large, i.e. 
below the scale Aqcd — 200 MeV, quarks are strongly coupled and always confined into these 
color neutral hadrons. As a consequence of this strong coupling perturbative approaches are 
impossible. The quarks are deconfined at temperatures or densities above Aqcd- Well above 
this scale the strong interaction exhibits the phenomenon of asymptotic freedom T, where QCD 
becomes weakly coupled. In this regime all hadrons vanish and quarks and gluons form a state 
called the quark-gluon plasma or simply quark matter [2]. 

Quarks also feel the electroweak and the gravitational force. The first one is described by a 
quantum field theory, which comprises electromagnetic and weak interactions. The theory for 
the latter, general relativity, describes gravitation in terms of the curvature of space-time and 
is a classical field theory. Due to their relative weakness as compared to the strong force it is 
often justified to neglect the gravitional and the electroweak against the strong interaction and 
describe the considered matter by QCD only. As a counter example related to this work, where 
QCD alone is not sufficient, one may allude to neutron stars. These are the compact remnants 
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of supernova explosions of type II OlSEj. They have extremely large masses comparable to the 
mass of our sun, but only very small radii of several kilometers. The density in their inner cores 
are possibly sufficient to deconfine the quarks. Since these stellar objects have to be electrically 
neutral and /9— equilibrated, electroweak processes have to be considered. Furthermore, neutron 
stars are held together gravitationally. In order to describe their bulk properties adequately, one 
therefore has to account for gravity as well. 

In the present work, however, only the strong interaction will be accounted for and, more- 
over, only its weak coupling regime will be considered. The aim is to derive an effective action 
for strongly interacting matter in the deconfined phase, i.e. for quark matter [6;. This is moti- 
vated by the occurence of several well seperated momentum scales at high temperatures and/or 
baryonic densities, cf. Sec. \1.'A\ The derivation is performed in Sec. 12.11 The field of application 
comprises low energy phenomena in weakly coupled quark matter. By construction this action 
will only contain low energy quark and gluon modes as explicit degrees of freedom. At large 
temperatures such modes are, e.g., long wave-length collective excitations [71IH1I1|. At large 
quark densities and low temperatures only those quarks, which are located close to the Fermi 
surface, are relevant degrees of freedom. The derived action can be adjusted to different high 
temperature and/or density regimes by a suitable choice of the projection operators, which sep- 
arate high energy from low energy quark and gluon modes. These projection operators contain 
two independent cutoff parameters, Aq for quarks and Agi for gluons. In Sec. l2.2l it is shown that 
well-known effective theories for strongly interacting quark matter can be assigned to special 
choices of these cutoff parameters in this general effective action. In chapter |21 the phenomenon 
of color superconductivity in quark matter will serve as a physical application to demonstrate 
the usefulness of the derived action in a semi-perturbative context. 

Before going into more details, however, it is first necessary to provide an introductory 
overview of the basic principles of QCD and its phase diagram. This will be done in the remainder 
of this section. In Sec. 11.21 the phenomenon of color superconductivity will be introduced. In 
Sec. ll.3l the general concepts of effective theories are explained and some caveats with respect to 
quark matter are pointed out. Finally, the specific features of the effective approach developed 
in this work are discussed. 

QCD: Asymptotic freedom, symmetries, and the phase diagram 

QCD is a non-Abelian gauge theory with the gauge group SU{Nc)c, where the index c refers 
to color and Nc is the number of quark colors. The quarks correspond to the fundamental 
representation of SU{Nc)c, whereas the gluons are the gauge bosons and belong to the adjoint 
representation of SU{Nc)c- Apart from the color index, 1 < c < Nc, quarks also carry a fiavor 
quantum number, 1 < / < Nf. Suppressing all quark indices for simplicity the Lagrange density 
for QCD reads [TOl HH IH US] 

Cqcb = {il^'D^ - m) V - ^ Fr F^, + C^, + £ghost ■ (l-l) 

Since quarks are fermions with spin 1 /2 their wave functions ij) are 4A'^c -^/-dimensional spinors. 
The Dirac conjugate spinor is defined by 'ip = V'^70) where are the Dirac matrices and m is 
the current quark mass, which is a diagonal matrix in fiavor space, m = rriij 5ij. The six known 
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flavors are called up (n), down (d), strange (s), charmed (c), bottom (6), and top (t). Their 
masses are ordered as niu — nid <C nig <^ rric <C rrih <C mt. Explicitly, the current quark masses 
appearing in vCqcd for the three lightest quark flavors are ~ 5 MeV, rud — 10 MeV, and 
rus — 100 MeV |14j . The masses for the three remaining quark flavors are so much larger that 
these flavors will not play any role in the following. 

The covariant derivative is defined as = 5^ — Ta, with the strong coupling constant 
g. The vector fields A'^ represent the gluons, where the adjoint color index a runs from 1 to 
N"^ — 1. The Nc X Nc matrices Ta are the generators of the local SU{Nc)c gauge symmetry. 
Throughout this work I use the representation Ta = Aa/2, where Aq are the Gell-Mann matrices. 
The gluonic field strength tensor is defined as 

Fr = O^^A: - d-'A^a + gfabc < A-^, , (1.2) 

where fabc are the structure constants of SU{Nc)c- They are defined through ifabcTc = [Ta,Th]. 
In the non-linear terms in Eq. (|1.2|) proportional to fabc the non-Abelian character of QCD 
manifests itself: the gluons carry color charge and couple to themselves. This ultimately gives 
rise to the numerous non-trivial features of this theory. The terms £gf and £ghost in Eq. (|1.1|) 
contain gauge fixing terms and the contribution of the Faddeev-Popov ghosts. Since they are of 
no further relevance for this work I will not specify them here. 

The QCD Lagrange density given in Eq. (|1.H) is needed to calculate the grand canonical 
partition function of QCD ^ 

Z{T,V,fi) = J Vi^VtpVA^ exp{5QCD[V5,V',^]} (1-3) 
with the tree-level action of QCD at quark chemical potential 



I (£ + /xAA) 
Jx 



:i.4) 



The quark chemical potential is associated with the net quark number conservation. The number 
density operator of the conserved net quark number is = tp'joip. From Eq. (|1.3|) one can derive 
thermodynamical quantities such as pressure, p{T,fi), entropy density, s{T,fj.), and particle 
number density, n{T,fi), 



piT,^^) = T ^ '^^ 



dV 



n(r,/x) = 



[1.5) 



In particular, the quark number density n{T, //) for massless quarks and at T <C /x is proportional 
to the third power of /i, n ~ /i^. Note that in the thermodynamical limit, V — > oo. In Z{T,fi) 
becomes an extensive quantity and as such is simply proportional to V. Therefore, the depen- 
dence of the pressure on the volume cancels out. In order to calculate the expectation value 
of any given operator O in the grand canonical ensemble one employs the following averaging 
prescription 

^^^i^T^Ai^ O exp{5QCD[V^,V',^]} • (1-6) 

Since QCD is a quantum field theory the vacuum is treated as a polarizable medium. Con- 
sequently, the result of any experiment that measures the effective color charge in the coupling 



14 



Chapter 1. Introduction 



constant = g"^ /^tt of two strongly interacting quarks will depend on the energy scale of the 
experiment. This dependence on the scale Q is expressed by the /?— function of QCD jl() | ll H [T^. 
In the so called MS-scheme at the three loop level one finds for the scale Q 



Q 



das 
dQ 



2 4 



where 



/3o 



11 



-iV 
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19 



N 
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9 ^ 27 



(1.7) 



:i. 



Since in nature the number of quark flavors is Nf = 6 it follows that /? < and QCD is indeed 
an asymptotically free theory as already mentioned above, cf. left diagram in Fig. 11.11 For a 
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Figure 1.1: The experimental and theoretical dependence of as{Q) on the considered energy scale Q = fi 
(left figure). Experimental values for as{Mz) (right figure), |14| . 

more quantitative analysis one has to solve the renormalization group equations. At three-loop 
order one finds ^S] 



asiQ) 



/5olog(Q2/A2^CD) 



2/3i log [log(QVA2jcD 
3! 



+ 



log log(QVA|cD) 



1 



log(QVA!3CD) /3o^log2(QVA2jcD) 
P2P0 5' 



+ 



:ii 



(1.10) 



This solution may be used to compare the values of as{Q) and Aqcd from different experiments. 
With Eq. H1.10() one may first calculate the scale Aqcd by matching it to the measured as{Q) 
and then use it to extrapolate Og to one standard reference scale, say to the mass of the Z 
boson, Mz = 91.1876 it 0.0021 GeV. This of course depends on the number of flavors used in 
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Eq. H1.1U() . For energies not much larger than Mz one may restrict oneself to Nf = 5 active 
flavors and finds for Aqcd the value jl4j 

A^^ = 217+^^ MeV . (1.11) 

MS -^■^ ^ ^ 

For as{Mz) one obtains Jl] 

a^(Mz) = 0.1187 ±0.002 , (1.12) 

cf. right diagram in Fig. 11.11 Decreasing the scale Q, the coupling as{Q) grows very large for 
Q Aqcd- Such a divergency of the coupling is sometimes called a Landau pole. In the case of 
quantum electrodynamics (QED) the Landau pole appears at very large energy scales, whereas 
the coupling is perturbatively small at low energies. 

Turning to the symmetries of the QCD Lagrangian given in Eq. Hl.l|l one observes that 
£qcd is invariant under Lorentz transformations and translations in space-time. Furthermore, 
it respects charge conjugation, parity and time-reversal invariance. Another symmetry is its 
invariance under gauge transformations, which are generally local, i.e. dependent of space-time, 
Uix) G SU{N,), 

ip{x) U{x)ip{x) , A^(x) ^ U{x)A^'{x)U^{x) + iU{x)d^'U^ (x) , (1.13) 

where A^^ = A'^Ta- A further symmetry is the global (i.e. independent on space-time) quark 
flavor symmetry of £qcd- To elaborate this, one may decompose the quark flelds with respect 
to their chirality 

1 ± 75 

tp = ^r + ^l , V'r/ = Vr/ ^ , Vr,i = ^ , (1.14) 

where Vr^e. are the chirality projectors. One observes that in £qcd only the mass term mixes 
quarks with different chirality 

■0* niij il)^ = -4)]. rriij ipl + i;} ruij ipi , (1-15) 

where use was made of the orthogonality of Vr and and of Vr,e7o = loPt.r- Hence, assuming 
zero quark masses for all flavors, m = 0, £qcd becomes chirally symmetric, i.e. invariant under 
global chiral U{Nf)r U{Nf)i transformations given by 



A,e Urf ll^rl , Url = CXp 



a'i.^.Ta G U{Nj)r,, . (1.16) 



V 



a=0 



Here, the are the parameters and Ta the generators of U{Nf)r/ with Tq ~ 1. While the 
kinetic and interaction parts of >Cqcd are generally invariant under arbitrary transformations 
of U{Nf)r U{Nf)i, the mass term is only invariant in the case that the parameters and 
af of the given transformations Ur and Ui are equal, a" = a^. One therefore introduces the 
vectorial group U{Nf)v C U{Nf)r (8) U{Nf)i with elements that fulfill Ur = Ui, which is often 
denoted as U{Nf)y = U{Nf)j.+i. These vectorial transformations leave >Cqcd invariant if all 
masses are equal. In nature, however, only the two lightest quarks have approximately equal 
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masses, which reduces the vector symmetry to an approxhnative U{2)v symmetry. Analogously 
one recovers an approximative U{3)v symmetry, if one considers energies much larger than nis. 
The subset U{Nf)A = U{Nf)/U{Nf)y comprises all transformations fulfilling Ur = Uj. Hence, 
the symmetry of these so-called axial transformations is always explicitly broken by the quark 
masses. 

Since any unitary transformation can be decomposed into the product of a special unitary 
transformation and a complex phase one may finally write U{Nf)r (8) U{Nf)i as SU{Nf)r (8) 
SU{Nf)i (g) U{l)v (8 U{1)a- The phase factors due to U{l)v correspond to an exact symmetry 
(independent of the quark masses) of QCD, corresponding to the conservation of baryonic charge, 
i.e. of quark numbers, and is often denoted as U{1)b- The U{1)a symmetry is explicitly broken 
by an anomaly at the quantum level of QCD jJTj. In hot and/or dense matter the instantons 
corresponding to this anomaly are screened JH] so that the U{1)a symmetry may become 
effectively restored again in matter. 

Below some critical density and temperature of the order of magnitude of Aqcd the true 
ground state of the vacuum is populated by the so called chiral condensate. In Fig. 11.21 the 
corresponding region is labelled by xSB. Often it is also referred to as the hadronic phase. It is 
described by the order parameter 



Its structure in flavor space connects right and left handed quarks. In the chiral limit of zero 
quark masses the difference among the quark flavors vanishes and = 6^^^, i.e. the order 
parameter of the chiral condensate becomes diagonal in flavor space just as the mass matrix 
discussed above. It follows analogosly that the chiral condensate spontaneously breaks the (ap- 
proximate) axial part of the chiral symmetry. Since the broken symmetry is global, Goldstone's 
theorem |1U1 llll [T^ applies and massless Goldstone bosons are created. Their number is given 
by the number of broken generators. Since in the considered case the chiral condensate breaks 
the SU{Nf)A symmetry, Nj-1 Goldst one bosons are generated. In the case of A'^^ = 2 (con- 
sidering only up and down quarks) this corresponds to the generation of the three pseudoscalar 
pions. These are not exactly massless, since the symmetry was already explicitly broken by the 
up and down quark masses and therefore it was only an approximate symmetry. Generally, one 
calls the massive Goldstone bosons corresponding to the breaking of approximate global sym- 
metries pseudo-Goldstone bosons. For Nf = 3 also the heavier strange quark is involved, which 
makes the spontaneously broken symmetry even more approximate as in the case of Nf = 2. 
Consequently, the created pseudo-Goldstone bosons are more massive. They correspond to the 
pseudoscalar meson octet made of pions, kaons, and the eta meson. Since all these mesons 
consist of a quark-antiquark pair, the quark number symmetry U{l)v and the gauge symmetry 
of electromagnetism ?7(l)em remain intact. 



Inside this hadronic phase the ground state of (infinite) nuclear matter with baryonic density 
no — 0.15 fm^'^ is located at T = and fj, = 308 MeV. The value of the chemical potential is 
given by the mass of the neutron, rriN = 939 MeV, reduced by the average binding energy per 
nucleon in nuclear matter, Bq = 16 MeV/ A One obtains /is = 939 MeV - 16 MeV = 923 




(1.17) 



XSB: SU{3)c0SU{3)v0U{1)b(^U{1) 



em 



(1.18) 
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Figure 1.2: The phase diagram of strongly interacting matter. 

MeV = 3^. The factor 3 arises from the fact that a nucleoli is a baryon with baryonic charge 1 
and as such consists of three (up and down) quarks each carrying baryonic charge 1/3. These 
quarks are dressed by their strong interaction with the chiral condensate, so that the quark 
masses become proportional to the expectation value of the chiral condensate. This effectively 
gives them much larger constituent masses as compared to the respective current quark masses 
m in >Cqcd- 

The line emerging from the point of infinite nuclear matter separates the gaseous from the 
liquid phase of hadronic matter and corresponds to a first order phase transition [2011^. At its 
endpoint at T ~ 10 MeV it becomes of second order exhibiting critical phenomena similarly to 
normal water HI HSl [21 Hi UHl ill • 

The evolution of strongly interacting matter in the early universe corresponds to the red 
line at zero net baryonic density, /i = 0, in Fig. 11.21 For small values of fj, the numerical 
lattice simulations of QCD [281 129j predict that the transition from the quark-gluon phase to 
the hadronic phase is a crossover [3()| I31j . These analyses fail so far in the regime ^ > T: 
lattice simulations are based on Monte Carlo importance sampling and rely on the probabilistic 
interpretation of the weight in the path integral. For fj, > the fermion determinant becomes 
complex, which causes the famous sign problem |281l29j . For small values of ^ a Taylor expansion 
technique has been proposed to investigate the critical endpoint of the chiral phase transition 
|32j . For recent experimental and theoretical progress on the properties of the confined and 
deconfined phases and the transitions among them see jSHllH^- As this field is not yet settled, 
the location of the line of the chiral phase transition and its critical endpoint in Fig. 11.21 are 
only approximate. 

Also the boundaries around the color superconducting regime ^USSHHEIEZI labelled by CSC 
are essentially unknown. Only in the limit of asymptotically large quark chemical potentials 
where QCD becomes weakly coupled it can be rigorously proven that the true ground state 
is color superconducting at sufficiently low temperatures. In this limit even the specific color 
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superconducting phase is known. It is the so-called CFL phase, which breaks again the chiral 
symmetry, cf. Sec. 11.21 At chemical potentials and densities of physical relevance, however, the 
the situation becomes less transparent. The densities inside neutron stars are bounded by the 
condition of hydrostatic stability to be n < 10 no [HZI- This corresponds to < 500 MeV, cf. the 
green line in Fig. 11.21 where QCD is not weakly coupled anymore. It is believed, however, that 
inside the cores of some neutron stars color superconductivity might occur and have observable 
effects on their properties |381 1391 HO] . The critical temperature and density for the onset of 
color superconductivity can only be estimated to be roughly of the order Hc ^ 400 and Tc < 50 
MeV. Furthermore, the search for the most dominant color superconducting phase for given T 
and fi is still ongoing wherefore the presumably complex structure inside the CSC region in Fig. 
11.21 is essentially unknown. These issues will be elaborated further in the following section. 

1.2 Color superconductivity 

The discovery of color superconductivity goes back to the late 1970's (33]. However, wider 
interest in the phenomenon of color superconductivity has only recently been generated by the 
observation that, within a simple Nambu-Jona-Lasinio (NJL) - type model ,41_ for the quark 
interaction, the so-called color-superconducting gap parameter (p assumes values of the order 
of 100 MeV 1321 • The gap parameter (j) appears in the dispersion relation of the gapped, i.e. 
Cooper paired quarks in the color superconducting medium 



where k is the modulus of the 3-momentum of the quark. The index r specifies the color and 
the flavor of the considered quark. In the special case that A^. = the quark is ungapped. It 
follows that gapped quarks carry a non-zero amount of 6 even on the Fermi 

surface, k = fi. Therefore, to excite (generate) a pair of gapped quarks costs at least the energy 
2\/A^</'- Such an excitation process can be interpreted as the break-up of a Cooper pair and the 
required energy as the Cooper pair binding energy [l^l |S1 EHl BHl HZI • 

Gap parameters of ~ 100 MeV would have important phenomenological consequences for 
the physics of neutron stars |381 1391 140j . It is therefore of paramount importance to put the 
estimates from NJL-type models |48| 149| 150| 1511 152j on solid ground and obtain a more reliable 
result for the magnitude of the gap parameter based on first principles. To this end, the color- 
superconducting gap parameter was also computed in QCD [521 [3^ ESI |5ni EZI • In this section a 
short review on color superconductivity is given including the mechanism of Cooper pairing in 
QCD, the different symmetry breaking patterns for various color superconducting phases, and 
some implications on the observables of neutron stars. At the end of this section a schematic 
version of the QCD gap equation in the weak coupling limit will be presented. It is suitable to 
identify and power-count the various terms that arise from the different gluon sectors. 

As already explained, at high densities the quarks are deconfined and free to move individ- 
ually. Since quarks are fermions, Pauli blocking forces them to build up a Fermi sea, where 
only quarks at the Fermi surface with Fermi momentum kp = — are able to interact 
and exchange momenta. For /i 3> m the typical momentum scale of these quarks is given by /U. 
Hence, for fj, ^> Aqcd the strong interaction becomes weakly coupled and single gluon exchange. 




(1.19) 
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cf. Fig. 11.31 becomes dominant. The quark-quark scattering amplitude for one-gluon exchange 
is proportional to ^| IH7| 

Z2 ^"^/i' = — {Sii'Sjj' - 5ij>5i>j) + {5ii>5jj> + 5ij>5i>j) . (1.20) 

a=l 

The indices i,j are the fundamental colors of the two quarks in the incoming channel, and 
the respective colors in the outgoing channel. If one interchanges two color indices in the incom- 
ing or in the outgoing channel the sign of the first term changes while the second term remains 
intact. The minus sign in front of the first, antisymmetric term indicates that this channel is 
attractive, whereas the other, symmetric term is repulsive. In terms of representation theory. 




Figure 1.3: Fcynman diagram for one gluon exchange among two quarks. 



Eq. ()1.2U|) corresponds to the coupling of two fundamental color triplets to the antisymmetric 
color antitriplet and the symmetric color sextet 

[3]c®[3]c = [3]?e[6]^ (1.21) 



Hereby, the 9-dimensional direct product of two 3-dimensional representations of SU{3)c has 
been uniquely decomposed into a direct sum of one 3-dimensional and one 6-dimensional repre- 
sentation. The term [3]" corresponds to the first term in Eq. (|1.2U|) and [6]^ to the second. 

The presence of an attraction among quarks is crucial. According to Cooper's theorem j43j 
any attractive interaction destabilizes the Fermi surface in favor of the formation of Cooper pairs 
[HI HSl 1^ |17j. Also for n ~ 400 — 500 MeV, where the theory is not yet in its weak coupling 
limit, ctsilJ') ~ 1, one expects the anti-triplet channel to be attractive. However, then it would 
be mediated by instantons 02 rather than by one-gluon exchange. 

As in the case of the chiral condensate, the Cooper pairs spontaneously reduce the symmetry 
of the system. Normal quark matter without any condensate is invariant under arbitrary 



NQ: SU{3)c0SU{Nf)v0SU{2)j0U{l)B0U{l)em (1.22) 

transformations, where also the global spin group SU{2)j is included with the total quark 
spin J = S + L. Nf is the number of quark flavors with approximately equal masses. The 
residual symmetry of the system in the presence of a Cooper pair condensate is given by the 
symmetry of the order parameter i.e. that of the Cooper pair wave function. Already from 
the symmetry breaking pattern of a specific color superconducting phase one can draw some 
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important conclusions for its key properties. In this context, however, it is necessary to remember 
the Anderson-Higgs mechanism I59j . It is roughly speaking the transfer of Goldstone's 
theorem onto gauge symmetries: if a gauge symmetry is spontaneously broken then those gauge 
bosons, which correspond to the broken generators, acquire masses. In this way they absorb 
the degrees of freedom, which would have been occupied by the Goldstone bosons if the broken 
symmetries had been global. One jargonizes this by saying that these gauge bosons "eat up 
the would-be-Goldstone-bosons and therefore become massive." (Strictly speaking, a gauge 
symmetry must first be fixed before it can be broken jHOj- After its fixing the gauge the symmetry 
has become global and the occurring breaking is explicit. The corresponding Goldstone bosons 
depend on the choosen gauge and are therefore non-physical. The masses of the gauge bosons, 
however, are independent of the gauge and are physical.) 

One canonical example for the Anderson-Higgs mechanism is the Weinberg-Salam model for 
the electroweak interactions [61_, which has a SU{2)j U{1)y gauge symmetry. The index / 
stands for the isospin and Y for the hypercharge of the left-handed leptons and quarks. One 
introduces the so-called Higgs field which spontaneously breaks the gauge symmetries to 
a residual C/(l)/+y. Physically one interprets this residual group as the gauge group U{l)em 
of electromagnetism with one massless gauge boson A, the photon. The three massive gauge 
bosons that emerge are identified with the and the Z boson of the weak interaction. Since 
the remaining gauge group U{1)i+y consists of joint rotations in isospin and hypercharge space, 
the Z and the A bosons must be generated by a "rotation" in the space of the gauge boson of the 
original U{1)y group with one of the gauge bosons of the original U{l)i group. Those "original" 
gauge bosons on the other hand become meaningless in the phase of the broken symmetries, for 
details see e.g. [101 CH US] . 

In the case of color superconductivity the Cooper pairs always carry color charge. Depending 
on the considered colorsperconducting phase, they therefore reduce the color gauge symmetry 
to either SU{2)c and give masses to five gluons, or break the color gauge symmetry completely. 
Then all eight gluons become massive. The respective gluon masses are called color Meissner 
masses and lead to the color Meissner effect: the respective magnetic gluon fields are expelled 
from the color superconducting medium. This is analogous to the more familiar electromag- 
netic Meissner effect, which occurs in electric superconductors. In that case the U{l)em gauge 
symmetry is spontaneously broken and the photons attain a Meissner mass. Since Cooper pairs 
of quarks are also electrically charged, a color superconductor may also be an electromagnetic 
superconductor depending on the breaking of the electromagnetic gauge symmetry in the con- 
sidered phase [72]. Analogously, the baryon number symmetry may be broken, since Cooper 
pairs have nonzero baryonic charge. A breaking of this symmetry is connected to the onset of 
superfluidity. 

In the following, the requirement of antisymmetry of the Cooper pair wave function will be 
used to analyze different symmetry breaking patterns. This approach, however, cannot give any 
information on the breaking of the U{1)b and C/(l)em symmetries. This as well as the values 
of the Meissner masses require more intensive calculations, which shall not be presented here. 
As shown above, the order parameter is an element of the representation [3]" of SU{3)c and is 
therefore antisymmetric in color space. In the case that the Cooper pair is a spin singlet, J = 0, 
the wave function must be an element of the antisymmetric representation [l]j of SU{2)j. If, 
however, it is a spin triplet, J = 1, the corresponding representation of SU{2)j is [3]j, which 
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is symmetric. Therefore, its representation with respect to SU{Nf)v must be antisymmetric if 
J = 0, and symmetric if J = 1. 

Starting with J = it follows that the Cooper pairs have to consist of two quarks with 
different flavors in order to fulfill antisymmetry. Assuming that only two flavors, i.e. up and 
down, participate in the pairing, Nf = 2, one finds as the corresponding antisymmetric 
representation, since analogosly to the spin one has in flavor space 

[2]v ® [2]v = m ® [3]y ■ (1.23) 

In color and flavor the most general ansatz for the order parameter reads 

^{f = e^Jke^' , (1.24) 

where the indices i,j are the color and /, g the flavor indices of the two quarks and the quantities 
eijk and e-^^ are totally antisymmetric tensors of rank 3 in color and and of rank 2 in flavor space, 
respectively. It follows that the order parameter is a (complex) 3- vector in the fundamental color 
space. By a global color rotation one may transform the order parameter $fe to point into the 
(anti-)3 direction in color space. Then only those quarks carrying (transformed) colors 1 and 
2 participate in the pairing. This demonstrates that in this phase the SU{3)c gauge symmetry 
is broken down to SU{2)c. It follows that five of the eight gluons acquire a Meissner mass in 
this phase. The remaining three gluons do not interact with the anti-triplett Cooper pairs and 
remain massless. One finally finds that for this phase, which is called the 2SC phase (indicating 
that only two quark flavors participate in the pairing), the residual symmetry is given by 

2SC: SUi2),<E)SU{2)v^Uil)B+c<S)U{l)em+c (1.25) 

Here, from the outset only the isospin symmetry SU{2)v was assumed, which is found unbroken. 
Furthermore, the condensate is invariant under joint rotations of U{1)b and U{l)em as well as 
under joint rotations of SU{2)c and U(l)em (EZlinnj. Since the latter two are both gauge groups 
this implies that, as in the Weinberg-Salam model, the photon is only rotated but remains 
massless (iO . Therefore the 2SC phase is no electromagnetic superconductor. The presence of 
the symmetry group U{1)b+c suggests that this phase exhibits nontrivial features concerning 
superfluidity. In helium-3 one encounters similar symmetry breaking patterns where the particle 
number conservation group is locked with rotations in spin or angular momentum space j47j . 

Assuming quark chemical potentials much larger than the strange quark mass one may 
consider the case that up, down, and strange quarks equally participate in the pairing. Then, 
for Nf = 3, it is 

[3]y® [3]y = [3]^e[6]f. , (1.26) 

wherefore the desired antisymmetric representation of SU{3)v is [3]j. Consequently, in Dirac 
and flavor space the ansatz for the condensate reads 

^{f = e^,kef"''^', . (1.27) 

It has a similar structure as the order parameter of superfluid helium-3 with spin 5 = 1 and 
angular momentum L = 1, which breaks the global SO{3)s and S0{3)l symmetries. There 
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one encounters a multitude of different pliases [ITj- Iii the present case of SU{3)c and SU{3)v 
symmetry breaking, however, it turns out that the following ansatz for <I>^ is dominant |H4| I65j 

^k = 4^- (1-28) 

This ansatz is similar to the order parameter of the chiral condensate, cf. Eq. H1.17() . in the chiral 
limit of zero quark masses, where the different flavors have identical properties and consequently 
^ij _ This order parameter simultaneously breaks the two global flavor symmetries 

SU{3)r SU{3)i to the one global vectorial symmetry SU{3)v- In the present case the global 
vectorial symmetry SU{3)v and the local symmetry SU{3)c are simultaneously broken to one 
global SU{3)c+v symmetry. Accordingly, the system is invariant under "locked" rotations in 
color and flavor space. This phase is therefore called the color-flavor-locked or CFL phase. All 
in all, it is invariant under 

CFL: SU{3)c+v(^SU{2)j(E)U{l)c+e^,n ■ (1.29) 

In contrast to the 2SC phase all three colors and flavors participate in the pairing and one 
therefore expects all gluons to interact with the Cooper pair condensate and attain a Meissner 
mass. This is in fact true, since the remaining color symmetry is only global. Furthermore, 
since C/(l)c+em is a gauge group the photon remains massless. However, it is mixed with one 
of the gluons. Therefore, the CFL phase is no electromagnetic superconductor, i.e. the Cooper 
pairs are neutral with respect to the rotated electric charge |66| . Due to the broken U{1)b 
it is a superfluid. In rotating systems as neutron stars superfluidity may result in vortices, 
through which magnetic fields can enter the superconducting phase in so called flux-tubes [HZ] • 
These vortices stick to the star's crust and and carry some amount of the stars total angular 
momentum. When the rotational velocity of the star gradually decreases due to electromagnetic 
radiation these vortices eventually rearrange and transfer some of the angular momentum to 
the star. As a consequence the rotational frequency is expected to suddenly increase leading to 
observable "glitches" [68] . 

As mentioned above, the CFL phase is physically significant only for /i ^ m^. For such 
large chemical potentials one could in principle neglect all masses and assume invariance under 
chiral symmetry rotations, SU{3)r '^SU{3)e instead of the vectorial SU{3)v- In the CFL phase 
left and right handed quarks separately form Cooper pair condensates, which corresponds to 
the breaking pattern SU{3)r+c "X" SU{3)e^c = SU{3)r-i-£+c- Therefore, also in the limit of very 
large densities the chiral symmetry is broken, however, not by the quark masses or the chiral 
condensate but by the CFL phase |69] . 

In the case that the Cooper pair carries spin 1, J = I, it may consist of two quarks of 
the same flavor [HEl EDI ESI ■ simplicity one first considers the case that there 

is only one quark flavor in the system, Nf = 1. As an element of [3]j the order parameter is 
a 3-vector in momentum space with components a = x,y, z. The most general ansatz for the 
order parameter reads in spin and color space [SSI ESI ESI ES] 

e^,fe<I>L (1.30) 

which generally breaks the SU{3)c ® SU{2)j symmetry and allows for a multitude of possible 
phases [731 174j . Here, only two special phases shall be analyzed, which are interesting due to 
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their analogy to the 2SC and the CFL phases of spin-0 color superconductors. These are the 
color-spin-locked or CSL phase 

$^ = 5^$, (1.31) 

and the polar phase, where 

$^ = 5^3^^"^^ • (1.32) 

In the first the color gauge group SU{2>)c and the global spin group SU{3)j are simultaneously 
broken to the global SU{3)c+j group, which "locks" color and spin rotations. One has in total 

CSL: SU{3)c+j (^U{l)v ■ (1.33) 

As in the CFL phase all gluons and the photon become massive. In contrast to the CFL 
phase, however, the U{1)b and the C/(l)em symmetry are broken. Hence, the CSL phase is 
a superfluid and a electromagnetic superconductor. The implications of the electromagnetic 
superconductivity on the observables of a neutron star are discussed in |72| . 

The polar phase is similar to the 2SC phase with SU {3)c'^SU (2) j breaking down to SU {2)c0 
U{l)j. One finds for the residual symmetry group 

polar: SU{2)^ <g) U{l)j U{l)v U{1)b+c ^ U{l)em+c ■ (1-34) 

Consequently, in the polar phase five gluons become massive. The condensate leaves rotational 
invariance in momentum space only with respect to one fixed axis. Indeed, the color supercon- 
ducting gap parameter, see below, exhibits only an axial symmetry [Z2|- The C/(l)em remains 
unbroken, wherefore the polar phase is no electromagnetic superconductor. However, a many 
flavor system, where each flavor separately forms a polar spin-1 gap, can be shown to be elec- 
tromagnetically superconducting |7Hj . 

The physical significance of one-flavor color superconductors is that in neutron stars the 
Fermi momenta of quarks with different flavors may become sufficiently different to rule out 
the normal BCS pairing mechanism, which requires exactly opposite quark momenta. Such a 
mismatch of Fermi surfaces indeed occurs in neutron stars if one respects its electrical neutrality 
and /3-equilibrium jSZl EZi EHl EOl EOl EH ESI • Total electric and color charge neutrality has to 
be fulfilled in order to allow for the gravitational binding of the star. Electrical charge neutrality 
for a two flavor system neglecting electrons, e.g., requires « 2n„, since the electrical charge 
of an up quark is Qu = 2/3, whereas Qd = —1/3 for a down quark. With n^j ^ ~ /^u d it 
therefore follows fi^ ~ 2^/^//„ ~ 1.26^u- In this simple case of two flavors, the requirement of 
/3— equilibrium will not change this estimate much, cf. |37j . 

In spin-1 color superconductors quarks of the same flavor may pair, wherefore the BCS- 
pairing mechanism remains always intact. Alternatively, other color superconducting phases are 
proposed with different pairing mechanisms. In the so-called interior gap or breached pairing 
the quarks with the smaller Fermi momentum are excited to a higher energy state in order to 
form Cooper pairs with the quarks with larger Fermi momentum |831 184j . Also displacements 
[551 ES] (LOFF phases) and deformations of the Fermi surfaces jHZ] are discussed to allow for 
Cooper pairing. The breached pairing mechanism has been shown to lead to Cooper pairs with 
gapless spectra [HHl ES EOl ED , which however appear to be unstable [HI ESI El ESI EHl • 
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To decide which of the many possible phases is favorable at given T and one has to 
determine the phase with the largest pressure. Generally, the pressure grows with the gain of 
condensation energy of the Cooper pairs. This in turn depends on the number of quark degrees 
of freedom, which participate in the Cooper pairing, and on the magnitude of the respective 
color superconducting gap parameters. The gap parameter is determined from the gap equation. 
For a general derivation see J3] and references therein. In Sec. 13. li the gap equation is rederived 
within the new effective formalism. In the following it is sufficient to concentrate on the gap 
equation in the schematic form 



Cln^ (^U/3 1n(^)+a 



(1.35) 



to discuss the various contributions from different gluon sectors and their respective orders of 
magnitude. To derive it one has to assume zero temperature, T = 0, weak coupling, g <^ 1, and 
apply the mean-field approximation. The solution is 

^ = 26// exp (^-0 [1 + 0(5)] . (1.36) 

The first term in Eq. ()1.35|) is of leading order since, according to Eq. ()1.36j) . g"^ \v? / (j)) ~ 1. 
It originates from the exchange of almost static, long-range. Landau-damped magnetic gluons 
j53j . One factor \n.{^/(j)) is the standard BCS logarithm, which arises when integrating over 
quasiparticle modes from the bottom to the surface of the Fermi sea, / dq/eq ~ ln(///(/)). (In 
fact, as will be shown in Sec. 13.2.'^ integrating over quasiparticle modes inside a layer around 
the Fermi surface of an "intermediate" width is sufficient to build up a BCS log.) The second 
factor comes from a coUinear enhancement ~ \n{^/eq) in the exchange of almost static 

magnetic gluons. The coefficient C, determines the constant c in the exponent in Eq. p.36|) . As 
was first shown by Son |53j . 

c.^. ,1.37) 

The second term in Eq. (|1.35|) is of subleading order, g^ln(/x/0) ^ g <^ 1. It originates from 
two sources. The first is the exchange of electric and non-static magnetic gluons [HHI^ I ^ I^. 
In this case, the single factor \n{fi/(j)) is the standard BCS logarithm. The second source is 
the quark wave-function renormalization factor in dense quark matter j97| l9Sj. Here, the BCS 
logarithm does not arise, but the wave- function renormalization contains an additional ln(/i/eq) 
which generates a ln(fi/(p). The coefficient /5 determines the prefactor b of the exponent in Eq. 
H1.36() . For a two-flavor color superconductor, 

where Nj is the number of (massless) quark flavors participating in screening the gluon exchange 
interaction. The third term in Eq. (|1.35j) is of sub- subleading order, ~ g^. The coefficient a 
determines the 0{g) correction to the prefactor of the color-superconducting gap parameter in 
Eq. (|1.36|) . Since a has not yet been determined, the gap parameter can be reliably computed 
only in weak coupling, i.e., when the 0{g) corrections to the prefactor are small. 
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Weak coupling, g <^ 1, however, requires asymptotically large quark chemical potentials, 
fi ^ Aqcd- The range of /U values of phenomenological importance is, however, < 1 GeV. 
Although the quark density n is already quite large at such values of /i, n ~ 10 times the nuclear 
matter ground state density, the coupling constant is still not very small, 5 ~ 1. It is therefore 
of interest to determine the coefficient of g in the 0(g) corrections to the prefactor in Eq. H1.36() . 
If it turns out to be small, one gains more confidence in the extrapolation of the weak-coupling 
result to chemical potentials of order ~ 1 GeV. 

It is worthwhile mentioning that an extrapolation of the weak-coupling result (|1.36)) for a 
two-flavor color superconductor, neglecting sub-subleading terms altogether and assuming the 
standard running of g with the chemical potential fi, yields values of (p of the order of ~ 10 MeV 
at chemical potentials of order ~ 1 GeV, cf. Ref. |1^. This is within one order of magnitude of 
the predictions based on NJL-type models and thus might lead one to conjecture that the true 
value of (j) will lie somewhere in the range ~ 10 — 100 MeV. 

However, in order to confirm this and to obtain a more reliable estimate of </> at values of 
of relevance in nature, one ultimately has to compute all terms contributing to sub-subleading 
order. In the following section the approach to this problem by an effective theory will be 
motivated. 



1.3 Search for an effective approach 

Although possible in principle, the task to calculate the sub-subleading order contributions to 
the color superconducting gap parameter is prohibitively difficult within the standard solution of 
the QCD gap equation in weak coupling. So far, in the course of this solution terms contributing 
at leading and subleading order have been identified, cf. discussion after Eq. (|1.36|) . However, 
up to date it remained unclear which terms one would have to keep at sub-subleading order. 
Moreover, additional contributions could in principle arise at any order from diagrams neglected 
in the mean-field approximation |97| I99j. Therefore, it would be ideal to have a computational 
scheme, which allows one to determine a priori, i.e., at the outset of the calculation, which 
terms contribute to the gap equation at a given order. 

As a first step towards this goal, note that there are several scales in the problem. Besides 
the chemical potential n, there is the inverse gluon screening length, which is of the order of the 
gluon mass parameter rUg. At zero temperature and for A'^^ massless quark flavors [0], 

™' = ^/§?' (1-39) 

i.e., rUg ~ gfi. Finally, there is the color-superconducting gap parameter cp, cf. Eq. H1.36() . In 
weak coupling, g <C 1, these three scales are naturally ordered, 4> <^ gn <^ js. This ordering 
of scales implies that the modes near the Fermi surface, which participate in the formation of 
Cooper pairs and are therefore of primary relevance in the gap equation, can be considered to 
be independent of the detailed dynamics of the modes deep in the Fermi sea. This suggests that 
the most efficient way to compute properties such as the color-superconducting gap parameter 
is via an effective theory for quark modes near the Fermi surface. Such an effective theory has 
been originally proposed by Hong |100| llOlj and was subsequently refined by others |102| I103| 
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At this point it is worthwhile reviewing the standard approach to derive an effective theory 
jl()7|ll()^ll()9j . In the most simple case, one has a single scalar field, (p, and a single momentum 
scale. A, which separates relevant modes, ip, from irrelevant modes, ip, (j) = ip + ip. The relevant 
modes live on spatial scales L ^ 1/A, while the irrelevant modes live on scales / < 1/A <C cf. 
Fig. 11.41 In the derivation of the effective action, one is supposed to integrate out the microscopic, 

"relevant" modes — separation ^ "irrelevant" modes 
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Figure 1.4: The separation of scales. 

irrelevant modes. Usually, however, this is not done explicitly. Instead, one constructs all 
possible operators Oi composed of powers of the field and its derivatives, which are consistent 
with the symmetries of the underlying theory, and writes the effective action as |1U8| 

Se{i[^]= I Y.9^0^{^). (1.40) 
JX ■ 

The coefficients, or vertices, gi determine the interactions of the relevant modes ip. A priori, they 
are unknown functions of the single scale A, gi = gi{A). All information about the microscopic 
scale / is contained in these vertices. Since the microscopic scale / <^ L, the operators Oi are 
assumed to be local on the scale L. 

The effective action (|1.4()|) contains infinitely many terms. In order to calculate physical 
observables within the effective theory, one has to truncate the expansion after a finite number 
of terms. One can determine the order of magnitude of various terms in the expansion (|1.4fl|) via 
a dimensional scaling analysis, which allows to classify the operators as relevant (they become 
increasingly more important as the scale L increases), marginal (they do not change under 
scale transformations), and irrelevant (they become increasingly less important as the scale L 
increases). To this end, one determines the naive scaling dimension of the fields, dim((/j) = 5, 
from the free term in the effective action. Then, if the operator Oi consists of M fields ip and 
derivatives, its scaling dimension is dim(C'j) = 6i = M6 + N. The operator Oi is then of order 
~ L~^\ For dimensional reasons the constant coefficients gi must then be of order ~ A'^"''*, 
where d denotes the dimensionality of space-time. Including the integration over space-time, the 
terms in the expansion (|1.4()j) are then of order ~ (LA)'^~^\ Consequently, relevant operators 
must have 5i < d, marginal operators 6i = d, and irrelevant operators 6i > d. At a given 
scale L, one has to take into account only relevant, or relevant and marginal, or all three types 
of operators, depending on the desired accuracy of the calculation. To calculate a process at 
energy E = sA with s ^ 1 and error of order s"'*'^ one only needs operators Oi with dimension 
Si < D + n. The final result still depends on A through the coefficients gi{A). This dependence 
is eliminated by computing a physical observable in the effective theory and in the underlying 
microscopic theory, and matching the result at the scale A. 
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There are, however, cases where this naive dimensional scahng analysis fails to identify the 
correct order of magnitude, and thus the relevance, of terms contributing to the effective action. 
In the following, three examples shall be discussed briefly. For the first example, consider 
effective theories where, in contrast to the above assumption, the vertices gi are in fact non-local 
functions. Such theories are, for instance, given by the "Hard Thermal Loop" (HTL) or "Hard 
Dense Loop" (HDL) effective actions (3 El- In these effective theories, valid at length scales 
L ~ l/{gT) or ~ l/{gfj,), respectively, there are terms gn^^ in the effective action, which are 
constructed from a quark or gluon (or ghost) loop with n external gluon legs; A is the external 
gluon field with 5 = 1. The coefficients gn are non-local and do not only depend on the scale 
A < T, or < /U, but also on the relevant momentum scale 1/L ~ gT, or ~ g^. Naively, one 
would expect gn to belong to a local n-gluon operator and to scale like A'^~". Instead, it scales 
like L"'"^ [?]• For arbitrary n, the corresponding term gnA^ in the effective action then scales 
like L^, independent of the number n of external gluon legs. 

The second example pertains to the situation when there is more than one single momentum 
scale A. As explained above, for a single scale A and a given length scale L, the naive dimensional 
scaling analysis unambiguously determines the order of magnitude of the terms in the expansion 
(|1.4()|) . Now suppose that there are two scales, Ai and A2. Then, the vertices gi may no longer 
be functions of a single scale, say Ai, but could also depend on the ratio of A2/A1. Two scenarios 
are possible: (a) two terms in the expansion (|1.4()|) . say gnOn and gmOm-, with the same scaling 
behavior may still be of a different order of magnitude, or (b) the two terms can have a different 
scaling behavior, but may still be of the same order of magnitude. In case (a), all that is required 
is that the operators On and Om scale in the same manner, say L^^, and that g^ ~ Ag"'^, but 
gm ~ . If Ai ^ A2, gm ^ 9ni and thus the two terms are of different order of magnitude. 
In case (b), assume l/L <C Ai <C A2, with A1/A2 ~ l/(AiL) ~ e <C 1 and let take the fields 
to have naive scaling dimension 5 = 1. Then, at a given length scale L, a term gn^"", with a 
coefficient gn of order A2~", can be of the same order of magnitude as a term gm^"^, rn ^ n, 
if the coefficient gm ~ A^~'"(A2/Ai)'^ with k = d + m — 2n. Although the scaling behavior of 
the two terms is quite different as L increases, they can be of the same order of magnitude, if 
the interesting scale L happens to be ~ A2/AI. In both cases (a) and (b) the naive dimensional 
scaling analysis fails to correctly sort the operators Oi with respect to their order of magnitude. 

The third example where the naive dimensional scaling analysis fails concerns quantities, 
which have to be calculated self-consistently. Such a quantity is, for instance, the color- 
superconducting gap parameter, which is computed from a Dyson-Schwinger equation within a 
given many-body approximation scheme. In this case, the self-consistent solution scheme leads 
to large logarithms, like the BCS logarithm in Eq. (|1.35|) . These logarithms cannot be identified 
a priori on the level of the effective action, but only emerge in the course of the calculation |55j . 

In order to avoid these failures of the standard approach, in this work I pursue a different 
venue to construct an effective theory. I introduce cut-offs in momentum space for quarks, Aq, 
and gluons, Agi. These cut-offs separate relevant from irrelevant quark modes and soft from 
hard gluon modes. I then explicitly integrate out irrelevant quark and hard gluon modes and 
derive a general effective action for hot and/or dense quark-gluon matter. One advantage of this 
approach is that I do not have to guess the form of the possible operators Oi consistent with 
the symmetries of the underlying theory. Instead, they are exactly derived from first principles. 
Simultaneously, the vertices gi are no longer unknown, but are completely determined. Moreover, 
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in this way I construct all possible operators and thus do not run into the danger of missing a 
potentially important one. 

I shall show that the standard HTL and HDL effective actions are contained in the general 
effective action for a certain choice of the quark and gluon cut-offs Aq, Agi. Therefore, the 
new approach naturally generates non-local terms in the effective action, including their correct 
scaling behavior, which, as mentioned above, does not follow the rules of the naive dimensional 
scaling analysis. I also show that the action of the high-density effective theory derived by Hong 
and others |inn| llOU 11021 11031 11041 llOBj is a special case of the general effective action. In this 
case, relevant quark modes are located within a layer of width 2Aq around the Fermi surface. 

The two cut-offs, Aq and Agi, introduced in the present approach are in principle different, 
Aq 7^ Agi. The situation is then as in the second example mentioned above, where the naive 
dimensional scaling analysis fails to unambiguously estimate the order of magnitude of the 
various terms in the effective action. Within the present approach, this problem does not occur, 
since all terms, which may occur in the effective action, are automatically generated and can be 
explicitly kept in the further consideration. I shall show that in order to produce the correct 
result for the color-superconducting gap parameter to subleading order in weak coupling, one 
has to demand Aq < gfi <^ Agi < /i, so that Aq/Agi ^ g <^ 1. Only in this case, the dominant 
contribution to the QCD gap equation arises from almost static magnetic gluon exchange, while 
subleading contributions are due to electric and non-static magnetic gluon exchange. The fact 
that Aq <C Agi is not entirely unexpected: at asymptotically large densities, where the scale 
hierarchy is </> <C ^/i ~ mg <^ fi, the dominant contribution in the QCD gap equation arises 
from gluons with momenta of order ~ (m^0)^/^ [531 1541 while typical quark momenta lie in 
a shell of thickness ~ 2(/> <C {mg(f>y^^ around the Fermi surface. 

The color-superconducting gap parameter is computed from a Dyson- Schwinger equation for 
the quark propagator. In general, this equation corresponds to a self-consistent resummation of 
all one-particle irreducible (IPI) diagrams for the quark self-energy. A particularly convenient 
way to derive Dyson-Schwinger equations is via the Cornwall-Jackiw-Tomboulis (CJT) formal- 
ism |106j . In this formalism, one constructs the set of all two-particle irreducible (2PI) vacuum 
diagrams from the vertices of a given tree-level action. The functional derivative of this set with 
respect to the full propagator then defines the IPI self-energy entering the Dyson-Schwinger 
equation. Since it is technically not feasible to include all possible diagrams, and thus to solve 
the Dyson-Schwinger equation exactly, one has to resort to a many-body approximation scheme, 
which takes into account only particular classes of diagrams. The advantage of the CJT for- 
malism is that such an approximation scheme is simply defined by a truncation of the set of 
2PI diagrams. However, in principle there is no parameter, which controls the accuracy of this 
truncation procedure. 

The standard QCD gap equation in mean-field approximation studied in Refs. [Hll ESI [SE| 
follows from this approach by including just the sunset- type diagram, which is constructed from 
two quark-gluon vertices of the QCD tree-level action (see, for instance. Fig. 13.11 in Sec. 13.1.11) . 
I also employ the CJT formalism to derive the gap equation for the color-superconducting gap 
parameter. However, I construct all diagrams of sunset topology from the vertices of the general 
effective action derived in this work. The resulting gap equation is equivalent to the gap equation 
in QCD, and the result for the gap parameter to subleading order in weak coupling is identical 
to that in QCD, provided Aq < (//i <C Agi < fj,. The advantage of using the effective theory is 
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that the appearance of the two scales Aq and Agi considerably facilitates the power counting 
of various contributions to the gap equation as compared to full QCD. I explicitly demonstrate 
this in the course of the calculation and suggest that, within this approach, it should be possible 
to identify the terms, which contribute beyond subleading order to the gap equation. Of course, 
for a complete sub-subleading order result one cannot restrict oneself to the sunset diagram, 
but would have to investigate other 2PI diagrams as well. This again shows that an a priori 
estimate of the relevance of different contributions on the level of the effective action does not 
appear to be feasible for quantities, which have to be computed self-consistently. 

This work is organized as follows. In Sec. 12. Il l derive the general effective action by explicitly 
integrating out irrelevant quark and hard gluon modes. In Sec. 12.21 1 show that the well-known 
HTL/HDL effective action, as well as the high-density effective theory proposed by Hong and 
others, are special cases of this general effective action for particular choices of the quark and 
gluon cut-offs Aq and Agi, respectively. Furthermore, I power count one special subset of dia- 
grams, loops of irrelevant quark modes, analogously to the HTL power counting scheme j7j for 
the choice (/> <C Aq <C <C Lgi < fi. 

Section 13.11 contains the application of the general effective action to the computation of 
the real part of the color-superconducting gap parameter. In Sec. 13.21 1 consider the energy 
dependence of the gap function within the effective theory and compute its imaginary part for 
quarks near the Fermi surface. In Sec.l^I conclude this work with a summary of the results and 
an outlook. 

The units I use are h = c = ks = 1- 4- vectors are denoted by capital letters, = 
(fco,k), with k being a 3-vector of modulus |k| = k and direction k = k/A;. For the sum- 
mation over Lorentz indices, I use a notation familiar from Minkowski space, with metric 
gfj.u _ (Jiag^_|_^ _^ _^ although I exclusively work in compact Euclidean space-time with vol- 
ume V/T, where V is the 3-volume and T the temperature of the system. Space-time integrals 
are denoted as Jq^^ dr /yd^x = J-^. Since space-time is compact, energy-momentum space is 
discretized, with sums {T/V)J2k = Z^k- -^o^ ^ large 3-volume V, the sum over 

3-momenta can be approximated by an integral, {^/y)J2k — / '^^k/(27r)'^. For bosons, the sum 
over n runs over the bosonic Matsubara frequencies co^l = 2n7rT, while for fermions, it runs over 
the fermionic Matsubara frequencies = (2n + 1)itT. In a Minkowski-like notation for four- 
vectors, xq = t = —IT, kf) = —iujn^^. The 4-dimensional delta-function is conveniently defined 
as = ,5(r)(5(3)(x) = -i dix'^) 6^^\^). 
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Chapter 2 



A general effective action for quark 
matter 



2.1 Deriving the effective action 

In this section, I derive a general effective action for hot and/or dense quark matter. I start 
from the QCD partition function in the functional integral representation (Sec. 12.1.1)1 . I first 
integrate out irrelevant fermion degrees of freedom (Sec. I2.1.2() and then hard gluon degrees of 
freedom (Sec. I2.1.3j) . The final result is Eq. ()2.46|) in Sec. 12.1.41 I remark that the same result 
could have been obtained by first integrating out hard gluon modes and then irrelevant fermion 
modes, but the intermediate steps leading to the final result are less transparent. 

2.1.1 Setting the stage 

As discussed in the introduction, cf. Eq. ()l.llll.4j) . the partition function for QCD in the absence 
of external sources may be written as 




(2.1) 



where for convenience the pure (gauge-fixed) gluon action has been introduced 



Sa[A] = [-iFr(X)F;,(X)] +Sghost[^] 



(2.2) 



The partition function for quarks in the presence of gluon fields is 




(2.3) 



where the quark action is 




(2.4) 
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In fermionic systems at nonzero density, it is advantageous to additionally introduce charge- 
conjugate fermionic degrees of freedom, 

^c(^) = C^^(X) , ^c{X) = C , ^{X) = Ci>l{X) , i^{X) = il^liX) c , 

(2.5) 

where C = ^7^70 is the charge-conjugation matrix, C ^ = = C'^ = —C, C ^7jC = — 7^^; a 
superscript T denotes transposition. I may then rewrite the quark action in the form 

S,[A, = 1 ^ ^ g^HX, Y) ^{Y) + I ^ ^{X) t^A-iX) ^{X) , (2.6) 
where I defined the Nambu-Gor'kov quark spinors 

) , ^ = {A^c) , (2.7) 
and the free inverse quark propagator in the Nambu-Gor'kov basis 

e-,x,r).(l«Jl-f-^) , (2.8) 

with the free inverse propagator for quarks and charge-conjugate quarks 

[Gi]-\X, Y) ^ {ip^ ± ^70 - m) 5W(X - Y) . (2.9) 
The quark-gluon vertex in the Nambu-Gor'kov basis is defined as 

As I shall derive the effective action in momentum space, I Fourier-transform all fields, as well 
as the free inverse quark propagator, 

*(X) = -^^e^^-^*(i^), (2.11b) 



K 
K,Q 

Am = ^j:e-^^-'' Am ■ (2.11d) 



p 

The normalization factors are chosen such that the Fourier-transformed fields are dimensionless 
quantities. The Fourier-transformed free inverse quark propagator is diagonal in momentum 
space, too, 

G-yKO)--( ° ]6^'^ (212) 
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where [G^]-^{K) = If ± fi-fo - m. 

Due to the relations (|2.5j) . the Fourier-transformed charge-conjugate quark fields are re- 
lated to the original fields via ijjciK) = Cip^{—K), ipciK) = ilJ^{—K)C . The measure of the 
functional integration over quark fields can then be rewritten in the form 

Vi)Vi} = YldipiK) dipiK) = M Y[ di!{K)d7p{K)d^p{-K)d7p{-K) 

= Yl d^{K)dip{K)d^c{K)dtl;c{K) =J\f" H d-^{K)d'^{K) 

(K-K) {K-K) 

= V^V^ , (2.13) 

with the constant normalization factors M, A/"', M" . The last identity has to be considered as a 
definition for the expression on the right-hand side. 

Inserting Eqs. (|2.11|) - (|2.13|) into Eq. (|2.,''{j) . the partition function for quarks becomes 

Zq[A] = J P^P^exp 



Here, I employ a compact matrix notation, 



(2.14) 



^ [Qo' + gA) ^^Y. [Qo\K,Q)+gA{K,Q)\ ^{Q) , (2.15) 

K,Q 

with the definition 

AK, Q) ^ n^liK - Q) . (2.16) 

The next step is to integrate out irrelevant quark modes. 

2.1.2 Integrating out irrelevant quark modes 

Since I work in a finite volume V , the 3-momentum k is discretized. Let us for the moment 
also assume that there is an ultraviolet cut-off (such as in a lattice regularization) on the 3- 
momentum, i.e., the space of modes labelled by 3-momentum has dimension D < oo. I define 
projection operators Vi, V2 for relevant and irrelevant quark modes, respectively, 

^-1 = Pi ^ , ^'2 = ^2^, '^'i = ^70^170 , ^2 = ^70^270 • (2.17) 

The subspace of relevant quark modes has dimension A^i in the space of 3-momentum modes, 
the one for irrelevant modes dimension A'^2) with Ni + N2 = D. 

At this point, it is instructive to give an explicit example for the projectors Pi, 2- In the 
effective theory for cold, dense quark matter, which contains the high-density effective theory 
Unni Unil Uni CnSl Uni IMI discussed in Sec. rrn\ as special case and which I shah apply in 
Sec. 13. II to the computation of the gap parameter, the projectors are chosen as 

ViiK,Q) ^ (^^^ e{A^-\k-kF\)d%, (2.18a) 

= I, A+ + A^e(|fe-fc^|-Aq) J ^^.Q- 

(2.18b) 
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Here, 

Ak = ^ [^k + 670 (7 • k + m)] , (2.19) 

are projection operators onto states with positive (e = +) or negative (e = — ) energy, where 
= \/k^~+~m? is the relativistic single-particle energy. The momentum cut-off Aq controls 
how many quark modes (with positive energy) are integrated out. Thus, all quark modes within 
a layer of width 2Aq around the Fermi surface are considered as relevant, while all antiquark 
modes and quark modes outside this layer are considered as irrelevant. Note that, for the Nambu- 
Gor'kov components corresponding to charge-conjugate particles, the role of the projectors onto 
positive and negative energy states is reversed with respect to the Nambu-Gor'kov components 
corresponding to particles. The reason is that, loosely speaking, a particle is actually a charge- 
conjugate antiparticle. For a more rigorous proof compute, for instance, ipciiK) = C (— X) 
using '0i(-K) = ^(-ivT) 70 A+k7o (for \k - < Aq) and 7o C[A+J^C7-i 70 = A,^. In Sec. lO 
I shall discuss other choices for the projectors "Pi, 2) pertaining to other effective theories of hot 
and/or dense quark matter. The following discussion in this section, however, will be completely 
general and is not restricted to any particular choice for these projectors. 
Employing Eq. 1)2. 17p . the partition function 1)2. 14|) becomes 

^M]= f n ^*n2^*nexp|i E ^nGn^^m] ■ (2.20) 
n=l,2 \ n,m=l,2 / 

From now on, ^'1,2, ^1,2 are considered as vectors restricted to the A^i^2-dimensional subspace 
of relevant /irrelevant 3-momentum modes. The matrices G^n^ ^ = 1)2, are defined as 

g-^iK, Q) = GoU^, Q) + gAnniK, Q) , (2.21) 

where the indices indicate that, for a given pair of quark energies fco, qo, the 3- momenta k, q 
belong to the subspace of relevant (n = 1) or irrelevant (n = 2) quark modes, i.e., is an 
{Nn X A'„)-dimensional matrix in 3-momentum space. The matrices Q^m-: n ^ m, reduce to 

g-i{K,Q)=gAnm{K,Q) , (2.22) 

since Qq^ is diagonal in 3-momentum space, i.e. Q^nrn = for n 7^ m. For a given pair of quark 
energies /cqi Qq-, Qnm. is a (A''„ x A^m)-dimensional matrix in 3-momentum space. 

The Grassmann integration over the irrelevant quark fields ^'2 can be done exactly, if 
one redefines them such that the mixed terms ~ Q^mi n ^ m, are eliminated. To this end, 
substitute 

T = ^-2 + ^22 Q21 ^1 , T = ^2 + ^i G12 S22 , (2.23) 
where ^22 is the inverse of ^22^, defined on the subspace of irrelevant quark modes. The result 



'22 

is 



Zg[A] = J V^iV^i exp 



i (^lY - Gu' G22 Q21) ^1 + ^ Tr, In ^22' 



(2.24) 



The trace in the last term runs over all irrelevant quark momenta K, and not only over pairs 
(K,—K), as prescribed by the integration measure, Eq. 1)2. 13|) . This requires an additional 
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factor 1/2 in front of the trace. A more intuitive way of saying this is that this factor accounts 
for the doubhng of the quark degrees of freedom in the Nambu-Gor'kov basis. Of course, the 
trace runs not only over 4-momenta, but also over other quark indices, such as Nambu-Gor'kov, 
fundamental color, flavor, and Dirac indices. I indicated this by the subscript "g". 
For a diagrammatic interpretation it is advantageous to rewrite 

Gn - Qu G22 Q21 = G^li + 9B , (2.25) 

where 

gB = gAii - gAi2 Q22 9A21 ■ (2.26) 

The propagator for irrelevant quark modes, ^22; has an expansion in powers of (7 times the gluon 
field, 

00 

Q22 = Go,22 E [A22 Qo,22r ■ (2.27) 

n=0 

This expansion is graphically depicted in Fig. 12.11 




Figure 2.1: The full propagator for irrelevant quarks. The right-hand side symbolizes the expansion 
(|2.27|l . The free irrelevant quark propagators Go,22 are denoted by double lines, the gluon fields A22 by 
curly lines. 



I = 1 + + ^ i + 

Figure 2.2: The diagrammatic symbol for the factor (1 + gAGo,22)^^ ■ 




Figure 2.3: The term VPi gB'i'i. A relevant quark field is denoted by a single solid line. 

Using this expansion, and suppressing the indices on A, Eq. ()2.26() can be symbolically 
written as 

gB = (l + gAgo,22r' gA , (2.28) 

which suggests the interpretation of the field i3 as a "modified" (non-local) gluon field. In the 
diagrams to be discussed below, the factor (1 -|- gAgo,22)~^ will be denoted by the diagrammat- 
ical symbol shown in Fig. 12.21 With this symbol, the expression ^1 gB can be graphically 
depicted as shown in Fig. 12.31 
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= 0^-0^-0-^-0:- 

Figure 2.4: The graphical representation of the term Trqln(722^ in Eq. H2.24II . 
Since 

In ^22' = In Go,l2 - E 9" [^0,22 ^22]" , (2.29) 

n=l 

the last term in the exponent in Eq. H2.24() also has a graphical interpretation, shown in Fig. 

m 

This concludes the integration over irrelevant quark modes. Note that the present treatment 
is (i) exact in the sense that no approximations have been made and (ii) completely general, 
since it is independent of the specific choice ()2.18() for the projection operators. The next step 
is to integrate out hard gluon modes. 



2.1.3 Integrating out hard gluon modes 



Combining Eqs. (|2.1|) . (|2.24j) . and (|2.26|) . the partition function of QCD for relevant quark modes 
and gluons reads 



(2.30a) 



S[A,^^,^,] EE SA[A] + ^^i{g^},+gB[A]}-i'i + ^Trglng^^^[A], (2.30b) 

where T>A = YlpdA{P). For the sake of clarity, I restored the functional dependence of the 
"modified" gluon field B and the inverse irrelevant quark propagator Q22 on the gluon field A. 
The gluon action in momentum space is 



Sa[A] = -i 5: A;1(Pi)[Ao-i][^^Pi,P2)A^(P2) 



Pl,P2 

1 9 
3! Wf^ 

1 f 9 



E ^pIpm,o K%{Pi,P2, P3) A2{Pi)A^,{P2)A2{Ps) 



Pl,P2,P3 
2 



4! VVfT3 



E ^K+PMP,,, KIs AZ{PM^{P2)A2{PM'd{PA) 

Pl,-,P4 



(2.31) 



Here, Aq^{Pi, P2) is the gauge-fixed inverse free gluon propagator. To be specific, in general 
Coulomb gauge it reads 



nu 1 



(4) 



-P2 ' 



A 



-1 



(P) 



Jab 



1 ^, 



(2.32a) 
(2.32b) 
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where ^c" is the Coulomb gauge parameter and = (0, p). The vertex functions are 

Vff^^{PuP2,P3) ^ ^ [{Pi - P2h 9aP + {P2 " P3)a Q^j + (^3 " Pl)(3 9^^] , (2.33a) 

^jade^ebc ^^^^^^^ _ _ (2.33b) 

The last term in Eq. (|2.31j) is the trace of the logarithm of the Faddeev-Popov determinant, 
with the full inverse ghost propagator W^^. The trace runs over ghost 4- momenta and adjoint 
color indices. 

Similar to the treatment of fermions in Sec. 12.1.^ 1 now define projectors Qi, Qi for soft and 
hard gluon modes, respectively, 

^1 = Qi^ , A2 = Q2A, (2.34) 

where 

Qi(Pi,P2) = e(Agi-pi)5g)p^ , (2.35a) 
Q.2{Px.P2) = G(pi - Agi) Jg^^p^ . (2.35b) 

The gluon cut-off momentum Agi defines which gluons are considered to be soft or hard, respec- 
tively. 

I now insert A = Ai + A2 into Eq. 1)2. 3U() . The integration measure simply factorizes, 
DA = T)AiT)A2. The action ^[yl, ^'i, ^'i] can be sorted with respect to powers of the hard 
gluon field, 

S[A, ^1,^1] = S[Ai,^i, ^1] + A2J[Ai, ^1,^-1] - ^ ^1, ^1] A2 + ^2, ^-1, ^-i] . 

(2.36) 

The first term in this expansion, containing no hard gluon fields at all, is simply the action 
(|2.30b|) . with A replaced by the relevant gluon field Ai. The second term, ^2^7! contains a 
single power of the hard gluon field, where 



5A 

The first contribution, 



= Jb[Au^i,^!i\ + JiooMi\ + MM . (2.37) 

A2=0 



Jr,[A^,^^,^^] = \^>Jg^\ ^1, (2.38) 
2 V M2/a2=o 

arises from the coupling of the relevant fermions to the "modified" gluon field B, i.e., from the 
second term in Eq. (|2.30b|) . With the notation of Fig. 12.21 all diagrams corresponding to A2JB 
can be summarized into a single one, cf. Fig. 12.51 It contains precisely two relevant fermion fields, 
and ^'i. The second contribution, vTioop; arises from the terms Tr^ln^^ and Tr^/ilnW^^ in 
Eqs. ()2.30b|) . (|2.31|) . The loop consisting of irrelevant quark modes as internal lines, coupled to a 
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+ 



Figure 2.5: The term A2J7b- The hard gluon field is denoted by a dashed Une, the soft gluon fields by 
wavy lines. 




Figure 2.6: The fermionic contribution to the term A2Jioop- 
ghosts with similar topology. 



There is an additional contribution from 



single hard and arbitrarily many soft gluons, is shown in Fig. 12.61 Finally, the third contribution, 
iTv, arises from the non-Abelian vertices, cf. Fig. 12.71 

The third term in Eq. 1)2. 36() is quadratic in A2, where 



6'^S[A,^i,^i 



SA2 6A2 



Ao22+n22[^l,^'l,1'l] . 



(2.39) 



A2=0 



Here, Aq 22 is the free inverse propagator for hard gluons. Similar to the "current" J^, cf. Eq. 
1)2. 37() . the "self-energy" 1122 of hard gluons consists of three different contributions, 

n22 [^1 , ^-1 , ^1] = [^1 , ^-1 , ^-i] + Hioop [^1] + Hv [Ai] , (2.40) 

which has a diagrammatic representation as shown in Fig. 12.81 The first two contributions on 
the right-hand side of Eq. 1)2. 4U() can be expanded as shown in Figs. 12.91 and I2.1(JI Figure 12.111 
depicts the three- and four-gluon vertices contained in the last term in Eq. (|2.4Uj) . For further 
use, I explicitly give the first term. 




5^B 
SA26A2 



(2.41) 



Finally, I collect all terms with more than two hard gluon fields A2 in Eq. ()2.36|) in the 
"interaction action" for hard gluons, 5/[^i, ^21^1; ^i]- I then perform the functional integration 
over the hard gluon fields A2- Since functional integrals must be of Gaussian type in order to be 
exactly solvable, I resort to a method well known from perturbation theory. I add the source term 
^42 J2 to the action (|2.30b|) and may then replace the fields A2 in Sj by functional differentiation 
with respect to J2, at J2 = 0. I then move the factor exp{Sj[Ai,6/6J2,^i,'^i]} in front of the 
functional yl2-integral. Then, this functional integral is Gaussian and can be readily performed 
(after a suitable shift of A2), with the result 



Z = yp^'iP^'iPAi exp|5[Ai,^'i,^i]- iTr^lnA22'| 
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Figure 2.7: The term A2Jv- 




O ^ A/ X 



Figure 2.8: The term A2TI22A2 according to Eq. (|2.4U|) . The first diagram on the right-hand side 
corresponds to the term A2ngA2. The second diagram is the fermion-loop contribution to A2nioop^2; 
there is an analogous one from a ghost loop. The last two diagrams correspond to A2ll\;A2. 



X exp < Sj 



OJ2 



exp 



1 



{J + J2) A22 {J + J2) 



J2=0 



(2.42) 



The trace over lnA22^ runs over gluon 4-momenta, as well as adjoint color and Lorentz indices. 
I indicate this with a subscript "17" . Note that this result is still exact and completely general, 
since so far the manipulations of the partition function were independent of the specific choice 
(|2.35|) for the projection operators Qi,2- The next step is to derive the tree-level action for the 
effective theory of relevant quark modes and soft gluons. 



2.1.4 Tree-level effective action 

In order to derive the tree-level effective action, I shall employ two approximations. The first 
is based on the principal assumption in the construction of any effective theory, namely that 
soft and hard modes are well separated in momentum space. Consequently, momentum conser- 
vation does not allow a hard gluon to couple to any (finite) number of soft gluons. Under this 
assumption, the diagrams generated by A2[J\oop + Jv)-, cf. Fig. 12. 6| 12. 7( will not occur in the 
effective theory. In the following, I shall therefore omit these terms, so that J = Jb- Note that 
similar arguments cannot be applied to the diagrams generated by ^2(nioop + IIv)^2i cf. Fig. 
I2.1U1 12.111 since now there are two hard gluon legs which take care of momentum conservation. 

My second approximation is that in the "perturbative" expansion of the partition function 
1)2. 42() with respect to powers of the interaction action S/, I only take the first term, i.e., I 
approximate e^^ ~ 1. At this point, this is simply a matter of convenience, since I do not need 
the higher-order terms in the expansion of e"^^ for the effective theories to be discussed in Sec. 
12.21 or for the calculation of the gap parameter in Sec. 13.11 However, one can easily reinstall 
them if required by the particular problem at hand. I note in passing that the approximation 
gS/ ^ -y becomes exact in the derivation of the exact renormalization group |117j . where one 
only integrates out modes in a shell of infinitesimal thickness. 

Diagrams generated by the higher-order terms in the expansion of e^^ are those with more 
than one resummed hard gluon line. Even with the approximation e^' ~ 1, Eq. 1)2.42(1 still 
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Figure 2.9: The term A2UBA2. 




Figure 2.10: The fermionic contribution to the term A2nioop^2- 

contains diagrams with arbitrarily many bare hard gluon Hnes, arising from the expansion of 

00 (_i\n 

InA^i = lnAo,L " E (Ao,22n22)" , (2.43) 

n=l 

and from the term Jb^tiJb in Eq. (|2.42|) . when expanding 

00 

A22 = Ao,22 E (-1)" (^22 Ao,22)" • (2.44) 
n=0 

With these approximations, the partition function reads 

Z = yp^'iP^'iPAi exp{5eff [^1,^1,^1]} , (2.45) 
where the effective action is defined as 

Seff[^i,^i,vl/i] = 5A[Ai] + i'fi{go-i\+5i3[^i]}*i + ^Tr,lng22'[^i] 
-^T¥3lnA22n^i,^i,^i] 

+ ^ 'I'l, ^-i] A22[^i, ^1, ^-i] Jb\Ax, ^-1, 'fi] . (2.46) 

This is the desired action for the effective theory describing the interaction of relevant quark 
modes, and soft gluons, A\. The functional dependence of the various terms on the 

right-hand side on the fields has been restored in order to facilitate the following 

discussion of all possible interaction vertices occurring in this effective theory. 

The diagrams corresponding to these vertices are shown in Figs. I2.12II2.1^ The three- and 
four-gluon vertices contained in 5yi[74i] are displayed in Fig. 12.121 In addition, ^^[^i] contains 
ghost loops with an arbitrary number of attached soft gluon legs. The topology is equivalent 
to that of the quark loops in Fig. 12.141 and is therefore not shown explicitly. The interaction 
between two relevant quarks and the "modified" soft gluon field, corresponding to ^'i. 
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Figure 2.11: The term ylallvAa. 




Figure 2.12: The three- and four-gluon vertices in S'yi[Ai], describing the self-interaction of soft gluons 
in Eq. !^U^ . 

is depicted in Fig. 12.131 This is similar to Fig. 12. 3( except that now all gluon legs are soft. 
Diagrams where an arbitrary number of soft gluon legs is attached to an irrelevant quark loop 
are generated by Trg In Q22 ' -'^^S- 12-141 This is similar to Fig. 12. 4| but now only soft gluon 
legs are attached to the fermion loop. The diagrams generated by the loop of a full hard gluon 
propagator, TrglnA22^, are shown in Fig. 12.151 The first line in this figure features the generic 
expansion of this term according to Eq. (|2.43j) . where the hard gluon "self-energy" insertion 
1122 ) cf . Eq. 1)2.40(1 , is shown in Fig. 12.81 The second line shows examples of diagrams generated 
by explicitly inserting 1122 in the generic expansion. Besides an arbitrary number of soft gluon 
legs, these diagrams also feature an arbitrary number of relevant quark legs. If there are only 
two relevant quark legs, but no soft gluon leg, one obtains the one-loop self-energy for relevant 
quarks, cf. the second diagram in the second line of Fig. 12.151 The next two diagrams are 
obtained by adding a soft gluon leg, resulting in vertex corrections for the bare vertex between 
relevant quarks and soft gluons. The first of these two diagrams arises from the n = 1 term in 
Eq. ()2.4.S|) . while the second originates from the n = 2 term. Four relevant quark legs and no 
soft gluon leg give rise to the scattering of two relevant quarks via exchange of two hard gluons, 
contained in the n = 2 term in Eq. I|2.43() , cf . the last diagram in Fig. 12.151 This diagram was 
also discussed in the context of the effective theory presented in Refs. |l()()[ll()l] . cf. discussion in 
Sec. 12.2.21 Finally, the "current-current" interaction mediated by a full hard gluon propagator, 

^22 Jbi Fig. I2.1(H contains also a multitude of quark-gluon vertices. The simplest one is the 
first on the right-hand side in Fig. 12.161 corresponding to scattering of two relevant fermions via 
exchange of a single hard gluon. 

An important question is whether the introduction of the momentum cut-offs Aq, Agi could 
possibly spoil the gauge invariance of the effective action (|2.46|) . This is not the case, since 
gauge invariance is already explicitly broken from the very beginning by the choice of gauge in 
the gauge-fixed gluon action (|2.2|) . I do not perceive this to be a disadvantage of my approach, 
since the computation of a physical quantity requires to fix the gauge anyway. The final result 
should, of course, neither depend on the choice of gauge, nor on Aq and Agi. Note that, up to 
this point, I was not required to specify the gauge in the effective action (|2.46|) . 

The effective action ()2.46() is formally of the form (|1.4U() . The difference is that Eq. (|2.46|) 
contains more than one relevant field: besides relevant quarks there are also soft gluons. It 
is obvious that in this case there are many more possibilities to construct operators Oj which 
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Figure 2.13: The term '^igBlAi] in the effective action ()2.46(l . 



Figure 2.14: The term Tr, hi ^732^^1] in the effective action (IT^ . 

occur in the expansion H1.40() . As pointed out in the introduction, it is therefore advantageous 
to derive the effective action (|2.46p by exphcitly integrating out irrelevant quark and hard gluon 
modes, and not by simply guessing the form of the operators Oi, since then one is certain that 
one has constructed all possible operators occurring in the expansion (|1.4()|) . 

As mentioned in the introduction, the standard approach to derive an effective theory, namely 
guessing the form of the operators Oi and performing a naive dimensional scaling analysis to 
estimate their order of magnitude, fails precisely when (a) there are non-local operators, or when 
(b) there is more than one momentum scale. Both (a) and (b) apply here. As I shall show below, 
the HTL/HDL effective action is one limiting case of Eq. ()2.46|) . and it is well known that this 
action is non-local. Moreover, as is obvious from the above derivation, there are indeed several 
momentum scales occurring in Eq. H2.46|l . Let us focus on the case of zero temperature, T = 0, 
and, for the sake of simplicity, assume massless quarks, m = 0, fj, = kp. To be explicit, I employ 
the choice (|2.18|) for the projectors Vi^2- In this case, the first momentum scale is defined by the 
Fermi energy /x. The propagator of antiquarks is ~ l/(A;o + fJ- + k). If Aq, Agj < //, the exchange 
of an antiquark can be approximated by a contact interaction with strength ~ 1///, on the scale 
of the relevant quarks, Lq » 1/Aq > 1/^, or of the soft gluons, Lgi » 1/Agi > 

The second momentum scale is defined by the quark cut-off momentum Aq. The propagator 
of irrelevant quark modes is ~ l/{kQ + fi — k). On the scale Lq of the relevant quarks, not only 
the exchange of an antiquark, but also that of an irrelevant quark with momentum k satisfying 
\k — fi\ > Aq is local, with strength ~ 1/Aq. However, suppose that the quark cut-off scale 
happens to be much smaller than the chemical potential, Aq <^ fi. In this case, antiquark 
exchange is "much more localized" than the exchange of an irrelevant quark, 1/n <C 1/Aq. 

The third momentum scale is defined by the gluon cut-off momentum Agj. The propagator 
of a hard gluon is ~ On the scale Lgi of a soft gluon, the exchange of a hard gluon with 

momentum p > Agi can be considered local, with strength ~ I/^gi- I shall show below, 
in order to derive the value of the QCD gap parameter in weak coupling and to subleading 
order, the ordering of the scales turns out to be Aq < (7// <C Agi < fi. Thus, antiquark exchange 
happens on a length scale of the same order as hard gluon exchange, which in turn happens on 
a much smaller length scale than the exchange of an irrelevant quark, l//i ^ ^ 
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Figure 2.15: The term Tr^ In A22^[Ai, 5*1, vjfj^] in the effective action (|2.46|l . The first hne corresponds 
to the generic expansion (|2.43l) . with "self-energy" insertions 1122, as shown in Fig. 12.81 The second hne 
contains some examples for diagrams generated when explicitly inserting the expression for n22- 
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Figure 2.16: The term Jb^tiJb in the effective action (|2.46(> . The thick dashed line is a full hard 
gluon propagator, i.e., it has the expansion l|2.43(l . The first diagram on the right-hand side of this figure 
results from the n = term of this expansion, while the next three diagrams originate from the n = 1 
term. Even a single insertion of a hard gluon "self-energy" 1122 gives rise to a variety of diagrams. Here, 
I only show the contributions corresponding to the first diagrams in Figs. 12.91 12.101 and the three-gluon 
vertex. The last diagram arises from the second term of the expansion shown in Fig. 12.21 
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2.2 Recovery of known effective theories 

In this section I show that, for particular choices of the projectors Vi^2 in Eq. ()2.17|) . several 
well-known and at first sight unrelated effective theories for hot and/or dense quark matter are 
in fact nothing but special cases of the general effective theory defined by the action (|2.46() . 
These are the HTL/HDL effective action for quarks and gluons, and the high-density effective 
theory for cold, dense quark matter. 

2.2.1 HTL/HDL effective action 

Let us first focus on the HTL/HDL effective action. This action defines an effective theory for 
massless quarks and gluons with small momenta in a system at high temperature T (HTL), or 
large chemical potential /x (HDL). Consequently, the projectors Vi^2 for quarks are given by 

ViiK,Q) = e(Aq-fc)4jQ, (2.47a) 
V2{K,Q) = eik-A^)6^^]Q, (2.47b) 

while the projectors for gluons are given by Eq. 1)2. 35(1 . (Note that, strictly speaking, the quarks 
and gluons of the HTL/HDL effective action should also have small energies in real time. Since 
the present effective action is defined in imaginary time, one should constrain the energy only 
at the end of a calculation, after analytically continuing the result to Minkowski space.) 

The essential assumption to derive the HTL/HDL effective action is that there is a single 
momentum scale, Aq = Agi = A, which separates hard modes with momenta ~ T, or ~ fi, from 
soft modes with momenta ~ gT, or ~ gfi. In the presence of an additional energy scale T, or 
H, naive perturbation theory in terms of powers of the coupling constant fails. It was shown by 
Braaten and Pisarski [7] that, for the n-gluon scattering amplitude the one-loop term, where n 
soft gluon legs are attached to a quark or gluon loop, is as important as the tree-level diagram. 
The same holds for the scattering of n — 2 gluons and 2 quarks. At high T and small the 
momenta of the quarks and gluons in the loop are of the order of the hard scale, ~ T. This gives 
rise to the name "Hard Thermal Loop" effective action, and allows to simplify the calculation of 
the respective diagrams. At large /i and small T, i.e., for the HDL effective action, the situation 
is somewhat more involved. As gluons do not have a Fermi surface, the only physical scale which 
determines the order of magnitude of a loop consisting exclusively of gluon propagators is the 
temperature. Therefore, at small T and large /i, such pure gluon loops are negligible. On the 
other hand, the momenta of quarks in the loop are ~ fx. Thus, only loops with at least one 
quark line need to be considered in the HDL effective action. 

In order to show that the HTL /HDL effective action is contained in the effective action 1)2. 46(1 , 
first note that a soft particle cannot become hard by interacting with another soft particle. This 
has the consequence that a soft quark cannot turn into a hard one by soft-gluon scattering. 
Therefore, 

gB[A,] = gAn ■ (2.48) 

Another consequence is that the last term in Eq. (|2.46|) . Jb^^tiJb^ vanishes since Jq is identical 
to a vertex between a soft quark and a hard gluon, which is kinematically forbidden. The 
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resulting action then reads 

Si^rgcT/M^, ^1, ^i] = Sa[Ai] + (GqI^ + gAu) + \ Tr, InQ^^^Ai] 

-^TrglnA^i[Ai,^i,^i]. (2.49) 

Using the expansion H2.29() one reahzes that the term Tr^ In^^^ generates ah one-loop diagrams, 
where n soft gluon legs are attached to a hard quark loop. This is precisely the quark- loop 
contribution to the HTL/HDL effective action. 

For hard gluons with momentum ~ T or ~ the free inverse gluon propagator is Aq 22 ~ 
or ~ fi'^, while the contribution IIioop to the hard gluon "self-energy" ()2.4U() is at most of the 
order ~ g^T'^ or ~ ^^/i^. Consequently, IIioop can be neglected and 1122 only contains tree-level 
diagrams, 1122 = + IIv. Using the expansion 1)2. 43(1 of TrglnA22^, the terms which contain 
only insertions of IIv correspond to one-loop diagrams where n soft gluon legs are attached to 
a hard gluon loop. As was shown in Ref. [7j, with the exception of the two-gluon amplitude, 
the loops with four-gluon vertices are suppressed. Neglecting these, one is precisely left with 
the pure gluon loop contribution to the HTL effective action. As discussed above, for the HDL 
effective action, this contribution is negligible. 

The "self-energy" Ilg contains only two soft quark legs attached to a hard quark propa- 
gator (via emission and absorption of hard gluons). Consequently, in the expansion (|2.43|) of 
TrglnA22^, the terms which contain insertions of IIv and Ilg correspond to one-loop diagrams 
where an arbitrary number of soft quark and gluon legs is attached to the loop. It was shown 
in Ref. [Jj that among these diagrams, only the ones with two soft quark legs and no four-gluon 
vertices are kinematically important and thus contribute to the HTL/HDL effective action. I 
have thus shown that this effective action, S'htl/hdLi is contained in the effective action (|2.49j) . 
and constitutes its leading contribution, 

'S'largeT//. = SuTL/UBL + higher Orders . (2.50) 

For the sake of completeness, let us briefly comment on possible ghost contributions. Ghost loops 
arise from the term Trg/jlnW"^ in ^^^[^i]. Their topology and consequently their properties 
are completely analogous to those of the pure gluon loops discussed above. 

I conclude with a remark regarding the HDL effective action. According to Eq. 1)2. 47() . at zero 
temperature and large chemical potential, a soft quark or antiquark has a momentum k ~ gn, 
i.e., it lies at the bottom of the Fermi sea, or at the top of the Dirac sea, respectively. These 
modes are, however, not that important in degenerate Fermi systems, because it requires a large 
amount of energy /cq ~ ^ to excite them. The truly relevant modes are quark modes with large 
momenta, /c ~ /i, close to the Fermi surface, because it costs little energy to excite them. A 
physically reasonable effective theory for cold, dense quark matter should therefore feature no 
antiquark modes at all, and only quark modes near the Fermi surface. Such a theory will be 
discussed in Sec. 12.2.21 

2.2.2 High-density effective theory 

An effective theory for high-density quark matter was first proposed by Hong llODj and was 
further refined by Schafer and others jlU2l I1U31 11U41 11U5| . In the construction of this effective 
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theory, one first proceeds similarly to the discussion in Sec. 12.11 and integrates out antiquark 
modes. (From a technical point of view, this is not done as in Sec. I2.1l bv functional integration, 
but by employing the equations of motion for antiquarks. The result is equivalent.) On the 
other hand, at first all quark modes in the Fermi sea are considered as relevant. Consequently, 
in the notation of Sec. 12.11 the choice for the projectors Vi^2 would be 

Vi{K,Q) = { - ) 44)^^ (2.51a) 



K 








K 


K 











V2iK,Q) = [ 1 4)^. (2.51b) 

Also, at first gluons are not separated into soft and hard modes either. After this step, the 
partition function of the theory assumes the form (|2.1|) with Zq given by Eq. 1)2. 24() . 

In the next step, one departs from the rigorous approach of integrating out modes, as done 
in Sec. 12.11 and follows the standard way of constructing an effective theory, as explained in the 
introduction. One focusses exclusively on quark modes close to the Fermi surface as well as on 
soft gluons. However, since quark modes far from the Fermi surface and hard gluons are not 
explicitly integrated out, the effective action does not automatically contain the terms which 
reflect the influence of these modes on the relevant quark and soft gluon degrees of freedom. 
Instead, the corresponding terms have to be written down "by hand" and the effective vertices 
have to be determined via matching to the underlying microscopic theory, i.e., QCD. 

In order to further organize the terms occurring in the effective action, one covers the Fermi 
surface with "patches" . Each patch is labelled according to the local Fermi velocity, = ^kp/fi 
at its center. A patch is supposed to have a typical size Ay in radial (k) direction, and a size A_|_ 
tangential to the Fermi surface. The momentum of quark modes inside a patch is decomposed 
into a large component in the direction of v^?, the particular Fermi velocity labelling the patch 
under consideration, and a small residual component, 1, residing exclusively inside the patch, 

k = nvF + \. (2.52) 

The residual component is further decomposed into a component pointing in radial direction, 
1|| = ■ 1), and the orthogonal one, tangential to the Fermi surface, 1_l = 1 — ly. The actual 

covering of the Fermi surface with such patches is not unique. One should, however, make sure 
that neighbouring patches do not overlap, in order to avoid double-counting of modes near the 
Fermi surface. In this case the total number of patches on the Fermi surface is ~ /i^/A^. 

In the following, I shall show that the action of the high-density effective theory as discussed 
in Refs. ^100, lOQ [ESI IMl &04, 105^ is contained in the effective action TUEh . To this end, 
however, I shall employ the choice ()2.18|) and (|2.35|) for the projectors for quark and gluon 
modes, and not Eq. l^TH for the quark projectors. As in Refs. [ITInl ITTvn ITTM ITTIHl [TTHl IT^ 
the quark mass will be set to zero, m = 0. One also has to clarify how the patches covering the 
Fermi surfaces introduced in Refs. jinn|lini1 fl02, 103, 10^ 1105j arise within the effective action 
Eq. 1)2. 46(1 . It is obvious that the radial dimension Ay of a patch is related to the quark cut-off 
Aq. One simply chooses Ay = Aq. Similarly, since soft-gluon exchange is not supposed to move 
a fermion from a particular patch to another, the dimension A_|_ tangential to the Fermi surface 
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must be related to the gluon cut-off Agi. Again, one adlieres to the most simple choice A_l = A 



Since Agi < /i, this is consistent with the matching procedure discussed in Ref. jl()2j . where the 
matching scale is chosen as A_l = \/2/i (which is only slightly larger than /i). The different 
scales Aq, Agi, and fi are illustrated in Fig. 12.171 The modulus of the residual momentum 1 in 
Eq. (|2.52|1 is constrained to / < max (Aq, Agi). 



2Ap 



2Ar 



Figure 2.17: A particular patch covering the Fermi surface. The tangential dimension, Aj^, is given 
by the maximum momentum transferred via a soft gluon, Agi, while the radial dimension. Ay is defined 
by the maximum distance of relevant quark modes from the Fermi surface, Aq. Also shown is a typical 
momentum transfer 1 via a soft gluon. 

In Nambu-Gor'kov space, the leading, kinetic term in the Lagrangian of the high-density 
effective theory reads 



1 V- T , iV -D 








*i(^,vf) , 



(2.53) 
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cf. for instance Eq. (1) of Ref. |lU3j . Here, the 4- vectors 

y^ = (l,v^) , F^ee(1,-v^) (2.54) 

are introduced. The covariant derivative for charge-conjugate fields is defined as = + 
igA'^T^ . The contribution (|2.53|) arises from the term {Qqii + gAiij in Eq. ()2.46|1 . In 
order to see this, use Vi = Vi to write 

^ \ Y. *i(i.VF)7oi ( ) *i(i,vj.) . (2.55) 

In the last step, I have approximated k ~ ^-|- v^? -I, which holds up to terms of order 0(Z^//i), cf. 
Eq. (|2.52|) . This is a good approximation if the modulus of a typical residual quark momentum 
in the effective theory is / <C max (Aq, Agi) ^ ^. I have also introduced the 4- vector L^^ = (ko, I) 



and, applying the decomposition (|2.52|) . I have written the sum over k as a double sum over vp 
and 1. The latter sum runs over all residual momenta 1 inside a given patch, while the former 
runs over all patches. With this decomposition, the spinors ^'i, ^'i are defined locally on a 
given patch (labelled by the Fermi velocity vp), and depend on the 4-momentum L. Note that 
a Fourier transformation to coordinate space converts V ■ L ^ iV ■ d. 

Now consider the term ^igAu^i. Since is not diagonal in momentum space, cf. Eq. 
1)2. 16() . in principle the two quark spinors ^'i, ^'i can belong to different patches. However, I 
have chosen the tangential dimension of a patch such that a (typical) soft gluon can by definition 
never move a fermion across the border of a particular patch, |k — q| <C Agj. Therefore, both 
spinors reside in the same patch and, to leading order, k ~ q ~ vp. With these assumptions one 
may write A^4^{K - Q)A+ V ■ Aa{K - Q) -fo A+, A+^JK - Q)A- ~ F • A,(K - Q) 7o A^ . 
Then, introducing the residual momentum 1' corresponding to the quark 3-momentum q and 
defining L'^ = (go, 1')) the respective term in the effective action becomes 

1 9 X- .1, (T \ ( V ■ Aa{L - L')Ta 

In coordinate space, the sum of Eqs. (|2.55|) and (|2.56() becomes Eq. 1)2. 53() . 

Subleading terms of order 0(1/ fi) in the high-density effective theory are of the form 

^0{i^) = 
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cf. Eq. (2) of Ref. |lU3j . Here, Dj_ = {0, {l — VFVp)-Y)}, and similarly for D^j_. The commutator 
of two gamma matrices is defined as usual, af^'^ = {i/2)['j'^ and F^^^Ta = [i/ g)[Dj^^, D^y\. 
As one will see in the following, this contribution arises from the term —g^ A12 G22 •A.21 ^1 in 
Eq. (ITTO . 

First, note that, with the projectors H2.18() . the irrelevant quark propagator Q22 contains 
quark as well as antiquark modes. In order to derive Eq. ()2.57|) . however, one has to discard 
the quark and keep only the antiquark modes. In essence, this is a consequence of the simpler 
choice (|2.51|) for the projectors Vi^2 in the high-density effective theory of Refs. IIOIL ll()2| 
I1U31 ri04, 105]. In this case the propagator Q22 may be simplified. A calculation quite similar to 
that of Eqs. (|2.55j) and (|2.5(i|) now leads (in coordinate space) to 

^-1 ^-1 A ( 2ii + iV-D \ 

- ^0,2^2 + 9A22 - 70 rs ^ 2^.-^V.Dc)^ ^'''^^ 

where rs acts in Nambu-Gor'kov space. This result may be readily inverted to yield 

Utilizing the projectors 1)2. 51() . one may also derive a simpler form for gAi2 and (7^21 • Consider, 
for instance, the term ^1 5^12 ^2- One follows the same steps that led to Eqs. (|2.55|) . i.e., one 
assumes that the spinors ^1 and ^'2 reside in the same patch, such that k ~ q ~ v^?. This allows 
to derive the identity A^ 4''{K - Q)K^ ~ ^t41{K - Q), where A''^ = {0, (1 - v^v^) • A'^}. 
Now introduce the 4- vectors L'^, L'^, as in Eq. (|2.56|) . which leads to 

(2.60) 

One may add a term ^_|_ to the diagonal Nambu-Gor'kov components, which trivially vanishes 
between spinors and ^'2. This has the advantage that, in coordinate space, 

- ( 'V ) • 

i.e., this term transforms covariantly under gauge transformations, and no longer as a gauge field. 
A similar calculation for 17.421 gives the result gA2i = gAi2- Combining Eqs. 1)2. 59() and 1)2. 61() . 
the term —g^ A12 G22 A21 ^1 corresponds to the following contribution in the Lagrangian, 

1^ :f , / ^, \ 1 ^ 1 f -iV-D W ^ J \ 

x^i{X,vf) . (2.62) 

Taking only the n = term, and utilizing 7/^7^^ = g>^'^ — ia^'^ , one arrives at Eq. ()2.57|) . Note 
that my definition for transverse quantities, e.g. Aj_ = {0, (1 — v^^v^?) • A}, slightly differs from 
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that of Refs. |1UL)1 llUlj . where Aj_ = — V^V ■ A. However, both definitions agree when 
sandwiched between spinors ^'1^2 and ^'2,1- 

At order 0(1//U^), besides the n = 1 term in Eq. ()2.62|) . there are also four-fermion interaction 
terms, cf. Eqs. (3-5) of Ref. [IQSi- In the effective action (|2.46j) . these contributions arise from 
the term Jb^22Jb which originates from integrating out hard gluons. (Since this is not done 
exphcitly in the construction of the high-density effective theory in Refs. |inf)| IIDH llf)2| 11031 
I1U41I105] . this term is not automatically generated, but has to be added "by hand".) To leading 
order, this term corresponds to the exchange of a hard gluon between two quarks, cf. the first 
diagram on the right-hand side of Fig. 12.161 If the quarks are close to the Fermi surface, the 
energy in the hard gluon propagator can be neglected, and Ao,22 ^ -'^/^gi- Since 1/Ag[ > l//i^, 
the contribution from hard-gluon exchange is of order 0{l/iJ?). Four-fermion interactions also 
receive corrections at one- loop order, cf. Fig. 5 of Ref. |inij . In Eq. (|2.46j) . they are contained 
in the term TrlnA22^, see the last diagram in Fig. 12.151 

Besides the quark terms in the Lagrangian of the high-density effective theory IWl ITnTl ITTM 
ll()3( ll()4| ll()5j . there are also contributions from gluons. The first is the standard Yang- Mills 
Lagrangian —{l/'i)F^^Fl^^, cf. Eq. (1) of Ref. jlU3j . This part is contained in the term S'Ai^i] 
in Eq. (|2.46)1 . cf. Eq. (|2.2|1 . The second contribution is a mass term for magnetic gluons. 



cf. Eq. (19) of Ref. [TTH|, Eq. (18) of Ref. IM], or Eq. (27) of Ref. [102], where nig is the gluon 
mass parameter (|1.39|) . This term has to be added "by hand" in order to obtain the correct value 
for the HDL gluon polarization tensor within the high-density efi^ective theory. In Eq. ()2.46|) 
this contribution arises from the n = 2 term of the expansion ^HTM of Trlng^2^. The gluon 
polarization tensor has contributions from particle-hole and particle-antiparticle excitations. 
The latter give rise to C^g ■ While this term arises naturally within my derivation of the effective 
theory, it does not in the high-density effective theory of Refs. IMl HHSl UMl HHSI 
because only antiquarks, but not irrelevant quark modes, are explicitly integrated out. Irrelevant 
quark modes can then only be taken into account by adding the appropriate counter terms. 

Sometimes, the full HDL action is added to the Lagrangian of the high-density effective 
theory, cf. Eq. (8) of Ref. |lU3j . This procedure requires a word of caution. For instance, 
an important contribution to the HDL polarization tensor arises from particle-hole excitations 
around the Fermi surface. Such excitations are still relevant degrees of freedom in the effective 
theory. However, in order for them to appear in the gluon polarization tensor they would 
first have to be integrated out. Therefore, strictly speaking such contributions cannot occur 
in the tree-level effective action. Of course, in an effective theory one is free to add whatever 
contributions one deems necessary. However, one has to be careful to avoid double counting. 
As will be shown in Sec. 13.11 the full HDL polarization tensor will appear quite naturally in an 
approximate solution to the Dyson-Schwinger equation for the gluon propagator, however, not 
at tree-, but only at (one-)loop level. 

It was claimed in Refs. |lflH 11021 1103] that a consistent power-counting scheme within the 
high-density effective theory requires A_|_ = Ay. In contrast, I shall show in Sec. 13.11 that a 
computation of the gap parameter to subleading order requires Aq = Ay <C A_l = Agi. This 
means that irrelevant quark modes become local on a scale Iq S> 1/Aq, while antiquark modes 
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become local already on a much smaller scale, Iq ^ l//i, cf. discussion at the end of Sec. 
12.11 As mentioned in the introduction, for two different scales power counting of terms in the 
effective action becomes a non-trivial problem. While the high-density effective theory of Refs. 
[10fl. .lfll..in2. 103, 104. 105. contains effects from integrating out antiquarks, i.e., from the scale 
1/fj,, the effective action H2.46|l in addition keeps track of the influence of irrelevant quark modes, 
i.e., from physics on the scale 1/Aq ^ l//i. Since all terms in the effective action (|2.46|1 are kept, 
one can be certain not to miss any important contribution just because the naive dimensional 
power-counting scheme is invalidated by the occurrence of two vastly different length scales. 

2.3 Towards a general effective Theory 

As explained in Sec. II. 3( the calculation of any physical quantity within the framework of an 
effective theory requires the knowledge of the orders of magnitude of the operators in the action of 
the effective theory. Only after having determined all operators needed for the desired accuracy 
the actual calculation can be performed jl071 11081 1109j . Therefore, in principle, all operators 
in the general effective action (|2.42j) yet await to be power-counted. This complete analysis is 
beyond the scope of this work and remains a formidable future project. 

In the following, however, I take a first step towards this goal and begin with the class of 
quark loops that appear in the effective action after integrating out irrelevant quarks in the case 
of large /x and small T, cf. Eq. (|2.46j) and Fig. 12.141 

2.3.1 Power counting quark loops at large fj, and small T 

In [Jj the power-counting rules for HTL/HDL were generally derived. Following this analysis, 
analogous rules are set up for loops of irrelevant quark modes at large /j, and small T. However, 
I choose the projection operators (|2.18|) instead to the HTL/HDL projection operators (|2.47j) . 
Consequently, the internal irrelevant quark modes in the considered loops are antiquarks and 
quarks far from the Fermi surface with \k — fj,\ > Aq. The typical momenta of relevant quark 
modes are located much closer to the Fermi surface, |A; — /x| ~ Aq ^ Aq ^ /u, where the upper 
index n abbreviates near to the Fermi surface. These near modes can be considered as the 
natural choice for the relevant degrees of freedom of a low-energy effective theory at large /i and 
small T (in contrast to the low-momentum modes of the HDL theory, cf. the short discussion at 
the end of Sec. 12.2.11) . Therefore, a careful analysis of such loops seems worthwhile. 

The cutoff for hard gluons may be choosen as Agi < fi, while the soft gluon modes in the 
effective action have typical momenta of order Agj <C Agi, cf. Fig. 12.171 It is clear from the outset 
that for Agj < and Aq — > the considered loops must reproduce the corresponding HDLs 
discussed in Sec. 12. 2. Tl Then, however, all quarks are integrated out and one ends up with a pure 
gauge theory for gluons. It would be more intriguing to keep a thin layer of quark modes around 
the Fermi surface as explicit degrees of freedom. If one chooses Aq sufficiently smaller than Agi 
the difference between the considered loops and the HDLs will become arbitrarily small. 

An interesting ordering for all these scales could be 

~ A^ « Aq = ct>yM^-y « A^,i « Agi < ^ , (2.64) 
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where 1 > y » 5 and = (Sir/A) m^, cf. Sec. IHT^ Soft gluons being exchanged among 
relevant quarks would then have energies of order (p. For the momentum scale A^j = (p^/^M'^^^ 
these gluons would feel the effect of Landau damping [71IH]. For the choice y = 2/3 the condition 
Aq <C Ag is fulfilled. Then the considered loops of irrelevant quark modes would contain Landau 
damping already on the tree-level of the effective action. On the other hand, the effect of the 
Meissner masses would be suppressed by a factor (j)/p ~ {(p/M)'^/^ <^ 1 compared to the non- 
colorsuperconducting contributions in the loops (at least in the case of two external gluons) 

inmini- 

For T > Tc one could use this theory for the investigation of precursory effects before the 
actual onset of Cooper pairing. In the strong coupling regime at /i < 500 MeV numerical 
calculations based on the NJL-model (all gluons integrated out) showed the occurrence of low- 
energy collective modes in the form of strong diquark fluctuations with momenta close to the 
Fermi surface |11H 1112] . This is in agreement with the weak coupling analysis in |113j where, 
after extrapolating down to /i < 500 MeV, it was argued that color superconductivity could be 
a type-II superconductor with a second-order phase transition and strong pair fluctuations at 
intermediate baryonic densities. On the other hand, at asymptotically large densities in weak 
coupling gauge field fluctuations are found to dominate and drive the CSC transition first-order 
[33 I114| I115| I116j . The effective theory with the two different cutoffs proposed above might 
be useful in describing these fluctuations and estimating the relative order of magnitude of the 
diquark and the gauge field fluctuations. 

The effect of the gap in the loops of irrelevant quarks certainly would have to be incorpo- 
rated for Agj ~ (p. This would require the introduction of a bilocal source term for irrelevant 
quark modes, cf. Sec. 13.1.11 leading to non-trivial off-diagonal entries in Nambu-Gor'kov space 
and mixing quarks with charge-conjugate quarks. In the following, however, charge-conjugated 
quarks will not be used for simplicity. 

The quark propagator in the mixed representation 

The free propagator of massless quarks at chemical potential fx and temperature T is given by 

where = |/i — ek\ and the projectors on positive (e = -|-) and negative (e = — ) energy states 
defined as A^ = I (l + 6707 • . Furthermore, /cq = —i0n with the fermionic Matsubara fre- 
quency 0n = (2n -|- l)7rT. For the evaluation of the Matsubara sum in the loop it is advantageous 
|119j to use the propagators in their mixed representation 

G'o(r,k)^rEe-'=o^G'o(i^). (2.66) 

After performing the Matsubara sum in Eq. (|2.66j) in terms of a contour integral in the complex 
^Q— plane, one obtains for the range < r < /3 (l21j 

Go(r, k) = - 5] A^7o {e{ek N{ei)] e^^^^ + e(^ - ek)N{ei)e^'^^] (2.67) 

e=± 
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^ Aijo@{s[ek - ^])/,(e|)e— , (2.68) 



e,s=± 



where N(x) = [e^^'^ + 1) ^ is the thermal distribution function for fermions. For the actual 
power counting the replacement 

Q{s\k-ii\) e(s[A:-/u] - Aq) = (2.69) 

will be necessary. Before that, however, the form given by Eq. 1)2. 68(1 may be kept to simplify the 
comparison with known results. Moreover, for the following calculations it is very convenient to 
introduce 

~U{el) ^ e{s) - sN{el) , s = ±. (2.70) 

In the following, the more compact notations 

f! = fsiei), (2.71a) 
a = e^. = l/i - ej/cil (2.71b) 
= Q{s[ek-n]) (2.71c) 

are used. The fg 's obey the following relations 

fl + fts = 1, (2.72) 
fig = /^-^^^^ (2.73) 

fXm-ftin, = '-^ + SlN{e^)-S2N{e2). (2.74) 

While e discriminates positive from negative energy states (quarks from antiquarks), the sign s 
defines the direction of the evolution in imaginary time r (appearing in the exponent e"'^'^^'^ in 
the propagator): s = + denotes propagation forward and s = — backwards in imaginary time. 
Simultaneously, for modes with positive energies (e = +), s also splits the momentum space 
in one part above {s = +) and one below {s = —) the Fermi surface via 0^ = Q{s[k — ^]). 
Hence, positive energy modes above the Fermi surface always propagate forward in imaginary 
time, whereas those below always propagate backwards. For antiquarks, however, one finds 

e; ^ ei-s) , (2.75) 

which means that negative energy modes always propagate backwards in imaginary time inde- 
pendently of their momentum. 

Furthermore, = N{(-^) is the minimum of the thermal distribution functions of quarks and 
quark-holes, whereas = ex.p{\k — fi\f3) N{e^) corresponds to their maximum. Moreover, fZ = 
N{e^) ~ exp[— (A; -|- ij,)(3] is the thermal distribution of antiquarks (exponentially suppressed) 
and = exp[—{k + fi)l3]N{e^) ~ 1 that of anti-quarkholes (abundant), respectively. Eq. 1)2. 72|) 
reflects the normalisation of f^, i.e. every state with positive (negative) energy and momentum 
k is occupied, either by a(n) (anti-)quark or a(n) (anti-)quarkhole. 

In the range — /? < t < one obtains similarly 

Go{T,k) = ^ A|7oG:/f,e— ^^ (2.76) 

e,s=± 
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and therefore, for the whole range —(5<t<(5 

Go(r,k) = -sign(r) ^ A|7oe^4„(,),e— ^ . (2.77) 

e,s=it 

If -/? < r < /? and -/3 < r ± /? < /? are fulfilled, it follows with sign(r ± (3) = ± = -sign(T) 
Go(r±/3,k) = sign(r±/3) ^ A|7oe^/,':^„(,±^), e— 



sign(r) Ak7oe^/!,^„(,), e-"(^)^^^ e— 



e,s=it 



= -sign(r) Klo@tfX.ir)se-''-' 

e,s=it 

= -Go(r,k). (2.78) 

Hence, as generally required, the mixed quark propagator Gq{t,\<.) is antiperiodic with period 
P within its range of definition. 

The 4-momentum dependent part of the considered quark loops has the general form 



K 



(-1)^-^7 ^^sTYTrsU-r'-Gom , (2.79) 



where fi = (//i, • • • //tv) and Ki = K + Pi. The prefactor (—1)^"*"^ comes from the Feynman rule 
to take —Go as the internal quark propagator and to write an overall (—1) in front of each quark 
loop. The momenta of the external gluons are given by Pj+i — Pi with Pn+i = -Pi- With the 
inverse form of Eq. ()2.66|) 

Go{K) = J dTe^°^Go{T,k) (2.80) 



it follows 

j/^(Pi,---,P^) = (-l)^+^ / ^tYTtsU^^^ / dne''^^Go{n,k,) 



= (-l)"^""' / -^"^sXlr^ I dneP'^^^ Coin, k.) e'"^^-^^^. (2.81) 

i^^i ^=l i fcO 

With the relation [U] 



^°E^0= E (-l)"^5 U/3-^r, (2.82) 

i=l ) m=-oo \ i=\ ) 



2.3. Towards a general effective Theory 



55 



one can perform the integral over one of the r's trivially, say over tn. The fact that in the 
integrals < < /? leads to 



N-l 



0<TN = m/3- J2n< f3 



(2.83) 



i=l 



and consequently the sum over m runs from 1 to A'^ — 1 only. Using e^^"^^ = 1 the integration 
over TAT yields 



(-1)^+1 



r 



d^.e(p?-P°^)nGo(r,,ki) 



N-l 



N-l 



N-l 



'N-l 



m=l 



(2^ 



e,s 



-sn^n 





7V-1 



'Af-1 



iV-l ■ ^3 



■' \'^^) e,s m=l i=l ^ 



7V-1 



'N-l 



(2.84) 



where s = (si, • • • , sn) and e = (ei, • • • , cat), as well as 



q-^l-'-HN =z TV 



TV 



= Pi -P%- Si€i + SNGN 



N 



©1 ^ n©? 



(2.85) 

(2.86) 
(2.87) 



i=l 



are introduced. In the simplest case of two external gluons one has only m = 1 and therefore 

3 



j ^,Y.^^''''®1 j driftlt 



£2 „f^in 
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J (27r)3 ^ 1^1 



2 



The above Eq. ()2.88p is in accordance with the standard result for two external gluon fields 
Eq.(40) in jl24j . In order to simplify power counting arguments the energy with respect to the 
Fermi surface is introduced 

E^ = k-efi (2.89) 

Since Qg = @{s[ek — fj]) = @{seE'^) one has s = sign(e£"^) . Furthermore, e| = |// — ek\ = \eE^\ 
and therefore se = eE'^ = ek — fi. So, in the denominators one can substitute 

- siei + 5262 = -eiki + e2A;2 (2.90) 

and obtains a form, which is suitable for a power-counting analysis 

J (2^)^t^ p?-p^-ei/ci + e2fc2 

To recover a more familiar form one rewrites the numerators as 

£l_^ + siN{ei) - S2N{e2) = + eiN{E--) - e2N{E'^) (2.92) 

and finds after performing the sum over s 

which is in accordance with the standard result Eq.(44) in jl24j . It is also in agreement with 
Eq.(5.77) in jH). Moreover, in the limit /i ^ it reproduces the corresponding formular (A. 11) 
in [7] using 

AT 

lim G| = 5,,e = n'^^..e, . (2.94) 

In the cases of = 3 and 4 one obtains after rather lengthy calculations, cf. Sec. ^ Eqs. (|A.4|) 
and 



(27r)3-^ - «pU_pU_g2^2 + g3^3 

^ + siN{ei) - s^N{e^) ^ + siiV(ei) - S2A^(e2) 



(2.95) 
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and 



P2-P4- ^2^2 + 64^4 



P3-P4- ^3^3 + 64^:4 



S2^ + siNiei)-S2Nie2) 



pI 



P2 



eiki + 62/^2 



„0 

P4 



ei/ci + 64^4 



+ 



^•2 



P3 



5262 + 53^3 



^ + siN{ei)-S3N{e3) 



^ + siN{ei)-S2N{e2] 



^^3 



eiA:i + e2k2 



Power counting of quark loops at large fi and small T 

In the case of two external gluons one has Eq. ()2.9ip 

r d^k 



^ + siN{ei)-S2N{e2) 
- eifci + 62^2 



(2.96) 



(2.97) 



where e| = e^jBf^ and Xt""^ = TV, 



, which is 



-7-00 

e 



1 + 6162 ki • k2 , 

7;*° = eifei + 62^2 ' '>' = x,y,z, 

6'^ (1 - 6162 ki • k2) + ei62 (klki + ^1^2 



=x,y,z 



(2.98a) 
(2.98b) 
(2.98c) 



First note that after the replacement given in Eq. (|2.69() with T ^ Aq ^ /i one can estimate 
an upper boundary for the thermal distribution functions N{ei). In the case ej = + one has 
N{ei) < exp(— Aq/T) and for = — it is N{ei) < exp(— /u/T). Hence, all terms proportional to 
N{ei) are exponentially suppressed compared to those proportional to Sj, and will therefore be 
neglected. Then, only terms with si = —S2 contribute, where one mode in the loop propagates 
forward and one backwards in imaginary time. Making use of Eq. (|2.75|) one finds that all 
possible combinations are 



61 = - ^ si = - 
ei = +, si = + 



S2 = + ^ 62 = + (g/g+) 

S2 = - ^ 62 = + {q^ /q^) or 62 



{q/q' 



ei 



+ , Sl 



S2 



+ ^62 = + (rM) , 



where q denotes an antiquark and q^ a quark above {s = +) or below (s = — ) the Fermi surface. 
After having neglected thermal excitations for irrelevant quarks, the latter two can be identified 
with quark and quarkhole modes, respectively. Hence, one has 
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[^till^2 Q(fcl - /X - Aq) ^1^2 e(/C2 - M - Aq) 

(27r)3 1 +^ + + ^2 + A:i + A;2 

0(^2 - Ai - Aq)e(/i - Aq - fcl) - e(A;i - ^ - Aq)e(^ - Aq - /C2) 



++ 



/ci + k2 



d^k f^^,^, e(fc-/^- Aq) ^ 0(fc-M-Aq) 



(27r)3 [ +- gO + fc + |k + q| -+ _gO + _^ |k _ q| 



^.1M2 e(|k + q| - - Aq)9(/i - Aq - fc) - 9(fc - ^ - Aq)9(/X - Aq - |k + q|) 

-gO-A: + |k + q| J 

(2.99) 



The first two terms correspond to the quark/antiquark loop. After writing 0(A;2 — fJ- — Aq) = 
1— 0(/i— Aq— A;2) the 1 is identified as the ultraviolet-divergent vacuum contribution and has to be 
removed by renormalization. Since Aq <C one has Q{ii — Aq — ki) ~ Q{fi — ki), which effectively 
restricts the integration region to the Fermi sphere with the volume /x^. In the considered case 
of soft external gluons, q^,q ~ Agj <^ fx, one can approximate l/{±q^ + A; + |k + q|) « 1/2A;. 
Furthermore, one may expand T^'^^'^ in powers of ^ yielding 

T^oo = 25e,e, + ol^ , (2.100) 

T^' = T^^ = 26e,e,k' + 2ei5e,,-e,{S'' -kW)^ + o(^^^ , i = x,y,z, (2.101) 

= 26'^6e„-e,+eie2k'y + o(^^^ , i,j = x,y,z. (2.102) 

For ei = —62 the dominant contribution arises from T^''_ = T^j^ ~ 2(5^^ — k^k^), whereas all 
others are suppressed by at least one power of q/k. Then the considered integral is of order 
^3//i ~ fj.^. After multiplication with from the vertices it is ~ S'^/i^- Choosing Lgj ^ gfj, <^ fi 
this would be of the order of the tree-level propagator, which is ~ (Lgi)^. For smaller values of 
Lgi, however, the considered term would be above tree- level, since its magnitude is independent 
of the cutoff Agj. E.g. for Agj ~ g'^fj. this term would be one power of 1/g above tree-level. 

The third term in Eq. ()2.99|) contains the particle/particle-hole excitations. Considering 
again the choice of cutoffs given in Eq. (|2.64|) with Aq <C Agj one may approximate 

\k + q\-k PS q-k, (2.103) 
G(|k + q| -^-Aq)e(^-Aq-A;) 
-e(A;-/u- Aq)G(/i- Aq- |k + q|) « 5(A; - /i) q • k . (2.104) 

To obtain this estimate note that the two quark modes in the loop must be located on opposite 
sides of the Fermi surface (above and below the Fermi surface) and they must be situated within 
a shell of width Lgj around the Fermi surface since their relative momentum is given by the soft 
gluon momentum. Hence, A; ~ /U S> (7 ~ Lgj. In Eq. (|2.104|) additionally use was made of tha 
fact that due to Lgj ^ Aq the cutoff Aq in the 0-functions on the l.h.s. can be neglected, cf. 
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Figure 2.18: The integration region in the quark-quarkhole contribution in Eq. H2.99|l is the dark shaded 
area for q ^ L^j ^ Lq. The approximation in Eq. H2.104|) is equivalent to including the light shaded 
areas corresponding to relevant quark modes with |A; — /uj < Aq. 



Fig. 12.181 Note that in this way one will not recover the full result because the distribution 
functions N{€i) as well as the gap 4> have been dropped from the outset. Furthermore, for the 
choice Agj > Aq the approximation in Eq. H2.1U4|) would become invalid and the calculation of 
the quark-quarkhole loop more comphcated (but still feasible). Using Eqs. r2.1 ():-{l2.1 ()4|1 the 
quark/quarkhole excitation in Eq. (|2.99() reads 

/^VT^xrrir^, (2.105) 

which is of order //^ for all //i,^2- i-e. it is at the same order as the quark-antiquark loop and 
as the tree-level gluon propagator. In fact, Eq. ()2.64|) has turned out to be the right choice 
for the cutoffs to reproduce the HDL gluon self-energy, for details cf. |124j. For the choice 
fx ^ Agj 5/^ ^ Aq the result would nontrivially depend on the fraction x = Aq/A^j and 
converge to the HDL result in the limit x — > 0. 

Turning to the cases of three and four external gluons one first notes that, if not all internal 
lines correspond to particle and particle-hole modes the loop will be beyond tree-level. To show 
this one investigates the corresponding formulae Eqs. (|2.95() and ()2.96|) . As in the two-gluon 
case the thermal distribution functions may be neglected. Assume that the line with label 1 
is an antiparticle (si = — ) and all other lines are particles {sj = -|- , j > 1). Then, using the 
approximation Eq. (|2. 103)1 all energy denominators are of order Agj except for those inside the 
angular brackets, which are of order fj,. Combining the latter, the energy denominators standing 
directly in front of the angular brackets are canceled. Then in the case of three external gluons 
the integral multiplied by is of order g^n^/ij,"^ ~ (7^Ag[, which is suppressed by g compared 
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to the tree- level, which is gA^^. In the case of 4 external gluons the integral multiplied by g'^ 
is of order (5r^/z^/(/i^Agj) ^ g^ and thus is also suppressed by g compared to tree-level, which is 
g^. Analogously one verifies that also all other integrals containing antiparticles are suppressed 
by at least one power of g in comparison to tree-level. So, even though the integral region is 
of order fj.^ in the case that one mode corresponds to an antiparticle and all others are quarks 
above the Fermi surface no HDLs are produced. 

In this context it seems worth mentioning that in [Jj it was claimed that loops with external 
gluons are HTLs if all but one of the A^ — 1 occurring energy denominators carry energies with 
opposite sign (and therefore are of order Lg ~ gT). It was argued there that in this special case 
the distribution functions appear as a sum which is of order 1 in comparison to the difference 
of two distribution functions, which is only of order g, appearing when all energy denominators 
are of order gT. This does indeed compensate the one energy denominator which is of the order 
of the hard scale T. If one, however, sums up all terms of this type systematically as it has been 
done above in the present treatment, one finds that these terms will cancel, up to a remaining 
part which is suppressed by a factor of g compared to tree-level. Hence, also in the regime of 
high temperatures all energy denominators must be of order of the soft scale Lgj in order to 
produce a HTL. 

The loops containing only quark and quarkhole modes can simplified in virtue of approxi- 
mations Eqs. (|2.1U3|) and H2.1U4|) to 



1/^2/^3 



{Qi,Q2 



(2vr) 



3 ^+++ 



(9i + 4) - (qi + q2) • k [gO - q2 • k 



q2 • k 



qi • k 



qi • k 



(2.106) 



and 



J^i-^^(Ql,Q2,Q3 



7^1 •■■M4 



9i + 92 + 93 - (qi + q2 + qa) • k 



9i + 92 - (qi + q2) • k 



qi • k 



q2 • k 



- qi • k - q2 • k_ 



+ 



93 - q3 • k 



(q2 + q3) • k 



_ q2 • k 

92 + 93 - (q2 + q3) • k ^2 - q2 • k 



(2.107) 



After multiplying the expression for the loop of three external gluons with (/^ it is of order 
(/3/i^/Ag[ ~ fl'Agj, which is of the order of the corresponding tree-level diagram. Although in the 
angular brackets differences appear, no cancellations happen that effectively reduce the overall 
order of magnitude. E.g. the difference in the brackets in Eq. (|2.106|) vanishes for Q2 —Qi ■ 
In that case the energy denominator in front of the brackets disappears as well as preserving 
the overall order of magnitude of the loop. 

Analogously, one finds for the loop with 4 external gluons g^fJ-'^ / (Agj)'^ ~ g'^, i.e. it contributes 
also at tree-level. Again the various terms cannot cancel each other for any configuration of 
external gluon momenta {Qi, • • • , Q4). If, e.g., again Q2 — > —Qi , the term in front of the first 
angular bracket will keep the order of magnitude of the term. 
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One may conjecture as in ^7,, how the above discussions of quark loops with 3 and 4 external 
gluons can be generalized to loops with any number > 3 of external gluons. Quark loops with 
N external gluon legs with all internal lines corresponding to particle or particle-hole modes get 
support from a thin layer around the Fermi surface with volume proportional of /^^A^j. The 
integrand consists of A^ — 1 energy denominators which are all of order 1 / Agj and the trace over 

Dirac space Tff' which is dimensionless. Multiplying the loop with it is of order 



9'^ JX§N^ = 9''-'iKir'' [j^J = tree-level x (^^ j . (2.108) 

Hence, if Agj ~ g^ the loop is at the order of the corresponding tree-level diagrams, gf^~^(Ag[)^~^, 
and therefore has to be resummed into the free gluon propagator or into a gluon vertex, respec- 
tively. Therefore, at this scale the free gluon propagation and the propagation via particle/hole 
and particle/antiparticle excitations are of equal importance. At scales A^j < gfj, these loop 

diagrams even dominate the corresponding tree diagrams by a factor ~ (^gfi/A^^^^ . Gluons at 
scales Agj < gfj, prefer to propagate as quark/antiquark or as quark/quarkhole pair rather than 
to propagate freely. In other words, decreasing the scale A^j sufficiently below gfj, the free prop- 
agation becomes a small correction to the "loop propagation" and the tree-level gluon vertices 
become a small correction to the dominant 3 and 4 gluon loops. This is important when con- 
sidering the static limit P ^ for gluons, for example. These loops of irrelevant quark modes 
are therefore relevant operators in the effective action, cf. Sec. 11.31 

At harder scales A^j > g^ quark loops are small corrections to the free propagation and the 
3 and 4 gluon vertices as in vacuum perturbation theory. Note that in the case of hot quark 

matter with negligible chemical potential one has analogously HTL ~ tree-level x (^gT / A^-^ . 

All other types of quark loops will be suppressed by powers of Ag[//i. If Agj > gfj, they are 
below tree-level and therefore part of the perturbative expansion. For Agj < gfi, however, those 
terms which are suppressed by Agj//i still dominate the corresponding tree diagram by a factor 
of order ^r^/Agj. They are therefore still relevant. Those terms which are suppressed by (Agj//x)^ 
scale as tree-level and consequently are marginal. Diagrams, which are suppressed by (A^i/fi)^ 
are always beyond tree-level and therefore perturbative. 

Before applying the effective action Eq. (|2.46|) to some perturbative quantity one first has 
to extend this analysis to all the remaining diagrams contained in the effective action. In the 
following Chapter 131 however, it will be demonstrated how the effective action can be combined 
with the C JT-formalism jl06j in order to facilitate the calculation of the non-perturbative color- 
superconducting gap-parameter. 
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Application to 

Color superconductivity 

3.1 Calculation of the QCD gap parameter 

In this section, I demonstrate how the effective action Eq. (|2.46jl derived in Sec. l2.1l can be apphed 
to compute the gap parameter of color-superconducting quark matter to subleading order. To 
this end this effective action is combined with the CJT-formahsm, which will be introduced 
below. The resulting subleading order gap equation, Eq. (|3.64|) is well-known. However, as 
will be explicitly shown, the derivation given here is substantially more rigorous as the known 
derivations in full QCD. One of the most important advantages is that the self-consistency is 
required only for relevant fields and propagators. Furthermore, the expansion of the various 
terms in the gap equation in powers of Aq/Agi ~ g emerges as a powerful tool to rigorously 
identify leading, subleading and sub-subleading order contributions to the gap parameter. 

3.1.1 CJT formalism for the effective theory 

The gap parameter in superconducting systems is not accessible by means of perturbation the- 
ory; one has to apply non-perturbative, self-consistent, many-body resummation techniques to 
calculate it. For this purpose it is convenient to employ the CJT formalism |lU6j . The first step 
is to add source terms to the effective action (|2.46|) . 

5efr[yli,^'i,^'i] — > Ses[Ai,^i,'i'i]+JiAi + ^AiKiAi + ^{^iHi+Hi^i + ^i}Ci^i) , 

(3.1) 

where I employed the compact matrix notation defined in Eq. 1)2. 15() . Ji, Hi, and Hi are local 
source terms for the soft gluon and relevant quark fields, respectively, while Ki and ICi are bilocal 
source terms. The bilocal source /Ci for quarks is also a matrix in Nambu-Gor'kov space. Its 
diagonal components are source terms which couple quarks to antiquarks, while its off-diagonal 
components couple quarks to quarks. The latter have to be introduced for systems which can 
become superconducting, i.e., where the ground state has a non-vanishing diquark expectation 
value, {tpitpi) / 0. 



63 



64 



Chapter 3. Application to color superconductivity 



One then performs a Legendre transformation with respect to all sources and arrives at the 
C JT effective action IMl HHl 



Here, S'eff [A, is the tree-level action defined in Eq. H2.46() . which now depends on the 

expectation values A = (Ai), ^ = (^i), and = (^'i) for the one-point functions of soft 
gluon and relevant quark fields. In a slight abuse of notation, I use the same symbols for the 
expectation values as for the original fields, prior to integrating out modes. This should not lead 
to confusion, as the original fields no longer occur in any of the following expressions. 

The quantities D^^ and in Eq. (|,3.2|) are the inverse tree-level propagators for soft gluons 
and relevant quarks, respectively, which are determined from the effective action Sqs, see below. 
The quantities A and Q are the expectation values for the two-point functions, i.e., the full 
propagators, of soft gluons and relevant quarks. The functional r2 is the sum of all two-particle 
irreducible (2PI) diagrams. These diagrams are vacuum diagrams, i.e., they have no external 
legs. They are constructed from the vertices defined by the interaction part of Scs, linked by 
full propagators A, Q. The expectation values for the one- and two-point functions of the theory 
are determined from the stationarity conditions 



The first condition yields the Yang-Mills equation for the expectation value A of the soft gluon 
field. The second and third condition correspond to the Dirac equation for ^ and ^, respectively. 
The effective action (j2.46|l contains a multitude of terms which depend on A,^,^, and thus 
the Yang-Mills and Dirac equations are rather complex, wherefore I refrain from explicitly 
presenting them here. Nevertheless, for the Dirac equation the solution is trivial, since ^i, ^'i 
are Grassmann-valued fields, and their expectation values must vanish identically, ^ = (^i) = 
^ = (^i) = 0. On the other hand, for the Yang- Mills equation, the solution A is in general 
non-zero but, at least for the two-flavor color superconductor considered here, it was shown 
jl21l 112^ to be parametrically small, A ~ (fp'/{g^^), where 4> is the color-superconducting gap 
parameter. Therefore, to subleading order in the gap equation it can be neglected. 

The fourth and fifth condition (|3.3() are Dyson-Schwinger equations for the soft gluon and 
relevant quark propagator, respectively. 



r[A^',^,A,g] = S,Q[A,^,^]--TTg\n^ 
+ lTr,lng-i + iT:,(G-i 



--iTr,(z.-A-l) 
g-l) +r2[^,^',^,A,g] . (3.2) 



= 



_5T _5T _6T _ 5T 
jA~'5^~'5^~JK~'5g' 



(3.3) 



A-i 



(3.4a) 
(3.4b) 



where 



n 



8T2 



(3.5a) 



(5A^ 

<5r2 



(3.5b) 
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are the gluon and quark self-energies, respectively. The Dyson-Schwinger equation for the rele- 
vant quark propagator is a 2 x 2 matrix equation in Nambu-Gor'kov space, 

(3.6) 

where S+ is the regular self-energy for quarks and S~ the corresponding one for charge-conjugate 
quarks. The off-diagonal self-energies the so-called gap matrices, connect regular with 
charge-conjugate quark degrees of freedom. A non-zero corresponds to the condensation of 
quark Cooper pairs. Only two of the four components of this matrix equation are independent, 
say and ^>+, the other two can be obtained via [G-J-i+S" = C{[G+]-^+J:+}^C-\ 

= 7o[<l'^]^7o- Equation (|3.6|) can be formally solved for Q |123j . 

(3.7) 



-1 




where 



g±= + + (3.8) 



is the propagator describing normal propagation of quasiparticles and their charge-conjugate 
counterpart, while 

H± = - {[G^]-^ + S^) (3.9) 

describes anomalous propagation of quasiparticles, which is possible if the ground state is a 
color-superconducting quark-quark condensate, for details, see Ref. |13) . 
The tree-level gluon propagator is defined as 

, . (3.10) 

5A5A ^ ^ 

Since I ultimately evaluate the tree-level propagator at the stationary point of F, Eq. (|3.3j) . 
where ^' = ^' = 0, I may omit all terms in S^s, Eq. 1)2. 46() . which are proportional to the quark 
fields. The only terms which contribute to the tree-level gluon propagator are therefore 

Using the expansions 1)2. 27() . (|2.29l) . (|2.43|) . and 1)2. 44() . and exploiting the cyclic property of the 
trace, one finds 

In order to proceed, note that the Dyson-Schwinger equations (|3.4() are evaluated at the sta- 
tionary point of the effective action, where ^ = ^' = 0, A ~ 0. For ^ = 0, the first term yields 
the free inverse propagator for soft gluons, Aq^^^, cf. Eq. ()2.31|) . plus a contribution from the 
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Faddeev-Popov determinant, ((^^Tr^/jlnW ^ /6A6A)a=o- The contributions from the three- and 
four-gluon vertex vanish for ^ = 0. Furthermore, according to Eq. (|2.16jl . 

This is a matrix in fundamental color, flavor, and Nambu-Gor'kov space, as well as in the space 
of quark 4-momenta K,Q. It is a vector in Minkowski and adjoint color space (F carries a 
Lorentz- vector and a gluon color index), as well as in the space of gluon 4-momenta P. One 
evaluates {SG22/SA)a=q using the expansion (\2.'27^ . Only the term for n = 1 survives when 
taking ^ = 0. For ^' = ^' = 0, I have Us = 0, cf. Fig. 12.91 and I only need to consider 
1122 = Hioop + Hy. Then, the term V*-^-* = {6IIv/5A)a=o corresponds to a triple-gluon vertex, 
cf. Fig. 12.111 where two hard gluons couple to one soft gluon. The term (5nioop/'^^)A=o is a 
correction to this vertex: it couples two hard gluons to a soft one through an (irrelevant) quark 
loop, cf. Fig. 12.101 According to arguments well known from the HTL/HDL effective theory, 
this vertex correction can never be of the same order as the tree-level vertex V^^^, since the 
two incoming gluons are hard. I therefore neglect (6Il\oop/dA)A=o in the following. Similarly, 
V^'''' = {6'^II\;/6A6A)a=o is a four-gluon vertex, cf. Fig. 12. Ill where two hard gluons couple to 
two soft ones, and {6'^IIioop/SA5A)a=o is the one-(quark-)loop correction to this vertex, cf. Fig. 
I2.1U1 Applying the same arguments as above I only keep V*-^-*. Arguments from the HTL/HDL 
effective theory also show that to leading order one may approximate A22 — ^0,22- Finally, 
utilizing the same arguments I approximate 6A22/SA ~ — Ao,22V*-'^^Ao,22- Then, the inverse 
tree-level gluon propagator of Eq. (|3.12|1 becomes 

D-^ = Ao + ^ Tr, (^0,22 f ^0,22 f ) - ^ Tr^ (Ao,22 V^^) Ao,22 V^^')) + ^ TV^ (Ao,22 V(^) 



SA6A 



(3.14) 

A=0 



The second term represents an (irrelevant) quark- loop, while the third term is a hard gluon loop. 
The fourth term is a hard gluon tadpole. Finally, the last term in Eq. (|3.14|) corresponds to 
a ghost loop necessary to cancel loop contributions from unphysical gluon degrees of freedom. 
Note that in the effective theory loop contributions involving irrelevant quarks and hard gluons 
occur already in the tree-level action (|2.46|) . Therefore, such loops also arise in the inverse tree- 
level propagator 1)3. 14|) for the soft gluons of the effective theory. For the projection operators 
(|2.35|) and (|2.47j) the inverse tree- level propagator (|3.14)) is precisely the HTL/HDL-resummed 
inverse gluon propagator. For small temperatures, T <C ^, the contribution from the gluon and 
ghost loops is negligible as compared to that from the quark loop, 

D-^ ~ Aq ^1 + ^ Tr, (^0,22 f ^0,22 f ) . (3.15) 

The inverse tree-level quark propagator is defined as 

Q-i = _2 ^it^ ' ' J . 3.16 

6"^ 6"^ 
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For ^ = ^' = 0, the last term in Eq. (|2.46|) does not contribute to G^^ because it has at least 
four external quark legs, and the two functional derivatives 6/6^, 5/5^ amputate only two of 
them. The first and the third term in Eq. H2.46() do not depend on at all, therefore 

o-' = e„:!.+sfiM + ^^^|J^. (3.17) 

Using the expansion formulae (|2.43j) and (|2.44|) and the fact that 1122 depends on only 
through Ilg, I obtain 

fai5^=^,fA,,|!I|). (3.18) 

I have exploited the fact that this expression is evaluated at ^' = ^ = 0, i.e., terms with external 
quark legs will eventually vanish. The trace runs only over adjoint colors, Lorentz indices, and 
(hard) gluon 4-momenta. Since A22 is a hard gluon propagator, the contribution from 1122 to 
A22 may be neglected to the order I are computing, and I may set A22 — Ao,22- Furthermore, 
{6^Uis/5'^6'^)a=o = -5^ f ^0,22 f , cf. Fig. EH At 'J' = ^' = 0, A ~ I are left with 

= So,li - 9^ Trg (Ao,22 f ^0,22 f ) . (3.19) 

As was the case for the tree-level gluon propagator, also the tree-level quark propagator receives 
a loop contribution; here it arises from a loop involving an irrelevant quark and a hard gluon line. 
The term r^o,22r under the gluon trace remains a matrix in the quark indices, i.e., fundamental 
color, flavor, Dirac, and quark 4-momenta. 

I now proceed to solve the Dyson-Schwinger equations (|3.4|) for the soft gluon and relevant 
quark propagator. To this end I have to determine r2. Of course it is not feasible to consider all 
possible 2PI diagrams. The advantage of the C JT formalism is that any truncation of r2 defines 
a meaningful, self-consistent many-body approximation for which one can solve the Dyson- 
Schwinger equations (|3.4|) . In the truncation of r2 I only take into account 2-loop diagrams 
which are 2PI with respect to the soft gluon and relevant quark propagators A, Q, 

22 2 

r2 = Tr,,g (gtgtA)-^ Tr,,^ (g f ^0,22 f a) - ^ Tr,,^ (gtgt Ao,22) • (3.20) 

The traces now run over quark as well as over gluon indices. Consider, for instance, the term 
^ r ^ P. It is a matrix in the space of fundamental color, flavor, Dirac and quark 4-momenta, of 
which the trace is taken through Tr^. In addition, due to the two factors P it carries two Lorentz- 
vector, adjoint-color, and gluon-4-momenta indices. The trace Tr^ contracts these indices with 
the corresponding ones from the gluon propagator A. 




Figure 3.1: Diagrammatic representation of Eq. H3.2()|l . 
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The diagrams corresponding to Eq. ()3.2U() are shown in Fig. 13.11 The first two terms are 
constructed from the quark-gluon couphng ~ ^ gB^. Using Eq. (|2.26|) . one may either obtain 
an ordinary quark-gluon vertex ~ ^ ^ involving one soft gluon and two relevant quark legs, 
or a vertex ~ g'^ ^ AQ22 with (at least) two soft gluon legs and two relevant quark legs. To 
lowest order, I approximate Q22 — ^0,22, which neglects vertices with more than two soft gluon 
legs. Taking two ordinary quark-gluon vertices and tying them together to obtain a 2PI 2-loop 
diagram, I arrive at the first term in Eq. (|3.2U|) . or the first diagram in Fig. 13.11 Taking one 
of the two-gluon-two-quark vertices and tying the legs together, one obtains the second term 
in Eq. (|3.2()j) . or the second diagram in Fig. 13. H respectively. Finally, the third term/diagram 
arises from the last term in Eq. ()2.46|) . To lowest order, this corresponds to a four-quark vertex 
~ fl'^ ^' r ^ Ao,22^ r ^- Tying the quark legs together to form a 2PI diagram, one obtains the 
corresponding term/diagram in Eq. H3.2U() /Fig. 13.11 

The combinatorial factors in front of the various terms in Eq. 1)3. 2U() are explained as follows. 
In the first diagram, there are two ordinary quark-gluon vertices. According to Eq. (|2.46|) . each 
diagram comes with a factor 1/2. Moreover, since there are two vertices, the diagram is, in 
the perturbative sense, a diagram of second order, which causes an additional factor 1/2 p3] . 
Finally, there are two possibilities to connect the quark lines between the two vertices. In total, 
I then have a prefactor —(1/2)^ x 1/2 x 2 = —1/4, where the minus sign arises from the fermion 
loop. The second diagram arises from the two-quark-two-gluon vertex, which already comes 
with a prefactor —1/2 in Eq. (|2.46j) . It is perturbatively of first order, and there is only one 
possibility to tie the quark and gluon lines together, so there is no additional combinatorial 
factor (and no additional minus sign) for this diagram. Finally, the third diagram arises from 
the four-quark vertex, {1/2)J'qIS.q 22-Jb-, in Eq. (|2.46|) . This vertex comes with a factor 1/2 and 
is perturbatively of first order. However, there are two additional factors 1/2 residing in Jq, 
since Jb ~ (1/2)^ F^*, cf. Eq. (|2.38j) . Again, there are two possibilities to tie the quark lines 
together, so that, in total, I have a prefactor —1/2 x (1/2)^ x 2 = —1/4, where the minus sign 
again stands for the quark loop. 




Figure 3.2: Diagrammatic representation of FJ* , after decomposing quark lines into relevant and 
irrelevant, as well as gluon propagators into soft and hard contributions. 

At this point, it is instructive to compare Eq. 1)3. 20() . in the effective theory with p^^^ 
which one would have written down in QCD at the same loop level, r^*^^ would be equivalent to 
the first diagram of Fig. 13.11 but now the quark and gluon lines represent the full propagators for 
all momentum modes, relevant and irrelevant as well as soft and hard. In order to compare with 
F2 of the effective theory, I decompose the quark propagators into relevant and irrelevant modes, 
and the gluon propagator into soft and hard modes. One obtains the six diagrams shown in 
Fig. 13.21 The first three are precisely the same that occur in F2 of the effective theory, including 
the combinatorial prefactors. The last three diagrams do not occur in F2 of the effective theory, 
because they are not 2PI with respect to the relevant quark propagator Q and the soft gluon 
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propagator A. Nevertheless, they are still included in the CJT effective action of the effective 
theory, Eq. (|l'{.2j) : opening the relevant quark line of the fourth diagram, I recognize the loop 
contribution to the tree- level quark propagator G~^, cf. Eq. ()3.19|) . Now consider the fifth term 
in Eq. (|3.2jl : here, this loop contribution to is multiplied by Q and traced over, which 
yields the fourth diagram in r^*""^. Similarly, opening the soft gluon line of the fifth diagram, I 
identify this diagram as the irrelevant quark-loop contribution to the tree-level gluon propagator 
D~^, cf. Eq. 1)3. 15() . The third term in Eq. (|3.2^ . where this contribution is multiplied by A and 
traced over, then yields the fifth diagram of E^ ^. Einally, the sixth diagram resides in the term 
~ Tig In A22^ of the tree- level effective action Scs, cf. Eig. 12.1.'^ Therefore, in principle, the CJT 
effective action 1)3. 2() for the effective theory contains the same information as the corresponding 
one for QCD. However, while in QCD self-consistency is maintained for all momentum modes 
via the solution of the stationarity condition (|3.3|) . in the effective theory self-consistency is only 
required for the relevant quark and soft gluon modes. These are the only dynamical degrees 
of freedom in the CJT effective action; the irrelevant fermion and hard gluon modes, which 
were integrated out, only appear in the vertices of the tree-level action (|2.46j) . In this sense the 
effective theory provides a simplification of the full problem. 

3.1.2 Dyson-Schwinger equations for relevant quarks and soft gluons 

After having specified r2 in Eq. H3.2U() I am now in the position to write down the Dyson- 
Schwinger equations (|3.4j) explicitly. Eor the full inverse propagator of soft gluons I obtain with 
Eqs. (|n3a|), and (1^1) 

Tr^ (^0,22 f ^0,22 f ) + 2 Trg (g t Go,22 f ) + Tr, f f)] . (3.21) 

The first term in square brackets takes into account the effect of quark-antiquark excitations 
as well as quark-hole excitations far from the Eermi surface. The second term is the contri- 
bution from excitations where one quark is close to the Eermi surface (a relevant quark) while 
the second is far from the Eermi surface or an antiquark (an irrelevant quark). The relevant 
quark propagator Q can have diagonal elements in Nambu-Gor'kov space, corresponding to nor- 
mal propagation of quasiparticles, as well as off-diagonal elements, corresponding to anomalous 
propagation of quasiparticles, cf. Eq. 1)3. 7|) . However, in the second term in square brackets the 
latter contribution is absent, because ^0,22 is purely diagonal in Nambu-Gor'kov space, cf. Eq. 
1)2. 12() . This is different for the last term in square brackets, which corresponds to quark-hole 
excitations close to the Fermi surface. Both quark propagators have to be determined self- 
consistently and may have off-diagonal elements in Nambu-Gor'kov space. Consequently, the 
trace over Nambu-Gor'kov space gives two contributions, a loop where both quarks propagate 
normally, and another one where they propagate anomalously. Diagrams of this type have been 
evaluated in Ref. |124j and lead to the Meissner effect for gluons in a color superconductor. 

For the full inverse propagator of relevant quarks I obtain with Eqs. ()3.4b|) . (|3.5b|) . (|3.19j) . 
and (Hnnj) 

g-' = 

Go,u ~ [Tr^ (Ao,22 f e?o,22 f ) + Tr^ (a f ^0,22 f ) + Tr^ (Ao,22 E ^ f) + Tr^ (a f ^ f)" . 

(3.22) 



^ — ^0,11 ^ 2 



70 



Chapter 3. Application to color superconductivity 



The first two terms in square brackets do not have off-diagonal components in Nambu-Gor'kov 
space. They contribute only to the regular quark self-energy. The other two terms in square 
brackets have both diagonal and off-diagonal components in Nambu-Gor'kov space. The diagonal 
components contribute to the regular quark self-energy, in particular, the fourth term leads to 
the quark wave- function renormalization factor computed first in Ref. [125j . It gives rise to 
non- Fermi liquid behavior jHZI- The off-diagonal components enter the gap equation for the 
color-superconducting gap parameter. 

The system of Eqs. (|3.21|) and (|3.22j) has to be solved self-consistently for the full prop- 
agators of quarks and gluons. However, as was shown in Ref. |126| . in order to extract the 
color-superconducting gap parameter to subleading order it is sufficient to consider the gluon 
propagator in HDL approximation; corrections arising from the color-superconducting gap in 
the quasiparticle spectrum are of sub-subleading order in the gap equation. For the present 
purpose this means that it is not necessary to self-consistently solve Eq. (|3.2H) together with 
Eq. ()3.22() : one may approximate Q on the right-hand side of Eq. H3.21() by ^0,11 • In essence, this 
is equivalent to considering only the first term on the right-hand side of Eq. (|3.22j) when solving 
Eq. ()3.21|) . Of course, under this approximation the effect of the regular quark self-energy (lead- 
ing to wave-function renormalization) and of the anomalous quark self-energy (which accounts 
for the gap in the quasiparticle excitation spectrum) are neglected. 

With this approximation, and using Qq = Go,ii © ^0.22; one may combine the terms in Eq. 
(TTTrj) to give 




(3.23) 



Taking the gluon cut-off scale Agi to fulfill gn <C Agi < /i, soft gluons are defined to have momenta 
of order gfj,. I compute the fermion loop in Eq. (|3.23j) under this assumption (taking the soft gluon 
energy to be of the same order of magnitude as the gluon momentum) . I then realize that the soft 
gluon propagator determined by Eq. (|3.23j) is just the gluon propagator in HDL approximation. 
I indicate this fact in the following by a subscript, A = Ahdl- Armed with this (approximate) 
solution of the Dyson-Schwinger equation (|3.2H) I now proceed to solve Eq. (|3.22|) . I consider 
the two independent components [G"*"]"^ -|- S+ and in Nambu-Gor'kov space separately. 
Due to translational invariance it is convenient to define [G^]^^{K, Q) = {1/T)[G^\^^{K) (5^q, 
S+(A', Q) = (l/T)S+(if) 5^^]q and using Eqs. (fTTnj) . (|2T2|l . (l3T3|l . I obtain the Dyson- 

Schwinger equation for [G^]^"*^ -|- S^, 

[G+]-\K) + ^+{K) = 

\Gtiir\K) - 5' ^ E { [^0,22]^,^ {K-Q) + [Ahdl]^,^ {K-Q)^ ^,T^ GoMQ) IuT' 

^ E { [^0,22]^,^ {K-Q) + [AhdlI^T {K-Q)] 1,T^ Q^{Q) JuT" . (3.24) 

Note that the first sum over Q runs over irrelevant quark momenta, < (7 < /i — Aq and 
fi + Aq<q< 00, while the second sum runs over relevant quark momenta, /i — Aq < g < fi + A^. 
There is no double counting of gluon exchange contributions, since the hard gluon propagator 
Ao,22 has support only for gluon momenta |k — q| > Agi, while the HDL propagator is restricted 
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to gluon momenta |k — q| < Agi. To subleading order in the gap equation, I do not have to 
solve this Dyson-Schwinger equation self-consistently. It is sufficient to use the approximation 
~ CfQ^ii on the right-hand side of Eq. ()3.24|) and to keep only the last term which, as discussed 
above, is responsible for non-Fermi liquid behavior in cold, dense quark matter. The net result 
is then simply a wave- function renormalization for the free quark propagator Gq^^ |125j . 

70 - 7 • k , (3.25) 

where g = g/{3V2Tr) and = (37r/4)m^, with the gluon mass parameter rUg defined in Eq. 
H1.39() . Neglecting effects from the finite life-time of quasi-particles |123| . which are of sub- 
subleading order in the gap equation, the wave-function renormalization factor is 

Z(feo)= + . (3.26) 

Due to translational invariance, it is convenient to define ^~^{K,Q) = (l/T) ^~^{K) S^^^q and 
'r.~^{K,Q) = T3^{K) 6^^^Q, and the Dyson-Schwinger equation for ^^{K) becomes 

<^+{K) = 5^ ^ E { [^0,22]^,^ {K-Q) + [AudlZ (K - Q) } 7m (T»)^ H+(Q) 7.T^ • (3.27) 

Here, the sum runs only over relevant quark momenta, /x — Aq < q < fj, + Aq. This is the 
gap equation for the color-superconducting gap parameter within my effective theory. There is 
no contribution from irrelevant fermions, since their propagator is diagonal in Nambu-Gor'kov 
space. 

While the gluon cut-off was taken to be Agi < fi, so that soft gluons have typical momenta 
of order gfi, so far I have not specified the magnitude of Aq. In weak coupling, the color- 
superconducting gap function is strongly peaked around the Fermi surface |53l 1551 IBl] . For a 
subleading-order calculation of the gap parameter, it is therefore sufficient to consider as relevant 
quark modes those within a thin layer of width 2Aq around the Fermi surface. For the following, 
my principal assumption is < gfj, <^ Agi < As I shall see below, this assumption is crucial 
to identify sub-subleading corrections to the gap equation ()3.27|) . which arise, for instance, from 
the pole of the gluon propagator. Note that this assumption is different from that of Refs. 
jlUll I1U3] . where it is assumed that Aq ~ Agi. 

For a two-flavor color superconductor, the color-flavor-spin structure of the gap matrix is 

[ini 

^+{K) = J3T275 A,^ e(Aq -\k-fi\) m) , (3.28) 

where {J'i)ij = —ieij3 and {T2)fg = —i^fg represent the fact that quark pairs condense in the 
color-antitriplet, flavor-singlet channel. The Dirac matrix 75 restricts quark pairing to the 
even-parity channel (which is the preferred one due to the U{1)a anomaly of QCD). In the 
effective action ()2.46|1 . antiquark and irrelevant quark degrees of freedom are integrated out. 
The condensation of antiquark or irrelevant quark pairs, while in principle possible, is thus not 
taken into account; the bilocal source terms in Eq. ()3.1() only allow for the condensation of 
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relevant quark degrees of freedom. The condensation of antiquarks or irrelevant quarks could 
also be accounted for, if one introduces bilocal source terms already in Eq. (|2.6|) . i.e., prior to 
integrating out any of the quark degrees of freedom. While there is in principle no obstacle 
in following this course of action, it is, however, not really necessary if one is interested in a 
calculation of the color-superconducting gap parameter to subleading order in weak coupling: 
antiquarks contribute to the gap equation beyond subleading order |127j , and the gap function for 
quarks falls off rapidly away from the Fermi surface, i.e., in the region of irrelevant quark modes, 
and thus also contributes at most to sub-subleading order to the gap equation. Consequently, 
the Dirac structure of the gap matrix (|H.27|) contains only the projector Aj^ onto positive energy 
states. The theta function accounts for the fact that the gap function (f){K) pertains only to 
relevant quark modes. 

Inserting Eq. 1)3. 25(1 and the corresponding one for + as weh as Eq. (TOHl) . into 

the definition (|3.9|) for the anomalous quark propagator, one obtains 

= J3T275 Aq e(Aq - k - Ml) [q^/zt!S]'-el • ^^-^^^ 

One now plugs this expression into the gap equation 1)3. 27(1 . multiplies both sides with J3r275Aj^, 
and traces over color, flavor, and Dirac degrees of freedom. These traces simplify considerably 
since both hard and HDL gluon propagators are diagonal in adjoint color space, [Ao,22]ab = 
5ab ^0 225 [^HDL]afe = ^ab ^HDL- '^^^ result is an integral equation for the gap function (j){K), 

9 ^HV ^\ : AHiy /T^ 



1 7. 



(3.30) 

The sum over Q runs only over relevant quark momenta, |g — /i| < Aq. Also, the 3-momentum 
k is relevant, \k — ^\ < Aq. 

3.1.3 Solution of the gap equation to sub-leading order 

In pure Coulomb gauge, both the hard gluon and the HDL propagators have the form 

A°°(P) = A^(P) , A^\P) = , A'^{P) = {S'^ -p'p^)A^{P) , (3.31) 

where A^'* are the propagators for longitudinal and transverse gluon degrees of freedom. For 
hard gluons 

A^,22(^') = (3.32a) 

KMP) = (3-32b) 
while for soft, HDL-resummed gluons 
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with the HDL self-energies [5] 



nHDL(PO,P) 



■■■^HDL 
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- m„ 
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2p ^ 
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+ 1 



Po-P 



(3.34a) 
(3.34b) 



The HDL propagators (|3.33|) have quasiparticle poles at pQ = ±u!i^t{p), and a cut between 
Pq = —p and po = P 9 . The gluon energy on the quasiparticle mass-shell is always larger than 
the gluon mass parameter, ajg^tip) > i^ig, where the equality holds for zero momentum, p = 0. 



mi = (2n+l)jiT 



©c 




(a) (b) 

Figure 3.3: (a) The contour C in Eq. (|3.35(l encloses the poles of tanh[(7o/(2r)] on the imaginary qo 
axis, (b) Deforming the contour C and adding semicircles at infinity to enclose the poles of the quark 
propagator on the real qo axis. 

I first perform the Matsubara sum, using the method of contour integration in the complex 
qo plane |9| ITB]. 

^ E fi^o) tanh (^^^ f{qo) (3.35) 

where the contour C consists of circles running around the poles u>l^ = (2n+l)7rT of tanh[(7o/((2T)] 
on the imaginary qq axis, cf. Fig. 13.31 fa). Inserting the propagators ()3.32() and 1)3. 33(1 into Eq. 
(|3.30|) . I have to compute four distinct terms. The first one arises from the exchange of static 
electric hard gluons. Since Aq 22 (-f) does not depend on po = ko — qo, only the quark propagator 
gives rise to a pole of f{qo), cf. Fig. 13.31 (b). After deforming the contour and closing it at infinity 
as shown in Fig. 13.31 (b), one employs the residue theorem to pick up the poles of the quark 
propagator. 



74 



Chapter 3. Application to color superconductivity 



with iq = eq Z{eq). Here, I have used the fact that the quark wave-function renormaUzation factor 
is an even function of its argument, ^((70) = Z{ — qQ)^ cf. Eq. (|,'-{.2(i|) . An essential assumption 
in order to derive Eq. ()3.36|) is that the gap function (t){Q) is analytic in the complex qq plane. 
This assumption will also be made in all subsequent considerations. 




Figure 3.4: Same as in Fig. 13.31 but for magnetic hard gluon exchange. Now also the gluon propagator 
has poles at fco ± p in the complex go plane. These are further away from the imaginary axis than the 
poles tq of the quark propagator, because for my choice of quark and gluon cut-offs, Aq ^ Agi, I have 
:S Aq < Agi < p. 



By the same method one computes the second term in Eq. (|3.3L)p , corresponding to magnetic 
hard gluon exchange. This is slightly more complicated, since not only the quark propagator 
but also Ag22(-P) has poles in the complex go plane. The latter are located at pQ = itp, 
i.e., go = ^0 i P ) cf. Fig. 13.41 The external quark energy Uq is fixed and, prior to analytic 
continuation fco ^ + irj to the quasiparticle mass-shell, is equal to one particular fermionic 
Matsubara frequency, cf. Fig. 13.41 The residue theorem now yields four contributions, two 
from the quark and two from the gluon poles. Using tanh[(A;o ibp)/(2T)] = ibcoth(p/2T) and 
analytically continuing /cq — ^ Cfc + I find 



tanh ( ^ 
.2T 



4e„ 



+ coth ( — ) — 
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[qolz{qo)f 
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q) 



(efc - eq + irif - p"^ 
Z'^{p + efc) (pip + eA;,q) 



(efc + €q + iriY -p2 
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2'^(p-efc)0(efc -P,q) 



{p + h + ir])^ - e^Z'^{p + ik) {p - h - ivf - e«-Z'^(p - Cfc 



. (3.37) 



Since the gluon momentum is hard, p > Agi, and thus much larger than the quasiparticle energies 
ik, iq which are at most of the order of the quark cut-off Aq <C Agi, to order 0(Aq/Agi) one may 
neglect the terms (e^ it -|- irj)'^ in the energy denominators of the first term. Furthermore, in 
the second term one may approximate Z{p it e^) ~ Z{p) = 1 + 0{g^) and (p(p it e^, q) ~ (j){p, q). 
Note that the gap function is far off-shell for p > Agi ^ Aq > — gj. Then, to order 0(Aq/Agi), 
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one may also neglect ik , in the energy denominators of the second term. I obtain 



HQ) 



[qo/Z{qo)f 



-ltanhr4^ 
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(3.38) 



As a next step one estimates to which order the two remaining terms contribute to the gap 
equation H3.3U() . At T = 0, one may set the hyperbolic functions to one. I shall also ignore the 



difference between the on-shell and off-shell gap functions, and take (t){p, q) 



const.. For the purpose of power counting, one may restrict oneself to the leading contribution 
of the Dirac traces in Eq. (|3.3fl|) . which is of order one, cf. Eqs. (|3.55|) below. In order to obtain 
the leading contribution of the first term in Eq. i)3.38() . one may also set Z^(eq) ~ 1. The integral 
over the absolute magnitude of the quark momentum is J dqq'^, while the angular integration is 
/ dcosO = / dpp/{kq). Then, the first term in Eq. 1)3. 38(1 leads to the following contribution in 
the gap equation 



k 



dq ■ 



^+1 dp 2 
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^.1 P 
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I — 

9 



(3.39) 



where I approximated k q and employed the weak-coupling solution ()1.36|) to estimate 
ln(2Aq/(/)) ~ 1/g. Furthermore, for Agi < //, the angular logarithm is ln(2/i/Agi) ~ 0(1). Ac- 
cording to the discussion presented in the introduction, the contribution from hard magnetic 
gluon exchange is thus of subleading order in the gap equation. Note that the term arising from 
hard electric gluon exchange, Eq. (|3.36|) . is of the same order as the first term in Eq. 1)3. 38(1 . and 
thus also contributes to subleading order. The way I estimated the first term on the right-hand 
side of Eq. (|3.38l) is equivalent to just taking the hard magnetic gluon propagator in the static 
limit, Aq 22(-P) — which is correct up to terms of order 0(Aq/Agj). To this order, the 

propagator for hard magnetic gluons is thus (up to a sign) identical to the one for hard electric 
gluons. Since the ratio Aq/Agi ~ ^r/u/^ = 5, this approximation introduces corrections at order 
0{g'^(f)) in the gap equation, which is beyond sub-subleading order, 0{g'^(j)). 

Similarly, I estimate the contribution of the second term in Eq. 1)3. 38() to the gap equation 

M+Aq /-fc+g dp ^ An 



2 Y 
' k 



M-Aq 



dq I 



pZ 



Agl 



(3.40) 



i.e., for my choice Aq/Agi ~ g, this term contributes beyond sub-subleading order. Note that 
this estimate is conservative, as I assumed the off-shell gap function to be of the same order 
as the gap at the Fermi surface, (/)(p, q) ~ <j). However, I know [55] that, for energies far from 

even the on-shell gap function is suppressed by one power 



the Fermi surface, £„ 



of g compared to the value of the gap at the Fermi surface, (/)(Aq,q) ~ gcj). The off-shell gap 
function at qo = p ^ Agi Aq may be even smaller. In order to decide this issue, one would 



76 



Chapter 3. Application to color superconductivity 



have to perform a computation of the gap function for arbitrary values of the energy qq, and 
not just on the quasiparticle mass-shell, qo = iq. I note that for the choice Aq ~ Agi for the 
cut-offs |1U11I1U^^ . the ratio Aq/Agi is of order one and cannot be used as a parameter to sort the 
various contributions according to their order of magnitude. The expansion of the denominators 
in powers of Aq/Agi as seen on the right-hand side of Eq. H3.38() is then inapplicable. 




Figure 3.5: Evaluating the Matsubara sum for HDL-resummed gluon propagators, (a) The original 
contour C in Eq. H3.35|l . There is no circle around the point fco = go, where the corresponding term in 
the Matsubara sum has a cut arising from the HDL gluon propagator, (b) Deforming the contour C. (c) 
The contour Ccut running around the cut. (d) The contour Ci = C +Ccut which is closed at infinity. 

The third and fourth terms in the gap equation arise from soft, HDL-resummed electric 

and magnetic gluon exchange. Evaluating the Matsubara sum via contour integration in the 
complex go plane is considerably more difficult than in the previous cases, because the HDL 
gluon propagators A^^^ do not only have poles but also cuts. The analytic structure is shown 
in Fig. 13.51 (a). Besides the poles of the quark propagator at qq = zbe^, there are also those from 
the gluon propagator at qo = ko±tLi£^t{p). The cut of the gluon propagator between —p < po < p 
translates into a cut between — p < qQ < + p. Prior to analytic continuation, the gluon 
poles and the cut are shifted away from the real axis and located at the (imaginary) external 
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Matsubara frequency /cq- 

The Matsubara sum over qo is evaluated in the standard way, cf . Eq. (|3.35|) , with the caveat 
that the contribution at qo = fco, where the cut of the gluon propagator is located, has to be 
omitted. This is similar to the zero-temperature case where the Matsubara sum becomes a 
continuous integral along the imaginary qo axis and where one has to avoid integrating over the 
cut. Alternatively, the term qo = ko can be included in the Matsubara sum if one shifts the cut 
by some small amount ibie along the imaginary qo axis. The final result will be the same, as 
one still has to circumvent the cut by a proper choice of the integration contour. 

I now deform the contour as shown in Fig. 13.51 (b). and add and subtract a contour integral 
running around the cut. Fig. 13.51 (c). The integral over the contour C + Ccut can be closed at 
infinity, yielding the contour Ci shown in Fig. 13.51 (d). One obtains 
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(3.41) 



Evaluating the integral over Ci is rather similar to the case of hard gluon exchange: one just 
picks up the poles of the quark and gluon propagators inside the contour Ci. After analytic 
continuation ko ^ ek + if] one obtains 



^£c^.o^tanh(^) A^Hk(^)|^^/^(^^)p 



tanhl^^ ^'^''^^ 



-cothl';^^ 




^nuhi^k - e<? + ir], p) 0(65, q) + A^oL(^fc + + ir], p) 0(-eg, q) 



(3.42) 



Here, I approximated the quark wave-function renormalization factor Z(u;^ j zbe^) ~ 1-1-0(5'^). I 



also expanded the denominators of the quark propagator (efcibu;^ f-|-ir/) 

■m„, I may estimate un ^ 



For my choice of the cut-off Aq < gfi 
corrections of order 0{e'^/uj'j^ 

K 
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rrir, 



.u;|,[l+0(62/a;2j]. 
> Aq > eg, i.e., the 



where p ~ and e„ 



are small everywhere except for a small region of phase space 
q. (In principle, in the expansion of the denominators there are also 



linear terms, ~ ibefc/a;^^^, but these are very small everywhere for external momenta close to the 
Fermi surface. A; ~ /i.) Note that the gap function is again off-shell at the gluon pole, although 
not as far as in the case of hard gluon exchange, cf. Eq. (|3.38|) . The residues of the HDL gluon 
propagators at the respective poles are [HI 



p2 (jp. _|_ 2,ml 



3m^ 



[uj: 



(3.43a) 
(3.43b) 
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To very good approximation, one finds that Zt{uJt,p) — l/{2ut) for all momenta p. In the 
longitudinal case, the residue is very well approximated by Zi{uji,p) ~ L0i/{2p'^) for small mo- 
menta p < rUg, while for large momenta, rUg <C p, ~ exp[—2p'^/{3mg)]/p, i.e., it is 
exponentially suppressed |128j . 

These approximate forms allow for a simple power counting of the gluon-pole contribution 
in Eq. (|,3.42|) to the gap equation (|3.3fl|) . To this end, I approximate the gap function by its 
value at the Fermi surface, (p{zizuJi^t + efcjQ) — consider the limiting case T = where 

coth[uJi^t/{2T)] = 1. Then, the contribution from the longitudinal gluon pole is 



k 



9^ T I dqq 




g^<P^^g^ct^. (3.44) 
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In the first p integral, which only runs up to the scale m^, one may approximate ~ rrig, while 
in the second p integral, which runs from rUg to Agi < p, one may take u^i ~ p. To obtain the 
right-hand side of Eq. (|3.44|) I have set k q p, and I have employed my choice Aq < gp 
for the quark cut-off. This also allowed me to approximate logarithms of Aq/rug by numbers of 
order 0(1). With this choice for the quark cut-off, the contribution ()3.44|) is of sub-subleading 
order, ~ 0{g'^(t)), to the gap equation. 

With a more careful evaluation of the integrals, one could extract the precise numerical 
prefactor of the sub-subleading contribution ()3.44|) . Note, however, that further suppression 
factors may arise from the ofF-shellness of the gap function at (j){^^e,t+^k, q), which consequently 
would render this contribution beyond sub-subleading order. As noted previously, this issue can 
only be decided if 0(go;q) is known also off the quasiparticle mass-shell, and not only on-shell. 
I also note that the 1/p^ factor in the residue Zi is an artifact of the Coulomb gauge |128j . 
and does not appear in e.g. covariant gauge. One would have to collect all other terms of sub- 
subleading order to make sure that the complete sub-subleading contribution is gauge invariant 
and the term (|3.44|) is not cancelled by some other terms. 

Similarly, I estimate the contribution from the transverse gluon pole, 



g - dqq ^ 9 (p — g (p ^ ^ g (j) , (3.45) 

kJ^i-Aq J\k-q\2uj^ Jo Jmg Ojf "^3 

where I defined i = q — p- I approximated dpp ^ dtot since, for the purpose of power counting, 
to very good approximation one may take the dispersion relation of the transverse gluon equal to 
that of a relativistic particle with mass m^, uJt{p) — {p"^ +m'^y/'^ . I also used Aq < mg <C Agi and 
A; ~ g ~ /i. In conclusion, also the transverse gluon pole possibly contributes to sub-subleading 
order in the gap equation, with the same caveats concerning the off-shellness of the gap function 
as mentioned previously. 

Let us now focus on the integral around the cut of the gluon propagator in Eq. (|3.41j) . I 
substitute qo hy po = ko — qo = to and use the fact that tanh[go/(2r)] = — coth[w/(2T)]. Since 
the gluon propagator is the only part of the integrand which is discontinuous across the cut, I 
obtain after analytic continuation /cq ^ + ir] 

rf^OT^tanh — AjIdl(P) 



27rz Jc... '^2 \2TJ ^ [qo / Z {qoW - e^^ 

du - coth — — — — —2 />cut ^> P) > 3.46 
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where pl'^^{uj,p) = lmA^-^-^{u! + irj,p)/7r is the spectral density of the HDL propagator arising 
from the cut. Exphcitly, 
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(3.47b) 



In order to power count the contribution from the cut of A|jj-)l to the gap equation, it is sufficient 
to approximate the spectral density by |55] 
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(3.48) 



This form reproduces the correct behavior for uj <^ p. For uj < p, it overestimates the spectral 
density when p < rUg, while it slightly underestimates it for p > rUg. For the gap equation, 
however, this region is unimportant, since the respective contribution is suppressed by the large 
energy denominator (e^ — a; + irj)'^ — [Z{ek — uj) e^]^ ~ p'^ in Eq. 1)3. 46(1 . To leading order, one 
may set Z{ek — w) ~ 1. I also approximate (j){ek — w,q) ~ <j). Then, the uj integral can be 
performed analytically. (One may compute this integral with the principal value prescription; 
the contribution from the complex pole contributes to the imaginary part of the gap function, 
which I neglect throughout this computation to subleading order, cf. Sec. 13.21 ) This produces 
at most logarithmic singularities, which are integrable. I therefore simply approximate the uj 
integral by a number of order 0(1). Consequently, the contribution from Eq. 1)3. 46(1 to the gap 
equation is of order 
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(3.49) 



where I approximated the p integral by a method similar to the one employed in Eq. 1)3. 44() . For 
my choice Aq < (7|U, Eq. (|3.49)) constitutes another (potential) contribution of sub-subleading 
order to the gap equation. 

Finally, I estimate the contribution from the cut of the transverse gluon propagator. For 
all momenta p and energies —p < uj < p, a very good approximation for the spectral density 
(|3.47b|) is given by the formula 
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(3.50) 



This approximate result constitutes an upper bound for the full result H3.47b|) . The advantage of 
using this approximate form is that, interchanging the order of the p and uj integration in the gap 
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equation, the former may immediately be performed. Approximating Z(€k — w) ~ 1, neglecting 



the dependence of the gap function on the direction of q, and defining A = max(|A; 
T = the contribution to the gap equation is 
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(3.51) 



Here, I have used the fact that the particular combination of arctan's in the first line effectively 
cuts off the UJ integral at the scale u) ~ M. As usual, I have set k q ^. After estimating 
(/>(efc lb a;, Q') ~ one first performs the integral over u; 
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(3.52) 



The denominator under the logarithm in Eq. (|3.52|) does not contribute because of the sum over 
(72. In the region S, > (p one may neglect ik in the numerator under the logarithm and obtains 
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where the estimate on the r.h.s results after expanding the logarithm in powers of Cg/M. In the 
region ^ < (j) one has 
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(3.54) 



Consequently, the contribution 1)3. 51(1 is at most of sub-subleading order. 

The remaining term from the evaluation of the Matsubara sum in Eq. ()3.41|) is the contri- 
bution from the quark pole, i.e., the first line of Eq. 1)3.42(1 . This has to be combined with the 
subleading-order terms from hard-gluon exchange, i.e., from Eq. ()3.36|) and from the first line of 
Eq. ()3.38|) , in order to obtain the gap equation which contains all contributions of leading and 
subleading order. Before doing so, however, I also evaluate the Dirac traces in Eq. ()3.30j) . In 
pure Coulomb gauge, I only require 
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where I used p^ = (k — q)^ = k"^ + cp' — 2k q\^ - c^io eliminate k • q in favor of p^. Let us estimate 
the order of magnitude of the terms arising from the traces at the Fermi surface, k = p. Setting 
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fi + S,, where — Aq < ^ < Aq, one obtains 

Tr, (A+7oAq7o) = 2-^+o(^), (3.56a) 




(5U_^y)Tr,(^A+7,Aq7,j = " ^ + ^ (^^ j • (3.56b) 

As shown above, the contribution from hard-gluon exchange is at most of subleading order. 
Thus, for this contribution it is sufficient to keep only the leading terms in Eq. H3.56() . i.e., 
one may safely neglect terms of order 0(^^//i^) < 0(Aq/Agj) ~ 0{g^) or higher. Note that, 
since for hard gluon exchange p ~ ^ Agi, the terms /{2kq) cannot be omitted. However, 
since the magnetic gluon propagator is effectively ~ cf. Eq. ()3.38|) . i.e., (up to a sign) 

identical to the electric propagator, these terms will ultimately cancel between the electric and 
the magnetic contribution. This cancellation is well known, see for instance Ref. j76j . and is 
special to the spin-zero case. It does not occur in spin-one color superconductors where there is 
an additional exponential prefactor which suppresses the magnitude of the spin-one gap relative 
to the spin-zero case 

As is well known, electric soft-gluon exchange also contributes to subleading order in the 
gap equation. Thus, as in the case of hard-gluon exchange, one may drop the terms of order 
0(^^//i^) in Eq. (|3.56a|) . On the other hand, magnetic soft-gluon exchange constitutes the 
leading order contribution to the gap equation. I therefore would have to keep all terms up to 
subleading order, i.e., ~ 0{^/ n). Fortunately, the corrections to the result ()3.56b|l are of order 
0(^^//i^) ~ 0{g^), i.e., they are of sub-subleading order and thus can also be omitted. 

I combine Eqs. (lOHl) . (|S3HI), and the first hue of Eq. (|3l^ . and assume that the gap function 
is even in its energy argument, (/>(— eq,q) = (p{eg,q}, and isotropic, (p{eg,q} = (p{eg,q) = (pg. 
Then, on the quasiparticle mass-shell /cq = Cfc the gap equation (|3.3fl|) becomes 
^2 /-M+Aq I-. y'^if \ / f \ rk+q ( A 

= Lk k 4^ (^) 0, / {e(p - A,) ^ + e(A, - 
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(3.57) 

The next step is to divide the integration region in the p — q plane into two parts, separated by 
the gluon "light cone" |efc — seq\ = p. For my choice Aq ^ Agi the region, where \tk — seg\ < p, is 
very large, while its complement is rather small. In order to estimate the contribution from the 
latter to the gap equation, one may approximate the HDL gluon propagators by their limiting 
forms for large gluon energies, cf. Eqs. (|3.33j) . (|3.34j) . 

po » P : AfjDL(^) ^ , ^HDLiP) ^ \ ■ (3.58) 

p^ 

Following the power-counting scheme employed previously, the contribution from the electric 
sector is of order 
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This is a contribution of sub-subleading order, as long as one adheres to the choice Aq < gfi. 
Analogously, I estimate the contribution from the magnetic sector to be 
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(3.60) 



Consequently, all contributions from the region \ek — seg\ > p are of sub-subleading order, and 
the further analysis can be restricted to the region {ik — seq\ < p. In this region, it is permissible 
to use the low-energy limit of the HDL gluon propagator, which follows from Eqs. (|3.33|) . ()3.34|) 
keeping only the leading terms in the gluon energy, 
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Here, I only retained the real part of the transverse gluon propagator, since the imaginary part 
contributes to the imaginary part of the gap function, which is usually ignored. (In Ref. |55j it 
was argued that, at least close to the Fermi surface, the contribution of the imaginary part is 
of sub-subleading order in the gap equation.) With the approximation (|3.6H) . the gap equation 
1)3. 57() becomes 
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where I already performed the integration over hard gluon momenta p > Agi. The integration 
over soft gluon momenta can also be performed analytically. Formally, the terms ~ p^ /{Akq) 
give rise to subleading-order contributions, ~ Agj/(8fcg), but they ultimately cancel, since they 
come with different signs in the electric and the magnetic part. Other contributions from these 
terms are at most of sub-subleading order. Exploiting the symmetry of the integrand around 
the Fermi surface and setting A; ~ /i, I arrive at 




A6 



?2 _ ;2| 



(3.63) 



Here, I have neglected terms ~ — sig against 3m^ under the logarithm arising from soft 
electric gluons, and terms ~ (e^ — SEg)^ against M^(efc — seg)^ under the logarithm from soft 
magnetic gluons. 

Now observe that the gluon cut-off Agi cancels in the final result. 
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where b = 256TT^[2/{Nfg^)]^/^. This is Eq. (19) of Ref. [HHI, since f = g"^ /{IStt'^), with the 
upper hmit of the {q — /i) integration, 6, replaced by the quark cut-off Aq. 

The solution of the gap equation ()3.64|) is well known, and given by Eq. (|1.36|) . As was 
shown in Ref. ^Ei, the dependence on the cut-off Aq enters only at sub-subleading order, i.e., it 
constitutes an 0{g) correction to the prefactor in Eq. H1.36() . Therefore, to subleading order I 
do not need a matching calculation to eliminate Aq. 

The result 1)3.64(1 shows that the standard gap equation of QCD can be obtained from the 
effective action p.46|) . The above, rather elaborate derivation of Eq. ()3.64|) demonstrates that, 
in order to obtain this result, it is mandatory to choose Aq <ti Agi. This also enabled me to 
identify potential sub-subleading order contributions. However, I argued that, at this order, the 
off-shell behavior of the gap function has to be taken into account. This requires a complex 
ansatz for the gap function, which will be considered in the following Sec. 13.21 
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3.2 The imaginary part of the gap function 

The scalar gap function (^{K) contains the energy and momentum dependence of the (2SC) gap 
matrix ^^{K) = J^T2'y^ Aj^ ©(Aq — \k — ^|) (piK), which solves the gap equation Eq. (|3.3U|) 

<^+{K) = 5' ^ E ^ab(K - Q) l,{T^f H+(Q) . (3.65) 

Q 

Due to the dependence of the gluon propagator IS.{K — Q) on the external energy ko = —tOn = 
— i(2n+l)7rT the solution (j){K) must be energy-dependent itself. Hence, solving the gap equation 
self- consistently requires an energy dependent ansatz for the gap function. 

The function (t){K) solving Eq. ()3.65|) for some Matsubara frequency ko can be considered as 
the analytic continuation of the gap function (/)(a;,k) of physical and therefore real frequencies 
cj. Then the physical gap function 0(0-", k) has to have singularities along the real w-axis so that 
the analytic continuation 4'{K) is, as required, not a constant in complex energy /cq- To see this 
write without loss of generality 

(t>{K) = 4>{K) + ct>oiy) , (3.66) 

where the energy dependence of (j){K) is contained in (p{K). It is physical to assume, and it 
will be confirmed when solving the gap equation, that (piK) becomes local, i.e. independent of 
energy, for asymptotically large energies, i.e. that (^{K) — > and (f){K) </'o(k) for ko oo. 
Cauchy's theorem gives for any ko off the real axis 

A/u v\ I dz 0(2;, k) 



27ri z — ko 

c 

00 

du]' (t){uj' + ir], k) - (/)(a; - irj, k) ^ f (/"(^^k) 

27ri Lo' — ko I 2iTi z ' 

-00 c 

where the contour C circumvents all non-analyticities of </>(i^) on the real /cq— axis as depicted 
in Fig. 13.61 and C is a circle whose radius tends to infinity. I define the spectral density of (j){K) 
by 

p^{uj,k) = ^ [(t){uj + irj,'k) - (j){uj -ir],'k)] , (3.68) 
zvrz 

which contains all non-analyticities of (/>(a;,k). Using (j){z,'k)\zec' = 0o(k) one obtains 

00 

0(A:o, k) = f du ^^^^^ + ,/.o(k) , (3.69) 
J id — ko 
—00 

where the first term is identified as 4>{K). Eq. (|3.69|) holds for all ko surrounded by contour C, 
in particular for frequencies infinitesimally close to the real axis, ko = iv + ie with e > r/ > 0. 
For such frequencies one can employ the Dirac identity l/{x — ie) = V{l/x) + iT:5{x) to obtain 

4>{uj + ie, k) = 

= I du;' P^^'^'^^l +M^)=V f dw' + </'o(k) + iT^p^iu, k) . (3.70) 

J u) — ui — le J uj — u 
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Figure 3.6: The contour C circumvents all non-analyticities of 4i{K) on the real /cq— axis in the shaded 
area. The dots indicate the fermionic Matsubara frequencies, fco — —i(2n+ 1)ttT. 



Hence, for real (j)o(k) and p^{uj,'k) 

oo 

Re^{uj + ie,k) = V I ^^^^^^ + 0o (k) , (3.71) 

J UJ — Uj' 

lra.(j){oj + ie,'k) = 7r/)0(u;,k). (3.72) 
Combining Eq. 1)3. 68|) with Eq. ()3.72|) yields the identity 

Im Mijj + ie, k) = — \4>{'^ + k) — Mu) — ie, k)] , (3.73) 

i.e. non-analyticities along the real /co-axis are equivalent to a nonzero imaginary part of the 
gap function. Furthermore, one finds the dispersion relations for (f){uj + iry, k) 



oo 



Re(j){uj + ie,k) = -V dJ ^ (3.74a) 



-oo 

oo 



Im^(c. + ze,k) = -Ip /d^'^^^i^^^l±^, (3.74b) 

TT J U>' — UJ 

— oo 

i.e. lm(p{u; + ie,k) and Ke(p{iJ + ie,k) are Hilbert transforms of each other. Ret/) = 7i\lm(j)]. It 
follows that (j) ^ only if both Re (p and Im (p are nonzero. Consequently, the gap function (p 
is energy-dependent only if it has a nonzero imaginary part, Im0 = Imi;^ ^ 0, which in turn is 
generated by its non-analyticities along the real ko axis, cf. Eq. ()3.73|) . 
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This investigation accounts for the energy dependence of the gap function by including the 
non-analyticities of cj) in the solution of the gap equation. This finally leads to a complex gap 
equation, whose real and imaginary parts are coupled integral equations 

Re0 = Re0[0, A,S] , (3.75a) 
Im (t) = lm(j) [(/), A, S] , (3.75b) 

which have to be solved self-consistently. Of course, there are analogous DSEs for A and S, which 
in principle would have to be solved self-consistently simultaneously with Eqs. (|3.75|) . However, 
it has been shown jl26j that to subleading order one may replace A = Ahdl in Eq. H3.75a() . (This 
can be refined by taking A = Ahdl for soft gluon momenta and the free propagator A = Ao,22 
for hard momenta [IHIO Also S may be replaced by its one-loop approximation in Eq. H3.75a|) 
to subleading order |1251 198j . Furthermore, the contribution of ImS to Re</) in Eq. (|3.75a)) has 
been shown to be of sub-subleading order jl25j and therefore will be neglected in the following. 
However, the contribution of Imcj) to Hecp in Eq. ()3.75a() has not been systematically accounted 
for, yet. In the following it will be checked if lm(j) contributes leading or subleading order 
corrections to the real part of the gap function or not, i.e. whether the common simplification 
of Eq. (irzHa 

Re(/) = Re(/)[Re(?i, A,S] (3.76) 

is justified to subleading order or not. To this end it is admissible to make the same approxi- 
mations in Eq. p.75b|) as in Eq. H3.75a|) yielding 

Re ^ = Re (/>[(/), Ahdl, Ao,22, Re Si loop] , (3.77a) 
Im(/) = Im0[(/), Ahdl, Ao,22,ReSiioop] • (3.77b) 

This can be justified when considering Imcj) = luKp as the Hilbert transform of Re0 in its 
subleading order approximation. As such it must be insensitive to any sub-subleading order 
corrections to Re0. In other words, one assumes that contributions, which do not enter Re(/> 
directly via Eq. (|3.75aj) do not enter Kecj) indirectly via Eq. (|3.75b|) . As stated in the beginning 
the main contributions to Im0 are expected to arise from the energy dependence of the gluon 
propagator A{K — Q). 

The gap in the quasiparticle excitation spectrum is given by the modulus of the complex gap 
function 

|(/.| = y^(Re(/>)2 + (Im0)2 . (3.78) 

Hence, there are two possibilities how the physical quantity \<j)\ is aff'ected by Im(j) to subleading 
order. The first is that Imcj) enters Re(/> to subleading order through Eq. (|3.77aj) as already 
discussed. The second is that luKp appearing on the r.h.s. of Eq. (|3.78|) is of order Hecj). In 
|55j it is estimated that Im (p = on the Fermi surface and Im </> ~ 5 Re (j) exponentially close to 
the Fermi surface. This will be confirmed in the following. It follows that in these momentum 
regimes Imcj) contributes at most at order g'^Recp to \<p\, i.e. at sub-subleading order. Hence, 
the second of the two possibilities mentioned above can be discarded in this momentum regime. 
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However, at momenta |A; — /u| ~ (7 /u it is lm(j) ~ Kecj) and both the real and the imaginary part 
contribute at the same order of magnitude to \4>\. Also this will be confirmed in the following. 
The modulus appears in the complex gap equation only in the on-shell energy dispersion relation 



V'(fc-^)2 + |</,|2. 



(3.79) 



In the region = \k — fi\ ^ g fj,, where lm(p ~ Re0, it is ^ and therefore Cq ~ C- Hence, in 
the gap equation it is self-consistent to subleading-order to approximate |(/>| ~ Re(/) and write 



^(A:-/i)2 + (Re(/>g)2. 



(3.80) 



Consequently, in contrast to \(j)\ the on-shell energy can be affected by lui<j)q only through 
Re (j)q for all C- 

Before actually solving the complex gap equation it is worthwhile investigate how the sym- 
metry properties of the gap matrix ^~^(K) affect the gap function 4){K). Remembering ipc{K) = 
C-ip'^{—K), il^c{K) = ip'^{—K)C, and C = —C~^ = —C'^ it follows from the antisymmetry of 
the quark fields that 



K 



K 



K 



Hence, the gap matrix must fulfill 



(3.81) 



(3.82) 



Since C^^K^C ^ = [75A]'^j^]-^ and in the 2SC case [J3T2]"^ = JzT2 it follows for the gap function 



0(i^) 



-K). 



Assuming that the gap function is symmetric under reflection of 3-momentum k, 
fco, — k) one obtains with Eqs. 1)3.7113. 72|) 



Re 4>{ijj + ir/, k) 
Im (j)[u) + irj, k) 



Re (j){—oj + irj., k) , 
— Im (p{—uj + irj, k) 
-p^(-u;,k) . 



(3.83) 
•o,k) = 

(3.84a) 
(3.84b) 
(3.84c) 



Hence, one found that Re (j) is an even function in iv while Im (/> and p,^ are odd. 

Furthermore, one introduces the spectral representations for the gluon propagators. In the 
case of soft, HDL-resummed electric and magnetic gluon propagators, cf. Eqs. ()3.33)) . one has 



pPf (^, p) {6 [io - coeAp)] -H^- ^iAp)]} + Ptfi^^ P) (^(P -H)^ (3-85) 
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where pff are the contributions from the cut of the logarithm in the HDL self-energies H^^^, 
cf. (|3.,'Mj) . and already given in Eqs. (|3.47|) . The contributions from the gluon poles at a; = iuj£^t 
are given by the respective residues of the propagator, p^'f^{oj,p) = Zi^t{^,p), cf. Eqs. (|H.43j) . 
For the respective hard gluon propagators, cf. Eqs. ()3.32|) . one has Aq22(.^ ^ ^^iP) = 
and Ag 22(1^ i irj, p) = —l/[{u> ± irj)^ — p^] and therefore with Po 22('^' P) = [^0 22('^ + P) ~ 
Ao:22(^-^»?>P)]/(27ri) 

PoM^^P") = 0' (3.86a) 
PoM^^P") = sign(a;)5(u;2 -/) . (3.86b) 

With these spectral densitis, /Ohdl ^^'^ Po 22' ^^-y express the respective gluon propagators 
in their spectral representations 

00 „ 00 

A^{P) = -\+ f du^S^^ , A\P)= [ du.'-^. (3.87) 

P'^ J OJ — Pq J OJ — Pq 



3.2.1 Solving the complex gap equation 

In the following solution of the complex gap equation 1)3. 77j) it will be shown that its imaginary 
part, Eq. (I3.77bl) . can be decomposed as 

lm(t) = A + B, (3.88) 

where A contains only contributions, where the gap function is on the quasi-particle mass-shell, 
and B contains the rest. Then all contributions in the imaginary part of the gap function, 
which are generated by the energy dependence, i.e. by the non-analyticities of (p, are collected 
in B. Therefore, all terms in B have not been considered so far and neglecting B must lead 
to the known subleading order result for Yie(j){eq,q). In the following it will checked, whether 
this result is self-consistent. Self-consistency is fulfilled only if the contributions contained in 
B enter Yie(j){eq,q) beyond subeading order, i.e. at sub-subleading order. To this end it is first 
assumed that self-consistency is indeed fulfilled and then analyzed, if B really contributes beyond 
subleading order. This procedure comprises of the following steps. 

In order to obtain A and first the Matsubara sum over the internal frequency is per- 
fomed, cf. Sec. 13. 2. Tl and then all imaginary contributions are extracted and organised according 
to Eq. (|3.88|) . cf. Sec. 13.2.11 In the first part of Sec. 13.2.21 the order of magnitude of Imc/) will 
be determined by inserting the known leading order result for Re i;A(eq, q) into A, while the yet 
unknown term B has to be neglected first. Of course, this estimate for Im (p can be self-consistent 
only if reinserting it into B in Eq. (|3.88|) only yields sufficiently small corrections to Im (j), i.e. if 

lxnct)K,A. (3.89) 

This will be analyzed in the second part of Sec. 13.2.21 where this estimate for Im (j) is used to 
determine the order of magnitude of B. It will be shown that Eq. (|3.89|) is indeed fulfilled, 
i.e. that B constitutes a small correction to A. With these results the order of magnitude of 
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Re0(ek,k) is estimated in Sec. 13.2.31 via the dispersion relation Eq. H3.74a|) . i.e. by Hilbert 
transforming A and B. It is found that Re(/>(ek,k) ~ 7i[A] ~ (p for momenta k close to the 
Fermi surface, while Ti.[B] ~ g'^cp, i.e. Imc/) contributes beyond subleading order to Re (p{e\i, k) or 
according to Eq. (|3.76j) 

Re = Re ^ [Re . (3.90) 

In Sec. 13.2.31 it is shown that also (/>o(k) ~ and that the imaginary part of the gap function 
contributes to (/)o(k) at g'^cp, i.e. beyond subleading order. Consequently, 

<Po = M^(i<P] (3.91) 

to subleading order. Hereby it is shown that to subleading order Eq. (|3.77aj) becomes Eq. (|3.76j) . 
which means that the imaginary part of the gap function does not enter the real part at this 
accuracy. In other words, the real part of the gap equation is decoupled from its imaginary part. 
On the other hand, the imaginary part of the gap equation is not decoupled from the real part 
of the gap equation, cf. Eq. (|3.89j) . It turns out that for momenta close to the Fermi surface 
Im (j) can be calculated self-consistently from Re (f) alone, which is performed in Sec. 13. 2. 51 



Performing the Matsubara sum 

The complex gap equation (|3.65j) reads to subleading order, cf. Eq. (|3.41j) . 

where the quark wave-function renormalization factor Z{qQ) = [1 + 5^ ln(M^/gg)]~^ contains the 
effect of the normal quark self-energy ^(Q). As explained before, only its real part contributes 
at subleading order and therefore will be accounted for in the following, while its imaginary part 
will be neglected. Physically this is equivalent to assume the (normal) quasi-particle excitations 
to be infinitely long-lived, i.e. ignoring their decay due to scattering processes with quarks 
inside the Fermi sea. Mathematically it amounts to neglecting the cut of the logarithm when 
performing the Matsubara sum over qq. Consequently, at subleading order the only effect of 
Z((7o) in Eq. (|3.92|) is a shift of the poles of the propagator 

to go = ±^(eq) eq = ±eq- 

The Matsubara sum over go in E<1- (|3.65j) is performed via contour integration, cf. Fig. 13.71 

A^^-*(fco,p,q) ^ T J2 ^''\Q - K)AiQmQ) 
c 
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Figure 3.7: The contour C in Eq. H3.94(l encloses the poles of ta,nh[qo / {2T)] on the imaginary qq axis. 
The additional poles and the cut at qq — fco arise from the the gluon propagator, while the two poles on 
the real axis are due to the quasiquarks. The undetermined singularities of the gap function on the real 
go— axis are indicated by the shaded area. 



In order to introduce the spectral densities of the gap function and of the magnetic and longitu- 
dinal gluon propagators, the contour C in Eq. (|H.94|) is deformed as shown in Fig. \'A.H\ yielding 



M''\ko, p, q) = I P: JC''\qo) + I ^ \lC'^q, + - }C''\qo - zr?) 



2TTi 



2TTi 



ago 
2-Ki 



IC'^\qo + ko + irj) - IC''\qo + ko - ir,) 



e,tt 



= loo + Iq + Iko 

The explicit form of the integral /qo reads 



^^tanh(^)A^'*(Q-K)A(QmQ). 



(3.95) 



(3.96) 



For the integration along the contour Coo one parametrizes dqo = i\qQ\e^^d9 and find for 
oo 



tanh ( 

2T 



A\Q - K) 



1 



1 , A(Q) ^ ^ , 0(Q) 



(3.97) 
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qo=k( 




Figure 3.8: Deforming tlic contour C to introduce tiie spectral densities of A ' and <j), cf. Eq 



Consequently, 



dqolC\qo) ~ — , dqoK,\qo) ~ \ 



(3.98) 



and loo vanishes, /qo = . Using Eq. ()3.93|) the integral /q can be decomposed 

oo 



[A(go + iv, q)0(^o + ir], q) - ^(^o - ^ij^Qo - iv, q 

oo 

1 _ f dan I 



'0 + qj - 9>(go - q 



1 

'2I 



-00 



E 



'0 + i??, q) 



q <7=± 



27ri 2 



o-eq+»7 



- qj 



(3.99) 



where the notation V±i^ refers to the principal value prescription for the quark pole at qo = ztcq. 
In the first term on the r.h.s. of Eq. 1)3. 99() the difference of the retarded and advanced gap 
function may be replaced by the spectral density of the gap function, cf. Eq. (|3.68|) . In the 
second term one may substitute q'o = qo — creq and deform the contour of integration into an 
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infinitesimally small semicircle around the quark pole at = 0. With dqQ = irje^^dO this term 
becomes 

^ E ^ / I tanh (^^) A''\ae^ - k,, v)Z\e^) 4>{a~e^ + ir/, q) 

= E A^'*(aeq- A:o,p)^2(eq)(/>(cTeq + ir?,q) . (3.100) 



The third term in Eq. (|3.99() can be evaluated analogously to the second. The result will differ 
from the second term only by containing the advanced instead of the retarded gap function. 
Collecting all three terms one obtains using </>(a' + irj) + (j)[ijO — irf) = 2Re0(ti; + irf) 

lo = ^E-^^e-, / d,o|tanhfe)A^'*(go-feo,p)^^(^o)^^^ 

— oo 

'^h'"''^ (It) E ^''*(^eq - A:o,p) . (s.ioi) 

The first term arises from the non-analyticities of the gap function, (j){Q), and the second from 
the poles of the quark propagator A(Q), cf. Eq. (|3.93() . at qq = iteq. The last integral I^q reads 

ho = 

oo 
— oo 

+ / ^^'°*^'(^^) ^'''ilo + iV,P)MQo + ko,ci)^iqo + ko,cO 

-11 

+ / ^^coth(^^^)A^'*(go-ir?,p)A(go + A:o,q)</'((Zo + A;o,q). (3.102) 

V 

In the first term the difference of the retarded and advanced gluon propagator, A^'*(go ± iv^ P); 
may be replaced by the gluon spectral density. In the second and third term the integrals over 
qo can be performed after taking the limit coth(go/2r)/2 — > T/qo for ry — > and parametrize 
the contour as dqo = irje^^dO. One obtains for the second term 

.„ 

/?' n p df) 1 T 
4— ^ A^'*(0 + iv, p)A{ko,q)Hko,q) = --A''\0 + iTj, p)A{ko, q)</>(A;o, q) . 

TT 

(3.103) 

For the third term the result differs only in the occurence of the advanced gluon propagator 
instead of the retarded. Combining all three terms in Eq. (|3.1()2|) yields 



'fco = T^E^^O / dgo^coth(^|;)</.(go + A:o,q)Z2(go + M 



2eq ^ J 2 \2TJ qo-aici + ko 
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-rReA^'*(O + i?7,p)A(fco,q)0(A;o,q) • (3.104) 

The first term is due to the non-analyticities of the longitudinal and magnetic gluon propagators, 
A^'*. The second arises from the pole of coth(qo/2T) at go = corresponding to the large 
occupation number density of gluons in the classical limit, qo <^T. The latter contribution has 
been shown to be beyond subleading order and will be discarded in the following. 



The imaginary part of the gap equation 

After having performed the Matsubara sum one considers energies /cq close to the real axis, 
ko = u + ir] in order to employ the Dirac identity to split the complex gap equation (|3.92|) into 
its real and imaginary part. To this end one extracts the imaginary parts of the contributions 
Iq and cf. Eqs. 1)3.10113. 104(1 . For the imaginary part of Eq. (|3.92|) one finds using that Re0 
and are even in on 



3 



Im (f)(uj + irj, k) 



TTs {k, P, q) ImM {uj + ir], p, q) + Tr* {k, p, q) ImM + iri,p, q] 



where 



(27r)3 
lmM^'^{uj + ir], p, q 



(3.105) 



dqo 



90 



Z [qo 



(Tec 



tanh ( 1 + coth 
2T 



UJ 



IT 



Re (A(eq, q) Z^{i^) ^ cjp^'*(a; - ae^, p) 



cr=± 



tanh ( ^ 
2T 



+ coth 



UJ 



2T 



2T 



lmM^^\uj + ir], p, q) + lmM'^J^{uj + irj, p, q 



(3.106) 



Note, that in the above Eq. I|3. 105(1 uj + irj appears with a minus sign in the argument of the 
gluon propagators. The traces over Dirac space Trf'*(A;,p, g) are introduced in Eqs. ()3.55() The 
first term on the r.h.s of Eq. ((3.106(1 . ImA^^*, is due to the non-analyticities of 4>{K) aiid has 
been neglected in all previous treatments. Inserted into Eq. ((3.105(1 it constitutes the term 
;S[Imi;^] introduced in Eq. ((3.88(1 . The second term, ImA4^*, contains all contributions to (t){K) 
that have been considered so far to subleading order. Inserted into Eq. ((3.105(1 it constitutes the 
term ^[Rec/)]. As already announced it contains the real part of the gap function always on the 
quasiparticle mass-shell. It is instructive to shortly analyze this issue before proceeding further. 



The anomalous propagator E{Q) = Z'^ {qo) (j){Q) / {q^ 

/9H(w,q) = ^ [E{uj + ir], q) -E{uj -ir],q)] 
Z'^{uj + ir]) (f){uj + ir], q) 




has the spectral density 



Z'^{uj — ir]) (p{uj — ir], q) 



[Z(cj)eq]2 + sign{uj)ir] 



UJ" 



[Z{uj)e^Y 



sign(w)i?7 



Z\uj)V 



[uj + ir], q) 



[UJ 



ir],q) 



UJ" 



[Z{uj)e^Y 
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-i7rsign(a;) ^^(eq) 5 (w^ - [(piu + ii], q) + (j){uj - ii], q)] | 



(3.107) 



where the cut of Z{uj) has been neglected. It follows that /9^(a;, q) / leads to additional support 
of /9H(a;,q) around the quasiparticle pole oj = Iq. Furthermore, neglecting Imc/) is equivalent to 
approximating 

ps(a;,q) ~ -sign(u;) Re</.(eq + i?7,q) Z^{e^) 5(00^ - e^) (3.108) 

as it has been done in [SS], cf. Eq. (41). As a consequence, the gap function on the r.h.s. of the 
gap equation is always forced onto the quasiparticle mass-shell and the gap equation takes the 
standard form, cf. Eq. ()3.64() . The occurrence of the external energy on the r.h.s. due to the 
energy-dependent gluon propagators indicates that the solution still would possess some energy 
dependence although not provided in the ansatz Eq. ()3. 108(1 . This demonstrates explicitly the 
inherent inconsistency of this ansatz. 

It is interesting to note that in Eq. (|3.1U7|) Im 7^ does not shift the quasiparticle pole into 
complex plane, since tq = \/{q — ^)^ + \4>\''^- Therefore, no damping of H is caused. In order 
to restore the effect of damping due to lm.(j) one would have to include ImS, which is beyond 
the scope of this work. 

In the limit of small temperatures, T ^ 0, the hyperbolic functions in Eq. (|3.1U6|1 simplify, 
yielding for u; > 



Z^(eq) Re (/>(eq, q) /•*(w - gq, p) i 



ImA4g*j.^o(a; + ir], p, q) + ImA4^*j,^o(w + ir/, p, q) 



(3.109) 



Here, use was made of the oddness of in cf. Eq. (|3.84c() . 
3.2.2 Estimating the order of magnitude of A and B 

Following the strategy explained in Sec. 13.2.11 it is necessary to know the order of magnitude of 
the real part of the on-shell gap function Re(/)(eq,q) in order to estimate the magnitude of the 
parts A and B of Im cj), cf. Eq. H3.88|l . It is known that 0^ = Re (t){^VL-, k) is momentum dependent 

4>{x) ^ 4>f'' F{x) , (3.110) 
where </>o^^ is the value of the gap for momenta on the Fermi surface 



6fc = 266'o/iexp(-^) (3.111) 
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with g = g/{^\/^7r) and the constants (subleading and therefore in principle irrelevant for the 
present purposes) 



5/2 



b = 256 TT^ 



Nfg^ 



6'n 



7t' + 4 



'0 = exp 



The variable a; is a exponential measure for the distance of k from the Fermi surface 



X = g In 



2bfi 



k - n + ek 



(3.112) 



(3.113) 



For \k — fi\ ^ 4> one finds x = tx 12 + O(^), while for \k — [i\ ~ M it is x ~ 0{g). The function 
F{x) is to leading order given by F{x) = sin(x). Hence, for momenta exponentially close to the 
Fermi surface, |A; — /i| ~ the on-shell gap function is ^{e\^^ k) = 
to leading order. For momenta \k — [x\ ~ M, one has 0(ek, k) = g 
the on-shell gap function is sharply peaked around the Fermi surface. This exponential decay 
off the Fermi surface becomes apparent when considering intermediate distances from the Fermi 
surface (between 4> and M) defined by the variable scale 



■^^^[1 + 0(f)], i.e. constant 
'q^*-'. Hence, the real part of 



Ay = ^M^-y 



(3.114) 



with < y < 1. For |A; — ^| ~ Ay one has x = y7r/2 + 0{g) and F{y) ~ sin(y7r/2) ~ y. Hence, 
decreasing y (i.e. exponentially admixing the scale M to Ay and exponentially receding from the 
Fermi surface) decreases the magnitude of the gap approximately linearly. For y ^ g the gap 
function has decreased to g4>W''~-'. The intermediate scales may be also used to investigate the 



generation of the so-called BCS-log. Denoting Ai 
one may write 



' and Ai 







M and exploiting Ai <C Ag 



< An 



M 



9 



Ai 

2 

J eq 





+ 9' 



di 
en 



+ 9 



^0 
2 jdi 

en 



g^\n{V2 + l) ^o+g^ 



In 



9^ In 



+ g-^ In 



2A| 
Ai 



ln(V2 + l) 



+ 9 In 




(3.115) 



where use was made of Eq. (|3.111|) yielding the BCS-log, \n{M/(j)) ~ ^/g- It is shown that 
the BCS-log arises from integrating over intermediate scales, Ai < ^ < Ag. In the region 
Kg < ^ < M the BCS-log does not occur due to the exponential g under the logarithm. In 
addition the contribution from this region is suppressed by an extra factor g from (pq ~ gcjyQ. 
In the QCD gap-equation the gluon propagator has to be added to the integrand. (Then the 
region Ai < ^ < A^ is enhanced by an additional large logarithm due to almost static. Landau- 
damped magnetic gluons.) However, the above observation that the BCS-log is generated by 
intermediate scales Ai>y>g remains valid in the full gap equation. 
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Landau-damped gluons contributing to A 



For the purpose of power counting the various terms in Eq. H3.105() one may restrict oneself to 
the leading contribution of the Dirac traces Eqs. (|3.55|) . which is of order one. The integral 
over the absolute magnitude of the quark momentum is / dqq^, while the angular integration 
is / dcosO = J dpp/{kq). Furthermore, only the zero-temperature limit, T = is considered. 



allowing for the estimate .^^(cq) ~ 1. The contribution of M.J^j<^q in Eq. (|3.1fl9j) arising from 
Pcut(^ - eq,p) to ^ is 



^cut(W:k) 



5 

2 / 



Re0(eq,q) / dpp p^^^{uj* ,y> 



(3.116) 



where a;* = — eq < oj, 5 = min(a;, Aq) and A = max(|,^ — Cl)^*) with C, = \k — p\. Due to the 
condition \ < p < Agi it immediately follows that ^^ut = for u; > Agi + Aq ~ /x. Inserting the 
approximative forms 



M2 



Pcut(w*>P) - — 



UJ p 



vr p^ + lyW^u*)"^ 



1 



2M^ UJ* 



TT p (p2 + 3 77),2)2 



(3.117) 



the integration over p can be performed analytically. For energies to < Agi one finds for the 
transverse part 



-4cut(w,k) 



Re( 



arctan 



A3 

Agi 



M2, 



arctan 



M2 



(3.118) 



For all C ^ u < Agi it is Agj/M^cj* ^ 1 and the first arctan in the squared brackets may 

be set equal to 7r/2. Considering <C M the argument of the second arctan is X'^/M'^lu* <^ 
1 and the combination of both arctans is ~ 1. Increasing the energy to a; ~ M one finds 
X^/M'^lo* ~ M/{M — S^), which becomes large only for ^ M. However, since the integration 
over ^ stops here anyway this case does not have to be analyzed further. Consequently, also for 
u; ~ M one may estimate the arctans to be ~ 1. For uj M one has }? /M'^uj* ^ 1 and the 
arctans finally cancel. Considering C, < M one finds for u; <C M that }? /M'^uj* ~ M/uj* ^ 1 



and the arctans cancel. The same is true for uj ^ M since then A^/M^ 



{u/MY > 1. Only 



in the region u; ~ M it is A^/M^u;* ~ M/uj* ~ 1 and the arctans do not cancel. 

One first considers C, M and to ^ (p. Because in this special case the integral over does 
not yield the BCS-log one finds 



,k) 



(3.119) 



For larger energies a; ~ Ay with < y < 1 and <C M one substitutes S,{y') ^ Ay', d^/^ 
ln{(t)/M) dy', one finds with Eq. (I3.11()|l 



M 



dy sin 



Try 



g (p cos 



Try 



(3.120) 
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Following the discussion after Eq. (|3.118() the limit y — > of Eq. H3.12U() is valid also for C ~ M. 
Hence for u M and all < Aq it is ^*ut(a;,k) ^ g cp . For all other values of to and C the 
contribution Al^^{uJ,'k) is strongly suppressed. 

In the longitudinal sector one finds for the integral over the gluon momentum p 



X(A) 



dp 



1 



(p2 + X2)2 X 



arctan 



Agi 
X 



, A 

arctan | — 



A 



X2 + A2 



1/X, forA<X 
1/A, forA>X 



(3.121) 



where it is abbreviated X'^ = 3m2. Since C < A^^ ^ X, solely the magnitude of to decides 
whether A < X or A X is realized. It follows that in contrast to the transversal case the 
order of magnitude of A^^^ is not dependent on the choice of C- Energies a; ~ with < y < 1 
correspond to A < X, where the r.h.s. of Eq. (|3.121|) is ~ 1/X. Beginning with the special case 
a; ~ Ai one finds 



AUu;,k) ~ 5' / -Re 

J £a 



M 



(3.122) 



For uj ^ Ay with < y < 1 one has similarly 



y 

-4cut(^^'k) ~ g<P j dy' 



sm 



Try 



g (p cos 



7ry\ 

2 y \M 



M 
y 



gcj) dy' sin 



Try 



m) 



m) 



(3.123) 



Hence, in the considered energy regime ^^ut is suppressed by the factor {(p/M)^ as compared to 
i.e. <C ^cuf For energies a; ~ M the longitudinal and the transversal cut contribute 
at the same order, ^^^^ ~ ^ gcj). For much larger energies, M ^ w < (C is bounded 
by Aq) it is A = ~ a; » X and the r.h.s. of Eq. (|3.121|) is ~ 1/A ~ l/uj. It follows with 
(5 = Aq ~ Ao 



* 

^Lt(u;,k) ~ g cP J dy sin (^^) ^ ^ g cP J dy sin (^^) ~ 5^ (3-124) 
1 1 

and one found that ^^ut ^ -^cut this large-energy regime. The results for ^^^^ are summarized 
in Tables EU and O 



Undamped gluons contributing to A 

The contributions from the undamped gluon excitations to A read analogously to Eq. H3.116|l 



<5 2/i 

4'*i,(u;,k) ~ J^Re(^[e^,ci) J dpp //^.^{lu* ,p) 5[u;* - u;e,t{p)] ■ (3-125) 
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Ijj 


4> 


^l>j(>0 


M 








g(\) cos (2^) 







A'- 




5 '/'(m)" cos (2^) 


94> 


g4> 


•^cut 


5V 


ff</>cos(2^) 


g(t> 


94> 



Table 3.1: Estimates for AI'*^^ and Acut — -^cut + -^cut different energy scales and C, <C M. 



to 




Ai>j,>o 
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M < < /i 


A^ 













A'^ . 




5(/.(^)'cos(2^) 




94> 






50 U)' cos (2^) 


94> 


94> 



Table 3.2: Estimates for AI'*^^ and Acut = A^^t + -^cut different energy scales and C ^ M. 



Due to the 5— function and uj£^t{p) ^ n^-g this contribution will be nonzero only for energies 
IV > nig ~ M. The upper boundary ~ 2/i in the integral over p is due to the constraint on 
relevant quark momenta near the Fermi surface in the effective theory, ^ < Aq, where Aq ~ gfj, 
is the quark cutoff, cf. Fig. 13.91 Due to this kinematic limitation it immediately follows that 
•^poie ~ fo^' cij > 2/i + Aq ~ 2/i, since then uj* = u — > u>£^t always and the (5— function on the 
r.h.s. of Eq. (|,3.125j) is identically zero. 

For the transverse sector one may approximate for all momenta p 

/5^„i,(a;t(p),p) =.-^-1^ , (3.126) 



and LOt{p) — yp^ + m^. One finds after substituting dpp ~ dtot and assuming C ^ M 

Aq 2/Lt 

^max 

~ 5' / — Re0(eq,q)~5</', (3.127) 



where ^max = min[(a;^ — m'^)/{2uj), Aq]. In Fig. 13.101 the integration regions for ^ and ut are 
visualized. For ^ > ^max it is = u; — ^ < ^ + < ujt and the (5— function under the 
integral over ujt is always zero. The last estimate on the r.h.s. of Eq. ()3. 127(1 is valid for energies, 
which are at least uj > nig + Ay with y < 1, because then ^max is at least ky and the BCS- 
log is generated cancelling one power of g. For C ^ -^^ one has p > |C ~ CI and therefore 



uJt > y "1,2 + 1^ - CP. For Ai < ^ < Ag it is uJt - ^Jmj + C > V^nig, while for Ag < ^ < Aq it 
is UJt > nig (independent of the signs of ^ and Q. It turns out that for uj exponentially close to 
nig, i.e. for energies of the form uj ~ nig + Aj^ with y > 0, w* = — ^ = cannot be fulfilled 
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in the whole integration region Ai < ^ < Aq. Consequently, A^^i^ = in this energy domain. 
Only for energies uj > ^ + > \f2mg one finds again J^^\^ ^ g4>i which remains valid up 
to a; < 2jLt. 

In the longitudinal gluon sector one may approximate for gluon momenta p < nig 

pUUMpU) ^ (3.128) 



and iOi{p) ~ + "^g- = '"^3 + with < y < 1 one finds after substituting dpp 
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Figure 3.10: The integration regions of ^ and cut in Eq. 1)3.127(1 . Only those values of ^ contribute, for 

which the integral along ujt, starting at LUt = yrn^+^ and ending at ujt ~ 2/i, has a point of intersection 

with Lo* ^ Lo — ^, LUt = w*. Only for lo > rrig + Ay with y < 1 the integration region is sufficiently large 
to generate the BCS-log. The dotted lines correspond to different values of ut. 



duJi uj£ and assuming C <C M 

Aq \/2mg 



-4poie(^. k) ~ g"^ f — Re </>(eq, q) f duji ^ 5[u)* - uJi] 

J 6a J L0£ Tflg 





in 



~ I ^Re.^(6q,q) . (3.129) 

where the upper limit ^max = {^'^ — ''^g)/(2w) ~ Aj^ is obtained analogously to the transversal 
case. Since .^max ^ and uj > mg one may approximate the energy fraction under the integral 
by jiuP' — rrig) ^ u)/{lo — nig). Then the integral over ^ may be readily performed yielding 

Al,,,iu;,k)^gcj>^ gc^ (^Y , (3.130) 

^ u; -nig \(p J 

where the BCS-log has cancelled one power of g. Analogously to the transversal case, for 
C M one finds A^^^^^ = if w is exponentially close to rUg. For lo > Jrnj^ + ^ \pimg it 
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UJ 


<mg + Ki 




> V2mg 


nig <^ LU < 2fi 


A* , 

•^polc 
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A^ , 

•^polc 
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•^pole 
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Table 3.3: Estimates for y^p*j(, and y^poic = -^poic + "^poic different energy scales and C, <C M. 



UJ 


<mg + Ai 


^ mg + Ai>y>o 


> V2mg 


nig <^ LU < 2fi 


A^ , 

•^pole 








94> 
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A'' , 

•^polc 








94> 




"^polc 
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94> 



Table 3.4: Estimates for A'^^^^ and y^poio ~ -^poic + -^poic different energy scales and ( < M. 



follows ^poig ^ 9 4^ again. For much larger energies, uj ^> mg, gluon momenta p ^> nig have 
to be considered in the estimate of ^poie) cf. upper dotted lines in Fig. n.lOl For those the 
longitudinal gluon spectral density becomes exponentially suppressed 



Ppole(^^«(p).P) 



exp 



P 



(3.131) 



and one has for nig <^ uj < 2fi with wz(p) — p 

Aa 



Aq 2ft , 

-4pole('^>k) - 5(^y" — Re(/)(eq,q) J duJe exp i 



9"(P I -J exp 

Ai 



\/2n 



3m| 



3m2^ 



9 9 exp 



3m2 



(3.132) 



which is the continuation of the estimate given in Eq. H3.13U|) to large energies lo ^ rUg and is 
valid for all C ^ Aq. The results for .Ap*[g are summarized in Tables ESI and EiH 

Landau-damped gluons contributing to B 

In the Landau-damped gluon sector A4g*j.^Q in Eq. H3.1U9() contributes to B as 

Aq UJ Agl 

^fut(^>k) g'^ — dqoJ2 P</.(go,q) / dpp/9cut(^'.P) > (3.133) 

where = a; — (?o < and A = max(|^ — Cli"^')- Substituting the estimates of A{qo,(i) for 
p^{qo, q) one estimates B^'^^{ui, k) for different domains of oj and C- Due to the condition A < Agi 
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in Eq. and A{qo,q} = for qq > 2/x it immediately follows that B^^^^{u!,'k) = for 

uj > Agi + 2^ ~ 3fi. Inserting the approximative forms 117(1 for /)^^^ into Eq. 133(1 the 
integration over p can be performed analogously to Eqs. ((3.11813.121(1 . In the transverse case 
one finds 



aA{qo,q) 



crer 



arctan 



M2 



UJ' 



arctan 



A3 



M2 



UJ' 



(3.134) 



The analysis of the domains of w, C and ^ where the arctans in the squared brackets in Eq. 
((3.134(1 do not cancel is analogous the one after Eq. ((3.118(1 . For all ^, C < Aq and u; < 3/i it is 
Agj/M^cij' S> 1 and the first arctan in the squared brackets may be set equal to it/2. Considering 



(^,uj <^ M the argument of the second arctan is /Al'^uj' 
this is not ^ 1 and the combination of both arctans is 

finds that \^ /M'^uj' becomes large only for go — > a;, because then /M^uj' ~ 



~ i^jM'^J . As long as C < (Af2a;')i/3^ 
1. Increasing the energy to w ~ M one 

v~eVM'(w-go) » 1- 

However, since the integration over go stops here anyway this case does not have to be analyzed 
further. Consequently, also for a; ~ M one may estimate the arctans to be ~ 1. For uj ^ M one 
has A3/M2 uj' ~ [(w— go)/M]2 which is < 1 only for u) — M < go ^ w. For values of w outside this 
range, the arctans cancel. Considering C ^ one finds for uj <^M that X? /M'^uj' ~ M/uj' ^ 1 
and the arctans cancel. (In the region where l^^ — CI is sufficienty small the arctans do not cancel. 
However, since this requires < M, the BCS-log cannot be generated. Therefore, this special 
case can always be neglected.) In the region u; ~ M it is X? /M'^uj' ~ M/uj' ~ 1 and the arctans 
do not cancel. For uj ^ M the respective analysis as made for C, <^ M applies. 

Before proceeding further one proves that the generation of the BCS-log can be prevented if 
additional logarithmic dependences appear under the integral over ^ in the following form 



Ao 



^ln^i±^U/$ln^^ + ^^ 



(3.135) 



where < y < 1. Introducing the dilogarithm |129j 



u 

U,{x)^ l^ln(l-e) 



(3.136) 



which has the values — ^tt^ = Li2 (—1) < Li2 (x) < Li2 (1) = gvr^ for —1 < x < 1, one may 
express the first term on the r.h.s. of Eq. ((3.135(1 as 




Li2 (-1) +Li2 




Lio 




(3.137) 
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Since < Ai/Ay < 1 this term is of order 1 and no BCS-log has been generated in this term. 
The second term on the r.h.s. of Eq. (||-i.l,''{5|) is 



Ao 



Ao/Aj, 



Ay/Ao 



^In 

X 



IT 



+ Li2 



1 + e 

1-x 



Ay 

Ao 



Ay/Ao 



^In 

X 



1 + 



1 

Li.(l)-LM-l)+Li.(-^)-Li,(^) 



Lio 



Ao 



(3.138) 



In the second step one substituted x = 1/C with dx/x = Similarly to Eq. 1)3.1371) it is 

< Ay/Ao < 1. Hence, also this term is of order 1 and no BCS-log has been generated here, 
either, proving the above statement. 



Beginning with energies uj 
Assuming <C M one has 



one may use Eq. (|3.119|) to estimate A ^ A', 



cut 



9 



rIC 

J J L 



a 



90 - f^eq 



A3 



arctan 



M2 



UJ' 



Vl/3 



In 



arctan 



A3 



M2 



(3.139) 



The arctans in the squared brackets cancel for 3> M'^cp, which leads to the upper boundary 
Ai/3 in the inegral over ^. As discussed before the generation of the BCS-log was prevented by 
the additional logarithm under the integral. For C ~ the arctans in Eq. ()3. 139)1 would have 
cancelled yielding B^^^{cj),'k) ~ 0. 

For Lo ^ Ay with < y < 1 one conservatively estimates A ~ w4*^t ~ 5 cf . Eq. 1)3.120)) , 
and obtain similarly to Eq. ()3.139)) . assuming <C M 



^*,,(A„k) ~ g'ct^l^ I dqoY. 



a 





en 



arctan 



arctan 



A3 



eq- Aj, 



en + A^ 



(3.140) 



where again no BCS-log was generated. Here, the arctans cancel for A ^ M'^Ay, leading to the 
constraint < Ay/^. Again, for C ~ the arctans in Eq. ()3.139)) would have cancelled yielding 
^*„t(A„k)~0. ■ 

For energies uj ~ nig + Ay with < y < 1 the combination of the arctangents in Eq. ()3. 134)1 
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is ~ 1 for all C ^ ^q- Integrating over qq from to nig one finds similarly to Eq. H3.14U() 



9 



In 



eq - Ao 



eq + Ao 



9 



(3.141) 



The dominant contribution, however, is generated when integrating over qq from nig to nig + Ay 
where one has A ~ ~ 5 </> {M/4>)y, cf. Eq. (jHSII), 



Bi^^inig + Ay,k) ~ 



a 



mg+Ai 



y Aa 



(Tmo 



(3.142) 



For energies iv ^ nig but not exponentially close to niq, the integration over qq from 2mg to a; 
may be performed after estimating A ^ gcp, cf. Eqs. 1)3. 12413. 129|1 . yielding 



dqo 



a 



2mo 



5^0 



In 



^ - eg 



In 



uj- M + eg 
uj - M -en 



~ . (3.143) 

Hence, one found that generally for energies lo > nig it is 

^*,,(a;,k) ~ 5' 'A- (3.144) 

In the limit of very large energies, uj ^ M, it is }? / {M'^ui') = [uj — q^)^ jM'^ < 1 only for 
a; — M < go ^ i^- For values of go outside this range the arctans in Eq. (|3.134j) cancel. Inside 
this range one may estimate A ^ gcp and find to be strongly suppressed 



BUu;,k) ^g'4> f- f dqo ^ ^ 9'<P I di ( 



Ao a; 



dqo 3 



MV 

UJ J 



(3.145) 



LU-M 

In the case of longitudinal gluons one obtains 

Aq iu 



Ai oj~M 



^Lt(^,k) ~ g'f^fdqoY: 

J en v/ ^ I 
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(3.146) 



where X(A) is defined in Eq. ()3.121|) . Analogously to the analysis of A^^^ one finds for lu 



bU<P,^) ~ 9' 



Aq 4> 



dqo Y 
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<^ciJ ^±qo- c^eq M 
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(3.147) 
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Table 3.5: Estimates for ySf^t and Scut = ■Scut + ^cut different energy scales and ( ^ M. 
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Table 3.6: Estimates for yS^ut and Scut = Scut + '^cut different energy scales and ( < M. 



and similarly for to ^ Ay with < y < 1 



(7 A, 



Cq J rrl qo - 0-eq M 

'^-='= 



M 







eq - Aj^ 


eq + Aj, 


i Cq 
^ 



(3.148) 

For UJ ~ mg + Aj^ with < y < 1 it is to' Z{X) ~ 1 and one finds as in Eq. (|3.142j) 

Bl,{mg+Ay,k) ~ g^ct). (3.149) 

Also for energies u; > M the analysis is very similar to the transversal case since uj'Z{X) ~ 1, 
and one finds again 



BU^,k)^g^ 



(3.150) 



In the limit of very large energies, uj ^> Agi ~ /u, only the range to — Agi < qq < lo contributes, cf. 
Eq. (|3.133j) . As soon as w — Agi > 2mg one may estimate A ^ g (j) and obtain with lo' Z{X) ~ 1 



Jlq U) 

BU-^,k) g^cpf^ f dqoY: 



a;-Agi 
MAgl 3 , M 



± go - creq 



LU-Arr 



g <p — • 

Ijj 



(3.151) 



Hence, also becomes small in the limit of large energies, cf. Eq. H3.145|) . 
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Undamped gluons contributing to B 

In the undamped gluon sector the term TW^'^^q in Eq. (|3.1U9|) gives the contribution 

Aq uj 2/J 

Koiei^, ^)-9' f-fd%T. ^M^^q) dppp^i^^iu', p) 5W - ^,,(p)] , (3.152) 

which is nonzero only for energies w > rrig. Due to the restriction p <2pm. Eq. 1)3.152(1 it follows 
with similar arguments as for iSf^f^ that i3p'*[g(u;,k) = for uj > 4/i. For transversal gluons one 
finds using the same approximations as for A^^^^ beginning with energies u! ~ nig + Ai<y<o and 
assuming (" <C M 

A,, 

d^ f ^ cr^(90,q. 



Ao oj 2fi 



<^-\/ni|+F Ay/2 



A 



ff'</' , (3.153) 



where to guarantee uj — y + > the upper boundary of the integral over ^ is reduced 

from Aq ~ Aq to the scale ^JrUgKy ~ Furthermore, one estimated A ^ gcj), which is valid 

since (/o < A^ < nig for the considered energies. For uj > 2mg the integral over go also runs over 
values go > iT^g receiving contributions from ^p^ie, cf. Eq. (|3.13Up . As a consequence one finds 

^*,i,(u;,k)~g2^ (3 154) 

analogously to Eq. ()3.142|) . For uj » nig the additional contributions from 2mg < qq < to are 
only ~ as can be seen in the same way as in Eq. p. 143(1 and 13^^^^ ~ g'^cj). However, for 
UJ > 2fj, + 2mg the condition uj' = uji can be fulfilled only for > uj — 2^1 > 2mg, where A ^ §</), 
and one may estimate 

J i J ±1 go - O-Cq J J Qn 



'^-^ " ' Ai u)-2^l 



M2/i 3,M 



~ /0^^~5'</'— . (3.155) 

O)^ Ul 



For C ~ and uj ~ nig + Ai<y<o the condition u — \Jnfig + |^ — CP > leads to |C| — < 

C < Id + Ay/2. Then the integral over ^ finally yields B\^^^^ ~ g^4>{Ay/2/M) ~ g^cj) {(j)/M)y/'^ . 
For > 2mg and uj ^ nig the same analyses as in the case oi <^ M apply. 
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In the longitudinal sector starting with energies co = nig + Aj, with < y < 1 and momenta 
p ^ mg one may employ Eq. (IIH28I) . Assuming <^ M one obtains 



aA{qo,q) 



Ai 



dcoe 



Lot — TWt 



dqo 



aA{qo,q) {u - qof 



(3.156) 



where upper boundaries of the integrals over ^ and qo are analogous to the transversal case, cf. 
Eq. p.l53p . Since qo < to — ^J^"^ + f^l is much smaller than uj ~ rrig + A^^, one may neglect qo 
against u; on the r.h.s of Eq. (|3.156|) . Therefore, one may estimate A ^ g (j) and finally obtain 

■"pole 



9 



up' — rv?g 



ijj 



UJ — m„ 



9 



di 


In 


UJ — 




UJ — 




M 




.( 




)'■ 



(3.157) 



For larger energies uj > 2mg the upper boundary of the integral over qo will just exceed rUg, 
where it is ^ ~ ^p^ig, cf. Eq. 1)3.130^ . One finds analogously to B^^^^e ^^^^ this gives the main 
contribution 



^pole(^>k) 



Ao u) 



cT ^(go,q) 



Ai "ifl 



dujf. 



UJo 



dqo 



a 



ujj - mj 



5[uj' - UJl] 



qo 



{uJ - qof 



^ qo - o-eq 90 - mg {uj - qof - m 



(3.158) 



where one exploited {uj — qo)'^/[{u} — qoY ~ ^-g] ~ 1 (since always uj — qo ^ mg) and estimated 
the integral over go as in Eq. (|3.142|) . For energies uj ^ 2mg one finds analogously to Eq. 1)3. 132^ 
for the contribution from the integration region 2m g < qo < uJ, where one may approximate 
Ar^g^cj), 



9^'P — / dqo V / duji exp -t^— | 

J J qo - o-eq J \ 3m^ 



5[uj' - UJl] 



2m 
Ao 



2mg 



qo - 



exp 



2fa; 



3m^ 
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Ijj 


< rn-g + Ai 




> 2mg 


rUg <^ UJ < Afi 


*-'pole 











^^ 

*-^polc 











^pole 







9^4> 


9'^§ 



Table 3.7: Estimates for /Bp'^j^, and Spoic at different energy scales and ^ ^ M. 



LO 


< mg + Ai 


^ mg + Ai>y>o 


Z'^rng 


rrig <^ LU < Afi 


^pole 





9^<^ {0'' 






pole 









9'<P i^Y 


^pole 












Table 3.8: Estimates for B^^^^ and Spoio ~ B^^^^ + B^^^^ at different energy scales and C ^ M. 



2mg ^ ^' 

Hence, in the considered large energy regime, u; 2mg the main contribution to B^^^^^^ comes 
from Eq. H3.158() . and it is B^p^Xc ~ 9^^- Similarly to the transversal case, this holds up to 
energies u; > 2/x + 2mg since then it is always go > Irrig leading to 

Bj„i,(u;,k)~/</< (^-j . (3.160) 

Analogously to the transversal case one finds that for C ^ ^ ^ind uj ~ ?7Zg + Ai<j,<o the condition 
ijj — ^Jfri^g + \i — CP > leads to \C,\ — < ^ < ICI + \/2- Then the integral over ^ finally 
yields ^^^jg ~ g^(j){M/(j))y {Ay/2/M) ~ g^cj) {M/(p)y/'^ . For to > 2mg and w > the same 
analyses as in the case of ^ <ti M apply. The results of this subsection are summarized in Tables 
OandEHl 



3.2.3 Estimating n[A] and n[B] and 0o 
Hilbert transforming A and B 

Having determined the order of magnitude of Im (j) at various characteristic energy scales one is 
now in the position to (qualitatively) Hilbert-transform Imcp for ^ <C M and to find the order 
of magnitude of Rec^ and the corrections due to i3. To this end one splits the integral over to in 
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the dispersion relation Eq. H3.74a|) as 



Re,^(ek,k) = V [ dc^M^ = p [ du Y 



/O0(w,k) 



uj — uek 



V 



Ai Ag Aq 2mg 2fi A/l 



+ + 



+ + 



Ai Ag r?i„+Ai 2mg 2^ 



u) — fiek 



(3.161) 



Since the integral is additive, one may choose for each region the leading contribution. At the 
first scale < a; < Ai it is ~ -^iut ~ S'^'/'i cf. Eq. (|3.119|) . and one has the contribution 



Ai 



uj — crek 



(3.162) 



At the scale Ai < w < Ag one has ~ A\.^^_ ^ d^P^ cf. Eq. (|3.12Up . With the substitution 
du/uj = ln(0/M) dy one finds 



vjdu^Y. 

A. 0-=± 



UJ — crek 



Ao 



9(t> J ^~5'Aln(^) Jdy^cp. (3.163) 



The contribution from A^^^ ~ gip^^'/M)^ at the same scale, cf. Eq. 1)3. 123^ . can be shown to 
be much smaller. One has 



r fduY: 

I <T=± 



UJ — crek 



M 



dy 



-V 



M 



M 



(3.164) 



At the scale Ag < u; < Aq one has ~ A\.^^ ~ w4^ut ^ g<l), cf . Eq. (|3.12Ul 13.123(1 , and finds 



Ao 

vJdujY. 



Aut('^,k) 



Ao 



duj 



UJ 



(j) dy gcj) . 



(3.165) 



For energies uj > nig + Ay with < y < 1 one has ~ -^poie 9 't' {M/cj))^, cf. Eq. (|3.13U|) . and 



finds with duj = ln{(f)/M) Ay dy 

rrig + Ao 



V 



j dujY. 



mg+Ai 



Integrating over 2mp < uj < 2p with p,^ ~ ^p^ig g (p, cf. Eq. (j3.127j) . one obtains 

V f duj — — ^ — ^ ^ qd) f — ~ (J (/> In ^ '-^ q d> . 
J ±1 w-ack J UJ \MJ 

2mg 2mg 



(3.167) 
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Finally, integrating over 2^ < lo < 4fi with ~ 'Bpoie ~ d^^PiM/uj), cf. Eqs. H3.151I3.155|) . one 
obtains 

p/do; V^^^^i^^~/0M/^~/0^~5V. (3.168) 

As expected, Rec/) ~ (j). Furthermore, I find several potential sources for sub-subleading order 
correction from to 7{[B] to Kecf). The first comes from B^^^, which is ~ 5^^cut fo^ w ~ 
with < y < 1, cf. Eq. (|3.14Up . After Hilbert transformation it yields a contribution of order 
g'^cj) to Keep, cf. Eq. (|3.163j) . and is therefore of sub-subleading order. At this point, also from 
•^poie ~ ff^^poie ^ sub-subleading order contribution seems possible, since the corresponding 
contribution from A^^i^ is ~ (/), cf. 1)3. 166(1 . As the latter, however, combines with (pQ to a 
subleading order term, cf. Sec. I3.2.4( it would be interesting to investigate if also B^^i^ finds 
an analogous partner to cancel similarly. If not, it would contribute at sub-subleading order to 
Re4>. Furthermore, one found that S ~ c/^ ~ g'^cj) for uj> M, cf. Eqs. (13.1 4413.1 5013.1 5413.1 5SI) . 
From the estimate in Eq. (|3.167|) one deduces that the corresponding contribution of TC[B] to 
Re0 is of sub-subleading order. In the next section it will be analyzed at which order Imcj) 
contributes to the local part of the gab function, (j)Q. 



The contribution of Im (p to (pQ 

The gap equation for </>o(k) is obtained by considering the integrals Iq and Ik^ , cf. Eqs. (|3.1U1I3.1U4() . 
in the limit kQ — > oo. Since p < 2|U the gluon spectral densities p^'^{qo,p) are nonzero only for 
Qo ^ 2^. Consequently, the integral over qq in Eq. H3.1U4() is effectively bounded by ±2/u. Then, 
due to the energy denominator under the integral, /kg tends to zero as /cq ^ oo. In the sec- 
ond term on the r.h.s. of Eq. (|3.1U1|) one has eq < Aq ~ gfi. Therefore, one may neglect the 
imaginary part of the gluon propagator A^'*(ibeq — kQ,p), cf. Eqs. (j3.33j) . for ko ^ 2iJ, > p . In 
the transverse case the gluon propagator becomes A* ~ l/A:o and the respective contribution 
vanishes in the limit ko oo. In the longitudinal sector one has —l/p^. Hence, the 

longitudinal contribution of the considered term does not vanish. In the first term on the r.h.s. 
of Eq. (|3.101|) the integral over only runs to values qo ^ n due to the presence of p^. In the 
limit fco — > oo again only the contribution from the static electric gluon propagator A^ — > —l/p"^ 
remains. Consequently, one finds for </>o(k) to subleading order 



00 (k) 



g f f dp 



3(27r)2 Up 



Tr^{k,p,q)< 0o(q) + Re (/)(eq + zr/, q) Z (eq) tanh 



2T 



-p I du^^P^^Z\Lo)tanh 



2T 



(3.169) 
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In the limit T ^ the hyperbohc functions simphfy, yielding after performing the integral over 
p and with Til{k,p,q) ~ 1 

Ao 

<^o(k) ~ 5'/yln(^){[0o(q)+Re<^(eq + ir?,q)' 

Ai 

7 1 

-V dojYl P4>'y^^^) \ ■ (3-170) 



Using Re (/'(sq + irj, q) ~ (/> one finds for the respective contribution in Eq. (|3.17Uj) 

Ao 

^V/f lnf^U</>, (3.171) 



Ai 



and hence (;/>o(k) ~ (p. Consequently, the integration over the first term in the squared brackets, 
(/)o(q), also gives a contribution of order (p to (/>o(k). The remaining term is the contribution 
from p0 and is identical to Eq. (|3.161|) up to an extra sign a arising from the hyperbolic tangent. 
One can conservatively estimate this term by appoximating ^ gtj) for < a; < 4/i and all 
Ai < ^ < Ao and adding Pcj, ^ g (j) {M/(p)y in the range u) ~ rUg+Ky , 1 > y > and Ai < ^ < Ag. 
One obtains 



A 



Ao 



g'(P -^j di In ~ ^ ~ (3.172) 



A 



and 



Kg mg+Ko Ag Q 



mg+Ai Ai 

As 



^Ai '^=± 
Ai 

where one has estimated ^rj(y/{oj — a^) ~ i/M, since ^ <C w ~ M. It follows that the 
contributions from p^ to 0o are of order g'^cf) and hence of sub-subleading order, cf. discussion 
after Eq. (131681). 

Finally, one has found that the contribution of Im (f) to Re k) = Re </>(ek, k) + (/'o(k) are 
in total beyond subleading order. 
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3.2.4 Reproducing Re (/)(ek + ^r^, k) to subleading order 

Although /9(^(a;,k) and </>o(k) have not been calculated (in principle possible) but only estimated 
in magnitude, the main goal of this analysis has been achieved. The imaginary part of the gap 
function is shown to be at most of order gRe(/)(ek,k) for ek < Aq ~ gp (then it is = Re0, 
cf. Eq. (|3.78j) ) and it enters Re(/)(ek,k) beyond subleading order. 

Neglecting all terms that have been identified to originate from p^, namely ImA4i3(u; + 
q, p) in Eq. (|3.106|) . and Hilbert transfoming the remaining term Im7W_4(u; + ir;, q, p) gives 
the gap equation for Rei?i'(ek,k) to subleading order 



Re 0(ek + irj, k) 



Adding 



00 (k) 



3 





r d^q z 


3 „ 


1 (27r)3 






— c 


» -=± 


1 r 

X - 

2 . 


tanh I — - 
V 2T 





■ Re (/)(eq + ii], q) 



a [Triik,p,q)p'{u;,p)+Trlik,p,q)p\u;,p] 



w - ek + aic 



-) 

2TJ 



Re (/)(eq + iri, q)^ tanh 



(3.174) 



(3.175) 



to Eq. (|3.174j) one reproduces the well known gap equation for Re 0(ek, k) to subleading order, cf. 

To see explicitly how 0o ~ </> and '^[^pQig] ~ 4> combine to a subleading order contribution, 
one first notes that all terms ~ coth are at most of sub-subleading order, cf. Sec. 13.2.11 Then 
Eq. (|3.174l) simplifies to 



Re(^(ek + i??,k) 



f d^q, Z2(eq) 



3 J (27r)3 2ec 



Re 0(eq + irj, q) - tanh ( ^ 



xV / dojJ2 



Tiiik,p, q)p'{u;, p) + Tii{k,p, q)p\u;, p) 
w - ek + O-Cq 



3(27r)2 

E 



Z2(e ) 

d^ — ^, ^ Re 0(eq + ir], q) tanh 



2er 



dpplTrl{k,p,q) 



2T 



+ ^HDL(ek-o-eq,p) 



2fi 



+ J dppTi:l{k,p,q) 22(ek - o-fq, p) 



+ Tr*(/c,p,g)A*HDL(ek - o-eq,p)| 



(3.176) 



where one used p^{LU,p) = for p > Agi in the effective theory, cf. Eq. ()3.86a() . Adding Eqs. 
1)3. 17513.176(1 . the 1/p^-term from the soft electric gluon propagator in Eq. I|3. 176(1 effectively 
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restricts the p— integral in 0o from Agi to After aproximating the hard magnetic gluon 
propagator as Ao22(ek — o"eq,p) = l/p^ + 0(Aq/Agi), one can combine it with the remaining 
contribution from (fiQ. Using Tr^(/c,p, q) — Tr* (A:,p, = 4 + 0(Aq/Agi), one finahy arrives at Eq. 
(124) of |98i. Hence, indeed (po + '^[^pQjj,] ^ g(j), i.e. contribute at subleading order to Ret/). 



3.2.5 Calculating Im 0(ek + «?7, k) near the Fermi surface 

In Sec. 13.2.21 the contributions A and B to Im (p have been estimated for different regimes of 
u and C, cf. Tab. 13.1113.81 In the case that both u and C are ~ Aj^ with some y satisfying 
< y < 1 one found ^*ut ^° the dominant contribution to Imc/), while in all other regimes 
of oj and C, different gluon sectors contribute at the respective leading order and therefore mix. 
Furthermore, for uj < 2mg also the contributions from B have to be considered, since there B is 
suppressed relative to A only by one power of g. 

In the following the imaginary part of the on-shell gap function Imi;^>(ek + ??7,k) near the 
Fermi surface, ^ ^ Ay with 1 < y < 0, at zero temperature, T = 0, is calculated. It follows 

lm<p{ek + i7],k) ~ ^%^)2 J y ^^(^q) (/)(eq + ir/, q) J dpp pl^^{uj* ,p)Til{k,p,q) 

Ai A 



— Z'^{eq)Re(p{eq + irj,q) J dp 



3(2^)2 7 ^ v^qy^.-v^v^q , y ^ p6 + (Af2^*)2 



A 



1 



9(2^)2 J ^ 



J ^z2(eq)Re</.(eq + ir/,q) 



y 

In [ — ] (p dy sm — — 



9(27r)2 \M J ^ J " V 2 



- ^ „ In 




187r2 








^ cos 









0COS (^^j +0(5^ 

+ 0{f) , (3.177) 

where one substituted ek = A^, and d^/$, = dy' ln{(j}/M) and used ln((/)/M) = — 37r2/(v^gi) 
and g = g/{3\/2Tr). Furthermore, it was sufficient to approximate Tr*(fc,p, g) ~ —2. The 
corrections due to the term —p'^/2fi'^ are easily shown to be suppressed by a factor ~ {M/ fi)'^ ~ 
g2 compared to the final result in Eq. ()3.177() . To estimate the contributions from the term 
— {k"^ — (f')'^ / {2 k qp^) ~ — 2(^ — C)2/p2 one first conservatively sets C = 0- At the lower boundary 
of the p— integral, p > A, this term is ~ 1. One finds for the integral over p 



M2 lo* ^2 / A; 



A 



15 
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The integral over ^ then can be estimated to be at most of order 

4 iy 

y , iy' y 4y' / 4> \ 15 / 6 \ 15 Ay 

h' - m i^r <h' i^r ^^^^ h^t- '-™> 

Hence, also this contribution is suppressed by a factor of g compared to the result given in Eq. 

mm- 



Chapter 4 

Summary and Outlook 



In this work I have presented a formal derivation of a general effective action for non-Abelian 
gauge theories, Eq. (|2.4(i)) . This was motivated by the occurence of well-separated momentum 
scales in hot and/or dense quark matter. To this end I first introduced cut-offs in momentum 
space for quarks, Aq, and gluons, Agi. These cut-offs separate relevant from irrelevant quark 
modes and soft from hard gluon modes. I then explicitly integrated out irrelevant quark and 
hard gluon modes. The effective action (|2.46j) is completely general and, as shown explicitly in 
Sec. I2.2.T1 after appropriately choosing Aq and Agi, it comprises well-known effective actions as 
special cases, for instance, the "Hard Thermal Loop" (HTL) and "Hard Dense Loop" (HDL) 
effective actions. I also demonstrated, cf. Sec. 12.2.21 that the high-density effective theory 
introduced by Hong and others |1UL)1 IIUIL I1U21 11U31 11041 llL)5j is contained in the effective action 

In Sec. 12.31 it is argued that in order to obtain a complete effective theory, in principle, all 
diagrams appearing in the effective action ()2.46|) have to be power-counted. This remains as a 
future project. A first step towards this goal, however, is performed for the special case of an 
effective theory for cold and dense quark matter, i.e. for quarks with momenta |A; — /u| ~ around 
the Fermi surface and soft gluons with momenta of order A^j <C fi. To this end, one specific class 
of diagrams occurring in the effective action (|2.4(i)) . loops of irrelevant quarks with external 
soft gluon legs, is power-counted for the cutoff parameters fulfilling </) <C Aq ^ Agi ~ /i and 
the projection operators given in Eqs. (|2.18|) and (|2.35j) . It is shown that for the considered 
momentum regime these loops are of the order of the corresponding bare vertices (bare gluon 
propagator and gluon vertices) times the factor {gfi)"^ /{A^^)'^ . They are therefore "relevant" 
operators in the effective action, as their magnitude increases when the soft gluon scale Agj 
decreases. They contain Debye screening and Landau damping (as known from HDLs), whereas 
the Meissner effect for the magnetic gluons is negligible as long as Agj ^> 0. A possible field 
of application could be the precursory effects to color superconductivity at temperatures above 
the (non-perturbative) onset of Cooper pairing, T > Tc ^ (p. 

As it is known from the HTL /HDL power-counting scheme there is another important class of 
diagrams, i.e. those with 2 external quark and N — 2 external gluon legs. It would be interesting 
to analyze how their orders of magnitude change as the cutoffs are shifted. Especially the 1-loop 
corrections to the quark-gluon vertex could contribute to the gap parameter at sub-subleading 
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order. 

In Sec. 13.11 1 showed how the QCD gap equation can be derived from the effective action 
(|2.46|) . The gap equation is a Dyson-Schwinger equation for the anomalous part of the quark 
self-energy. It has to be solved self-consistently, which is feasible only after truncating the set 
of all possible diagrams contributing to the Dyson-Schwinger equation. Such truncations can 
be derived in a systematic way within the general Cornwall-Jackiw-Tomboulis (CJT) formalism 
jlU6j . Here, I only include diagrams of the sunset-type, cf. Fig. 13.11 in the CJT effective ac- 
tion, which gives rise to one-loop diagrams (with self-consistently determined quark and gluon 
propagators) in the quark and gluon self-energies. 

In principle, the CJT effective action 1)3. 2|) for the effective theory contains the same in- 
formation as the corresponding one for QCD. However, while in full QCD self-consistency is 
maintained for all momentum modes via the solution of the stationarity condition ()3.3|) . in the 
effective theory self-consistency is only required for the relevant quark and soft gluon modes. 
These are the only dynamical degrees of freedom in the CJT effective action; the irrelevant 
fermion and hard gluon modes, which were integrated out, only appear in the vertices of the 
tree- level action ()2.46() . In this sense, the effective theory provides a simplification of any phys- 
ical problem that has to be solved self-consistently. How this facilitates the actual computation 
is exemplified explicitly with the solution of the color-superconducting gap equation. 

Usually, the advantage of an effective theory is that the degree of importance of various 
operators can be estimated (via power counting) at the level of the effective action, i.e., prior 
to the actual calculation of a physical quantity. This tremendously simplifies the computation 
of quantities which are accessible within a perturbative framework. On the other hand, the 
requirement of self-consistency for the solution of the Dyson-Schwinger equation invalidates any 
such power-counting scheme on the level of the effective action. For instance, perturbatively, 
the right-hand side of the gap equation ()1.35j) is proportional to g^. However, self-consistency 
generates additional large logarithms ~ ln(^/0) ~ I/5 which cancel powers of g. 

Nevertheless, it turns out that there is still a distinct advantage in using the CJT effective ac- 
tion (|3.2|) of the effective theory for the derivation and the solution of Dyson-Schwinger equations 
for quantities which have to be determined self-consistently, such as the color-superconducting 
gap function in QCD. This advantage originates from the introduction of the cut-offs which 
separate various regions in momentum space. Choosing Aq <C Agi they allow for a rigorous 
power counting of different contributions to the Dyson-Schwinger equation after expanding in 
Aq/Agi ^ 1. I explicitly demonstrated this in Sec. 13. H where I reviewed the calculation of the 
color-superconducting gap parameter to subleading order and in Sec. 13.21 where I considered 
the imaginary part of the gap function. 

In order to obtain the standard result 1)1.36(1 . it was mandatory to choose Aq < 5// <C Agi < 
cf. Fig. 14.11 This is in contrast to previous statements in the literature 10T\, 11021 ll()3j that a 
consistent power-counting scheme requires Aq ~ Agi. In particular, the choice Aq <C Agi has 
the consequence that the gluon energy in the QCD gap equation is restricted to values po < Aq, 
while the gluon momentum can be much larger, p < Agi. This naturally explains why it is 
permissible to use the low-energy limit ()3.61|) of the HDL gluon propagators in order to extract 
the dominant contribution to the gap equation (which arises from soft magnetic gluons). In 
previous calculations of the gap within the framework of an effective theory [lOH llfl2( 1103) , the 
low-energy limit for the HDL propagators was used without further justification, even though 
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for the choice Aq ~ Agi the gluon energy can be of the same order as the gluon momentum. 
The physical picture which arises from the choice < gfi <^ Agi < is summarized in Fig. 
14.11 Relevant quarks are located within a thin layer of width ~ Aq around the Fermi surface. 
Soft gluon exchange mediates between quarks within a "patch" of size ~ Agi inside this layer. 
The area of the patch is much larger than its thickness. Hard gluon exchange mediates between 
quark states inside and outside of the patch. 

Obviously, this picture, as well as all power-counting arguments, are rigorously valid only 
at asymptotically large values of the quark chemical potential, where ^ 1. In the physically 
relevant region, ^ < 1 GeV and 5 ~ 1, the scale hierarchy Aq < ^ Agi < /i breaks down. 
When all scales are of the same order, the patch on the Fermi surface becomes a sphere of the 
size of the Fermi sphere. 

In the course of the calculation, I was able to identify various potential contributions of 
sub-subleading order. However, I argued that, at this order, a self-consistent solution of the 
gap equation must take into account the off-shell behavior of the gap function. To this end, 
I investigated the imaginary part of the gap function, which (loosely speaking) generates the 
energy dependence of the gap function. Due to the energy dependent gluon propagator in the 
gap equation a self-consistent solution requires a energy dependent and therefore complex ansatz. 
Consequently, I considered a complex gap equation, which is pair of two coupled gap equations, 
one for the real and one for the imaginary part of the gap function. 
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Chapter 4. Summary and Outlook 



To solve it self-consistently, I first inserted the known leading order result of the real part of 
the gap function into the imaginary gap equation. In this way I estimated the order of magnitude 
of the imaginary part of the gap function for various energy and momentum regimes, cf. Tables 
I-III. The self-consistency of this approach was controlled by reinserting these estimates into the 
imaginary gap equation, cf. Tables IV- VI. Using these results it was checked at which order the 
imaginary part contributes to the real part of the gap function. It was found that it enters at 
sub-subleading order. Furthermore, several sub-subleading order terms in the real part of the 
gap equation arising from the imaginary part of the gap function have been identified. 

Finally, I calculated the imaginary part of the on-shell gap function for momenta exponen- 
tially close to the Fermi surface, cf. Eq. p.lZTf) . In this regime only Landau damped transversal 
gluons had to be considered. Also for the complex gap equation the two momentum cutoffs Aq 
and Agi proved to be a powerful means for the rigorous power counting of innumerable terms. 

For a complete sub-subleading order calculation it also appears to be necessary to include 
2PI diagrams beyond those of sunset topology in cf. Eq. 1)3. 20() and Fig. 13.11 As a first step 
it would be interesting to include the vertex corrections into the gap equation. 

The general effective action 1)2. 46() derived in this work puts the known effective theories 
for hot and dense matter on a more formal basis and gives better insight and control over 
the different approximations made. It connects previous results and puts them on a common 
footing. Thereby, it has the potential to provide a framework for systematic studies of various 
problems in the field of hot and dense systems. Besides an improvement of the result for the 
color-superconducting gap parameter beyond subleading order, I believe that it can serve as 
a convenient starting point to investigate other interesting problems pertaining to hot and/or 
dense quark matter. 



Appendix A 

Matsubara sums in quark loops 



In the following the Matsubara sum and the r— integrals in quark loops with three and four 
external gluon fields is performed, always starting from the general Eq. 



Three external glouns 

For three external glouns (A^ = 3) one has m = 1,2 and Eq. ()2.84|) becomes 



(27r) 



m=l 
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i=l 
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(2^ 
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6i ^^i'^i 
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6 J, Ti 



Using Eqs. (|T7^ and ((T7^ one finds 
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which leads to 
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(A.3) 



/3-Ti 

After having performed the integral over T2 one can combine three of the four remaining 
Ti -integrals by using fXjtl - fl^f-ss " f-s2 = -f-ss ^^'^ obtains 
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In the limit ^ — > and after using Eq. (|2.94|) this reproduces the corresoponding formula (A. 17) 



m 



13 



Four external Gluons 

In the case of four external gluons one has m = 1, 2, 3 and Eq. (|2.84j) becomes 

m=l i=l 



r rfili 3 3 ^ 
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Using Eqs. (|2.73p and l\2.T2\i and respecting < < /3 one finds 



-n4 + e-''''^<^0(Ti + T2-/3) + e E^J-/5 0(/3-n-r2) 



Hence, we have the fohowing r-integrals 



(A.6) 
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(A.7) 



/3-ri 2/3-Ti-r2 



Performing the r-integrals one by one and using Eqs. (|2.73|) and (|2.72j) after each integration to 
combine the corresponding terms one finally obtains the result 

d^k 
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J (2vr)3 - 
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P2-P3- ^262 + sses 



P1-P2- ^1^1 + •52e2 



P1-P3- siei + 5363 



Pi-pl- siei + 5464 
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In the limit ;U — > and after using Eq. (|2.94)) this reproduces the corresoponding formula (A. 19) 
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